The effect of stress concentration on the bending behavior of a ZrCuNiAl bulk metallic glass

The effect of stress concentration on the bending behavior of a ZrCuNiAl bulk metallic glass

Journal of Non-Crystalline Solids xxx (xxxx) xxx–xxx Contents lists available at ScienceDirect Journal of Non-Crystalline Solids journal homepage: w...

3MB Sizes 0 Downloads 74 Views

Journal of Non-Crystalline Solids xxx (xxxx) xxx–xxx

Contents lists available at ScienceDirect

Journal of Non-Crystalline Solids journal homepage: www.elsevier.com/locate/jnoncrysol

The effect of stress concentration on the bending behavior of a ZrCuNiAl bulk metallic glass Zhiliang Ninga,b,⁎, Yongjiang Huanga,b,⁎, Zhe Shenc, Haicao Suna,b, Weizhong Liangd, Jianfei Suna,b a

School of Materials Science and Engineering, Harbin Institute of Technology, Harbin 150001, PR China National Key Laboratory of Precision Hot Processing of Metals, Harbin Institute of Technology, Harbin 150001, PR China c Shenyang Engine Research Institute, Aero Engine Corporation of China, Shenyang 110000, PR China d School of Materials Science and Engineering, Heilongjiang University of Science and Technology, Harbin 150027, PR China b

A R T I C L E I N F O

A B S T R A C T

Keywords: Bulk metallic glasses Stress concentration Bending Fracture

The bending behaviors of a Zr55Cu30Ni5Al10 bulk metallic glass (BMG) samples with different notch radii are studied. The artificial notch will introduce stress concentration within the BMG sample during bending. That means, smaller notch radius leads to greater stress concentration around the notch area, and thus results in earlier failure and smoother fracture surface. Then, the sample with smaller notch radius possesses smaller fracture toughness. The relationship between the stress concentration and deformation feature has been discussed in terms of grease model.

1. Introduction

sample geometry dependence of bending ductility of BMG plates; thin plate usually exhibits ductility while thicker plate with the same chemical composition shows brittle feature under identical loading condition. Liu and Wang [28] demonstrated the shear band evolution under bending in a super-plastic Zr-based BMG. During service, BMG components are often suffered to a multi-axial complex stress state induced by artificial defects, such as holes, and notches [29–32]. As a BMG material with an artificial defect is loaded, a great stress concentration usually takes place near the defect. As a result, the stress distribution is inhomogeneous over the whole cross-section, which could exert a crucial role on the deformation of the BMG materials. Until now few research has been focused on the bending behaviors of BMGs after the introduce of stress concentration. In the present work, Zr55Cu30Ni5Al10 (at. %) alloy was selected as the model material due to its excellent glass forming ability. Notches with different radii, as stress concentrators, were introduced into the BMG samples. The effect of stress concentration induced by notch on the bending behaviors of the Zr-based BMG was studied. The fracture morphologies of the failed samples were examined after bending deformation. It is expected that the results obtained can provide a better understanding the bending deformation mechanism of BMGs, and thus extend their applications in complex multi-axial stress states.

During the last decades, bulk metallic glasses (BMGs), also termed as bulk amorphous alloys, have triggered intensive attention from materials science researchers due to their importance in both fundamental science and industrial uses [1–5]. However, at a temperature far below the glass transition temperature, the plastic deformation of BMGs usually localizes within few narrow shear bands of ~20 nm thick, which causes catastrophic failure and a lack of plasticity [6–8]. Some monolithic metallic glasses exhibit great plastic deformation capacity upon compressive loading [9,10] whereas almost all BMGs show disappointed tensile plasticity [11]. Such a brittle feature hinders seriously any practical applications of BMGs. Therefore, to understand the intrinsic deformation mechanism and thus to improve the ductility of BMGs have recently become one of the most important issues in BMGs research community [12–16]. As a kind of high-performance materials for engineering applications, knowledge regarding how BMGs behave under bending conditions is critical [17–21]. Therefore, the bending behaviors of BMGs deserve further attention. When a BMG alloy sample subjected to bending condition, the inherent stress gradient induced by multi-axial stress would favor the formation of multiple shear bands. This has aroused considerable research interests [22–24]. Zhang et al. [25] studied the multiplication of shear bands aided by free volume under three-point bending of a BMG. Conner et al. [26,27] examined the



Corresponding authors. E-mail addresses: [email protected] (Z. Ning), [email protected] (Y. Huang).

http://dx.doi.org/10.1016/j.jnoncrysol.2017.04.005 Received 5 January 2017; Received in revised form 15 February 2017; Accepted 11 April 2017 0022-3093/ © 2017 Elsevier B.V. All rights reserved.

Please cite this article as: Ning, Z., Journal of Non-Crystalline Solids (2017), http://dx.doi.org/10.1016/j.jnoncrysol.2017.04.005

Journal of Non-Crystalline Solids xxx (xxxx) xxx–xxx

Z. Ning et al.

Table 1 The maximum load, F, the maximum displacement, Dmax, stress concentration factor, average dimple height, h, and fracture toughness of the Zr55Cu30Ni5Al10 bulk metallic glass samples with different notch radii, ρ. Notch radius, ρ (mm)

Maximum displacement, Dmax (mm)

Maximum load, F(N)

Stress concentration factor, Kt

Dimple height, h (μm)

Fracture toughness, Pa·m1/2

0.2 0.3 0.5 1.0

0.3712 0.4465 0.4869 0.5531

868.42 959.86 1099.81 1226.78

2.808 2.319 1.823 1.315

5.56 6.38 6.54 7.09

38.42 43.27 46.88 51.56

2. Experimental procedures Ingots with a nominal chemical composition of Zr55Cu30Ni5Al10 were produced by arc-melting a mixture of the starting pure Zr, Cu, Ni, and Al elements, each with 99.5 wt% purity or better, under a Tigettered argon atmosphere. Prior to the preparation of the master ingots, pure Ti was melted enough to absorb the oxygen within the chamber, which is harmful to glass formation. To ensure a chemical homogeneity, each ingot was remelted at least four times with electromagnetic stir. Then as-cast BMG samples with a 3 × 30 × 75 mm3 dimension were prepared by drop casting the remelted alloys into a copper mold. The amorphous nature of the ascast samples was confirmed by X-ray diffraction (XRD) with Cu Kα radiation. Three-point bending samples with a 3 × 5 × 30 mm3 dimension were electrical discharge machined from the as-cast BMG alloy samples, and the top and the bottom surfaces of the samples were then carefully grinded and polished. Notches with a radius of 0.2 mm, 0.3 mm, 0.5 mm, or 1.0 mm were cut at one edge in the middle of the samples. Notch length is 2.5 mm. The outer appearance of the bending samples is shown in Fig. 1a. Bending tests were conducted on a computer-controlled Instron 5500R type machine at a displacement rate of 0.3 mm/min. The span was set to be 24 mm. During bending tests, the load-deflection curves were recorded continuously. After testing, fracture surfaces of all failed samples were examined in a scanning electron microscope, SEM (FEI Helios Nanolab600i) or a 3-D laser scanning confocal microscope (Phenom Pro X). In addition, a finite-element method (FEM) with ABAQUS software was adopted to simulate the stress distributions of the studied Zr-based BMG samples with different notch radii upon bending condition.

Fig. 1. (a) The outer appearance (from top to bottom shows the sample with a notch radius of 0.2, 0.3, 0.5, 1.0 mm, respectively) and (b) XRD patterns of three-point bending Zr55Cu30Ni5Al10 (at. %) bulk metallic glass samples.

3. Results and discussion Fig. 1b shows the XRD pattern obtained from the as-cast quaternary ZrCuNiAl alloy sample. It can be seen from Fig. 1b that the pattern of cast samples consists of only a series of broad diffraction maxima, without any detectable sharp Bragg peaks corresponding to crystalline phase. This confirms that the as-cast plate-shaped sample is amorphous in nature. Fig. 2 presents the representative load-deflection curves obtained from the room temperature bending tests of the studied Zr-based BMG alloy samples with different notch radii. As can be seen, for the studied Zr-based BMG alloy, the bending samples exhibit a similar elastic feature prior to failure, i.e., without any macroscopic plastic displacement. This indicates a typical brittle deformation feature for the studied BMG during bending. The maximum load, F and the maximum displacement, Dmax is 868.42 N and 0.3712 mm, 959.86 N and Fig. 2. The bending load-reflection curves of the Zr55Cu30Ni5Al10 bulk metallic glass samples with different notch radii (Color online).

2

Journal of Non-Crystalline Solids xxx (xxxx) xxx–xxx

Z. Ning et al.

Fig. 3. The stress field distribution near the notch during three-point bending deformation of the Zr55Cu30Ni5Al10 bulk metallic glass samples with different notch radii: a) ρ = 0.2 mm; b) ρ = 0.3 mm; c) ρ = 0.5 mm; and d) ρ = 1.0 mm (Color online).

Fig. 4. SEM micrographs of the whole fracture surface in the Zr55Cu30Ni5Al10 bulk metallic glass samples with different notch radii after three-point bending deformation tests: a) ρ = 0.2 mm; b) ρ = 0.3 mm; c) ρ = 0.5 mm; and d) ρ = 1.0 mm.

were used. It can be seen from Fig. 3 that the stress distributions in the notched Zr-based BMG samples are highly complex due to the existence of artificial notches during bending tests. The area near the notch tip exhibits a maximum stress, whereas, the area away from the notch has smaller stress. This suggests that U-shaped notch causes high stress level around the notch, indicating a stress concentration. It can be also noticed that smaller notch radius causes a higher maximum stress near the notch tip, demonstrating a greater stress concentration. The stress concentration factor, Kt, is defined as the ratio of the

0.4465 mm, 1099.81 N and 0.4869 mm, and 1226.78 N and 0.5531 mm, for the bending sample with notch radius ρ of 0.2 mm, 0.3 mm, 0.5 mm, and 1.0 mm, respectively, as listed in Table 1. These above results suggest that the smaller notch radius will result in earlier brittle failure during bending test, considering that the loading conditions are identical for all Zr-based BMG alloy samples. The numerical simulation results of the stress distributions for Zrbased BMG samples under a load of 100 N are displayed in Fig. 3a–d, in which an elastic modulus of 94.8 GPa and a Poisson's ratio of 0.366

3

Journal of Non-Crystalline Solids xxx (xxxx) xxx–xxx

Z. Ning et al.

Fig. 5. SEM micrographs of the fracture surface near the notch area in the Zr55Cu30Ni5Al10 bulk metallic glass samples with different notch radii after three-point bending deformation tests: a) ρ = 0.2 mm; b) ρ = 0.3 mm; c) ρ = 0.5 mm; and d) ρ = 1.0 mm.

maximum stress, σmax, to the nominal stress that would exist in the materials if the stress distribution keeps uniform, σnom; that is [33],

Kt =

σmax σnom

the smooth region is a region main consisting of vein-like patterns, marked by II in Fig. 5, similar with those of typical compressive fractured surfaces for BMGs, caused by the subsequent catastrophic failure. Here, we call it as transition region. Such vein-like pattern in the fractured surface of BMGs indicate ductile fracture feature, which can be ascribed to a substantial rise in temperature associated with the conversion of the stored elastic energy to local heat during final deformation stage [35]. These vein-like patterns on the fractography also clearly reveal a pure shear fracture process of a BMG. The liquid droplets, demonstrating a local great temperature rise [36–38], can be also observed in this region, as seen in Fig. 6. Lewandowski and Greer [39] experimentally estimated the local temperature rise based on a fusible coating and suggested that there exists a remarkable temperature rise of a few thousand Kelvin for a few nanoseconds within a shear band. These vein-like pattern and liquid droplet features observed on the fractured surface of the failed samples indicate that the bending deformation process is characterized by the shear band nucleation and fast propagation, which cause a catastrophic fracture. Interestingly, numerous radiate cores with different diameters can be also detected in the transition region of the fractured surface. A typical core is highlighted in Fig. 6e. Zhang et al. [11] suggested that, during tensile test, the veins radiate from cores and propagate towards the outside, causing the formation of a radiating core feature in the fractured surface. Zhang et al. [11] attributed this core feature to normal tension stress in the initial stage of fracture. Therefore, the transition region experiences a complex stress state of shear stress and normal stress. During the final stage of bending deformation, rapid fracture region forms in the fractured surface, see Fig. 7. Numerous dense and fine dimples dominate the rapid fracture region of the fractured samples with different notch radii. As the notch radius increases, the rapid fracture region becomes more and more uneven. Fig. 8 shows 3-D laser scanning confocal microscope images of the fractured samples. Fig. 8 inset presents the corresponding surface

(1)

Upon the Mode I bending conditions, the maximum stress concentration occurs near the tip of the notches, where σnom can be calculated by,

σnom =

3FS 2B (w − a )2

(2)

where w is the width of the sample, a is the notch length, F is the maximum load, S is the span, and B is the thickness of the sample. Substituting the simulated σmax values and the results calculated based on Eq. (2) into Eq. (1) yields the corresponding stress concentration factor for the BMG samples with different notch radii. The calculated Kt values are also listed in Table 1. As seen, for the studied Zr-based BMG, Kt increases from 1.315 to 2.808 as the notch radius ρ decreases from 1.0 mm to 0.2 mm. Fractography is an excellent methodology to understand the deformation behaviors of a material. Fig. 4 shows the whole fractured surface morphologies of the failed samples. The whole fractured surface of the BMG samples after bending deformation can be roughly divided into three regions. A smooth region can be found near the artificial notch for the all bending samples, as seen in Fig. 5 marked by I. This smooth region was caused by shear sliding [34]. Here, we call it as crack propagation region. The width of the smooth region is equal to the critical shear offset, an efficient parameter directly indicating the stable shear deformation ability [34]. The shear offsets of the samples with notch radius ρ of 0.2 mm, 0.3 mm, 0.5 mm, and 1.0 mm, is ~ 130 μm [Fig.5a], 235 μm [Fig.5b], 420 μm [Fig.5c], and 480 μm [Fig.5d], respectively. This mean that the sliding along shear bands is also an important bending deformation mechanism of BMGs. Close to

4

Journal of Non-Crystalline Solids xxx (xxxx) xxx–xxx

Z. Ning et al.

Fig. 6. SEM micrographs of the rapid propagation region of the fracture surface in the Zr55Cu30Ni5Al10 bulk metallic glass samples with different notch radii after three-point bending deformation tests: a) ρ = 0.2 mm; b) ρ = 0.3 mm; c) ρ = 0.5 mm; d) ρ = 1.0 mm; and e) the radiating patterns and liquid droplets (Color online).

and has numerous fine dimples. For the bending samples with other radii, the fractured surfaces exhibit the same features with 0.2 mm one [Fig. 8b–8d]. It can be also seen from Fig. 8 that, as the notch radius increases, the roughness of all the three regions in the fractured surface gradually increases. To establish the relationship between the fracture toughness and the 3-D fractured surface of the bent BMG samples, the fracture process has been quantitatively interpreted based on a “grease” model which was proposed by Takayama and Maddin [40]. Based on this model, as the glass in the fracture process zone can be considered as a viscous fluid due to great local heating associated with plastic deformation in the vicinity of the crack tip, a fluid meniscus will form in the crack tip [41]. Local stress causes a negative pressure inside the fracture process zone. The balance between the negative pressure and the surface tension of the viscous fluid would lead to the formation of a curved fluid surface, called a fluid meniscus [41]. Thus, the crack propagation can be

profile across the fractured surface. The color evolution from red to blue means the profile evolution from high to low. In order to quantitatively analyze the effect of notch radius on the bending behaviors of the studied Zr-based BMG alloy samples, the line profiles across the three regions, i.e., the crack propagation region, transition region, and rapid fracture region, are taken from the fracture surfaces, and shown in the inset of Fig. 8. It can be seen from Fig. 8 that the fracture features of three regions, especially the vein height, are different with each other. Taking the sample with a notch radius of 0.2 mm for example [Fig. 8a], its rapid fracture region consists of periodic peaks and valleys, and there is a large drop between peaks and valleys. The whole rapid fracture region shows a very rough feature. Besides the valleys with a relatively shallower depth, the transition region possesses some meniscus dimple-like features, and is a little smoother than the crack propagation region. Crack propagation region is smoother than both the transition zone and the rapid fracture region,

5

Journal of Non-Crystalline Solids xxx (xxxx) xxx–xxx

Z. Ning et al.

Fig. 7. SEM micrographs of the final fracture region of the fracture surface in the Zr55Cu30Ni5Al10 bulk metallic glass samples with different notch radii after three-point bending deformation tests: a) ρ = 0.2 mm; b) ρ = 0.3 mm; c) ρ = 0.5 mm;and d) ρ = 1.0 mm.

assumed to be a fluid flowing in a channel with a height of H. The shear dilatation of the viscous fluid would cause in the formation of numerous voids in the plastic deformation zone. As the cracks further propagate, the voids in the plastic deformation zone were merged with the main crack. Finally, numerous dimples with a height of h form in the fractured surface. Therefore, the channel height, H can be approximately estimated to be twice the vein height, h. According to the fluid meniscus instability, the relationship between the crack tip opening displacement (CTOD) and the critical stress intensity factor, i.e., fracture toughness, can be established as follows [41],

the artificial notch. This reveals that enhanced notch radius could lead to an improved fracture toughness of the glassy solids, as shown in Table 1. Therefore, the BMG sample with a larger notch radius fails sluggishly, and exhibits a higher maximum displacement and load, than that with a smaller notch radius, as shown in Fig. 2. Higher fracture toughness will result in a wider crack propagation region [Fig. 5], and higher roughness of the fracture surface for the sample with a smaller notch radius [Fig. 7, Fig. 8 and Table 1].

4. Summary

Kc =

CTODmσy E

(3) The effect of stress concentration induced by notch in the sample center on the deformation behaviors of a ZrCuNiAl bulk metallic glass has been studied at room temperature under bending loading. The notch radius of the bending sample exerts a crucial role on the deformation behaviors of the studied BMG. As the notch radius increases, both the maximum load and the maximum displacement prior to final failure increase. The fractured surface can be roughly divided into three regions, i.e., the crack propagation region, transition region, and rapid fracture region. The increase in the notch radius causes the gradual increase in the roughness of all the three regions in the fractured surface. Based on “grease” model, the relationship between the fracture toughness and the 3-D fracture surface of the bent BMG samples is quantitatively established. As the notch radius increases, the Kc value increases as well. Higher notch radius decreases the opening stress levels near the notch area, and promotes the plastic strain accumulation ahead of the notch, thereby causing higher fracture toughness of the studied BMG sample.

where m is a dimensionless constant depending on the material properties and the stress states, σyis the yield stress, and E is Young's modulus. The average vein height, h, is measured to be 5.56, 6.38, 6.54, and 7.09 μm for the bending sample with a notch radius of 0.2, 0.3, 0.5, and 1.0 mm, respectively, yielding Kc value of 38.42, 43.27, 46.88, and 51.56 MPa·m1/2 after substituting them into Eq. (3), as summarized in Table 1. Clearly, as the notch radius of the bending sample increases, the Kc value increases as well. From the above results, for the bending deformation of a BMG, a larger notch radius will result in a sluggish fracture, and higher fracture toughness. Next, this phenomenon will be interpreted based on the stress field around the notch. An increase in the notch radius leads to an increase in the plastic zone size ahead of the notch and causes this zone to rotate forward. Higher notch radius also results in a dramatic decrease in opening stress ahead of the artificial notch. Hydrostatic stress at all angles around the notch would reduce. Meanwhile, it also gives rise to a great increase in the plastic strain accumulation ahead of

6

Journal of Non-Crystalline Solids xxx (xxxx) xxx–xxx

Z. Ning et al.

Fig. 8. Three dimensional (3-D) laser scanning confocal microscope images of the fracture surfaces in the Zr55Cu30Ni5Al10 bulk metallic glass samples with different notch radii after three-point bending deformation tests: a) ρ = 0.2 mm; b) ρ = 0.3 mm; c) ρ = 0.5 mm; and d) ρ = 1.0 mm; The inset shows the corresponding surface profile. The column from left to right represents the rapid fracture region, transition region, and crack propagation region, respectively (Color online).

7

Journal of Non-Crystalline Solids xxx (xxxx) xxx–xxx

Z. Ning et al.

620–625. [19] C.C. Yu, J.P. Chu, C.M. Lee, W. Diyatmika, M.H. Chang, J.Y. Jeng, Y. Yokoyama, Bending property enhancements of Zr55Cu30Al10Ni5 bulk metallic glass: effects of various surface modifications, Mater. Sci. Eng. A 633 (2015) 69–75. [20] Y.J. Huang, D.J. Wang, D.P. Wang, Z.X. Zhang, J. Shen, Bending behavior of TiZrNiCuBe bulk metallic glass, J. Alloys Compd. 541 (2012) 359–364. [21] F. Haag, D. Beitelschmidt, J. Eckert, K. Durst, Influences of residual stresses on the serrated flow in bulk metallic glass under elastostatic four-point bending – a nanoindentation and atomic force microscopy study, Acta Mater. 70 (2014) 188–197. [22] F. Jiang, H.F. Wang, M.Q. Jiang, G. Li, Y.L. Zhao, L. He, J. Sun, Ambient temperature embrittlement of a Zr-based bulk metallic glass, Mater. Sci. Eng. A 549 (2012) 14–19. [23] C.X. Xie, Y.Z. Yang, S.Y. Zhong, S. Li, S.C. Deng, Formation, magnetic properties and bending deformation of Fe-based amorphous alloy without metalloids, J. Alloys Compd. 695 (2017) 877–880. [24] L.C. Zhang, F. Jiang, Y.L. Zhao, S.B. Pan, L. He, J. Sun, Shear band multiplication aided by free volume underthree-point bending, J. Mater. Res. 25 (2010) 283–291. [25] R.T. Qu, H.S. Liu, Z.F. Zhang, In situ observation of bending stress–deflection response of metallic glass, Mater. Sci. Eng. A 582 (2013) 155–161. [26] R.D. Conner, Yi Li, W.D. Nix, W.L. Johnson, Shear band spacing under bending of Zr-based metallic glass plates, Acta Mater. 52 (2004) 2429–2434. [27] R.D. Conner, W.L. Johnson, N.E. Paton, W.D. Nix, Shear bands and cracking of metallic glass plates in bending, J. Appl. Phys. 94 (2003) 904–911. [28] Y.H. Liu, W.H. Wang, Shear bands evolution in bulk metallic glass with extended plasticity, J. Non-Cryst. Solids 354 (2008) 5570–5572. [29] W.D. Li, H.B. Bei, Y.F. Gao, Effects of geometric factors and shear band patterns on notch sensitivity in bulk metallic glasses, Intermetallics 79 (2016) 12–19. [30] I. Singh, R. Narasimhan, Notch sensitivity in nanoscale metallic glass specimens: insights from continuum simulations, J. Mech. Phys. Solids. 86 (2016) 53–69. [31] J.X. Zhao, F.F. Wu, R.T. Qu, S.X. Li, Z.F. Zhang, Plastic deformability of metallic glass by artificial macroscopic notches, Acta Mater. 58 (2010) 5420–5432. [32] F. Gong, S.H. Chen, J.Q. Ran, Z. Yang, J. Ma, Tuning the performance of bulk metallic glasses by milling artificial holes, Mater. Sci. Eng. A 668 (2016) 50–54. [33] W.D. Pilkey, D.F. Pilkey, Peterson's Stress Concentration Factors, John Wiley & Sons, Inc, 2008. [34] F.F. Wu, W. Zheng, S.D. Wu, Z.F. Zhang, J. Shen, Deformation and fracture behaviors of Ti-based metallic glass under multiaxial stress state, Acta Mater. 60 (2012) 2073–2081. [35] C.T. Liu, L. Heatherly, D.S. Eaton, C.A. Carmichael, J.H. Schneibel, C.H. Chen, J.L. Wright, M.H. Yoo, J.A. Horton, A. Inoue, Test environment and mechanical properties of Zr-base bulk amorphous alloy, Metall. Mater. Trans. A 29 (1811) (1998) 1811–1820. [36] W.J. Wright, R.B. Schwarz, W.D. Nix, Localized heating during serrated plastic flow in bulk metallic glasses, Mater. Sci. Eng. A 319–321 (2001) 229–232. [37] G. Wang, Y.J. Huang, J. Shen, Novel TiCuNiCo composites with high fracture strength and plasticity, Mater. Design. 33 (2012) 226–230. [38] H. Zhai, H. Wang, F. Liu, Effects of Sn addition on mechanical properties of Ti-based bulk metallic glass composites, Mater. Design. 110 (2016) 782–789. [39] J.J. Lewandowski, A.L. Greer, Temperature rise at shear bands in metallic glasses, Nat. Mater. 5 (2006) 15–18. [40] S. Takayama, R. Maddin, Fracture of amorphous Ni-Pd-P alloys, Philos. Mag. 32 (1975) 457–470. [41] G. Wang, K.C. Chan, X.H. Xu, W.H. Wang, Instability of crack propagation in brittle bulk metallic glass, Acta Mater. 56 (2008) 5845–5860.

Acknowledgements This work was financially supported by the National Natural Science Foundation of China (NSFC) under Grant Nos. 51671067, 51371078 and 51671070. References [1] W. Klement, R.H. Willens, P. Duwez, Non-crystalline structure in solidified goldsilicon alloys, Nature 187 (1960) 869–870. [2] Y.J. Huang, J.C. Khong, T. Connolley, J. Mi, The onset of plasticity of a Zr-based bulk metallic glass, Int. J. Plast. 60 (2014) 87–100. [3] Y.J. Huang, H.B. Fan, D.J. Wang, Y. Sun, F.Y. Liu, J. Shen, J.F. Sun, J. Mi, The effect of cooling rate on the wear performance of a ZrCuAlAg bulk metallic glass, Mater. Des. 58 (2014) 284–289. [4] T.C. Hufnagel, C.A. Schuh, M.L. Falk, Deformation of metallic glasses: Recent developments in theory, simulations, and experiments, Acta Mater. 109 (2016) 375–393. [5] Y.J. Huang, X. Cheng, H.B. Fan, S.S. Guan, Z.L. Ning, J.F. Sun, Crystallization of a Ti-based bulk metallic glass induced by electropulsing treatment, J. Iron Steel Res. Int. 23 (2016) 69–73. [6] Y.J. Huang, W. Zheng, F.L. He, J. Shen, The temperature dependent dynamic mechanical response of a ZrCuNiAl bulk metallic glass, Mater. Sci. Eng. A 551 (2012) 100–103. [7] B.A. Sun, W.H. Wang, The fracture of bulk metallic glasses, Prog. Mater. Sci. 74 (2015) 211–307. [8] P. Thurnheer, F. Haag, J.F. Löffler, Time-resolved measurement of shear-band temperature during serrated flow in a Zr-based metallic glass, Acta Mater. 115 (2016) 468–474. [9] Y.J. Huang, J. Shen, J.F. Sun, Bulk metallic glasses: smaller is softer, Appl. Phys. Lett. 90 (2007) 081919. [10] C. Chen, J.L. Ren, G. Wang, K.A. Dahmen, P.K. Liaw, Scaling behavior and complexity of plastic deformation for a bulk metallic glass at cryogenic temperatures, Phys. Rev. E92 (2015) 012113. [11] Z.F. Zhang, J. Eckert, L. Schultz, Difference in compressive and tensile fracture mechanisms of Zr59Cu20Al10Ni8Ti3 bulk metallic glass, Acta Mater. 51 (2003) 1167–1179. [12] Y.J. Huang, J.C. Khong, T. Connolley, J. Mi, In situ study of the evolution of atomic strain of bulk metallic glass and its effects on shear band formation, Scr. Mater. 69 (2013) 207–210. [13] Z. Lu, W. Jiao, W.H. Wang, H.Y. Bai, Flow unit perspective on room temperature homogeneous plastic deformation in metallic glasses, Phys. Rev. Lett. 113 (2014) 045501. [14] J. Schroers, Bulk Metallic Glass. Phys. Today. 66 (2013) 32–37. [15] S.V. Ketov, Y.H. Sun, S. Nachum, Z. Lu, A. Checchi, A.R. Beraldin, H.Y. Bai, W.H. Wang, D.V. Louzguine-Luzgin, M.A. Carpenter, A.L. Greer, Rejuvenation of metallic glasses by non-affine thermal strain, Nature 524 (2015) 200–203. [16] R.T. Qu, Z.Q. Liu, G. Wang, Z.F. Zhang, Progressive shear band propagation in metallic glasses under compression, Acta Mater. 91 (2015) 19–33. [17] T. Yuan, G.Y. Wang, Q.M. Feng, P.K. Liaw, Y. Yokoyama, A. Inoue, Modeling size effects on fatigue life of a zirconium-based bulk metallic glass under bending, Acta Mater. 61 (2013) 273–279. [18] Y. Hu, H.H. Yan, J.F. Li, Y.H. Zhou, Bending plasticity of Zr55Al10Ni5Cu30 bulk metallic glass with monolithic amorphous structure, J. Alloys Compd. 688 (2016)

8