Elsevier

Mechanics of Materials

Volume 148, September 2020, 103494
Mechanics of Materials

Research paper
On the strain rate sensitivity of size-dependent plasticity in BCC iron at elevated temperatures: Discrete dislocation dynamics investigation

https://doi.org/10.1016/j.mechmat.2020.103494Get rights and content

Highlights

  • MDDP simulations were carried out on BCC iron samples to study the size effect under high temperature/ high strain rates.

  • Strong size effect is manifested irrespective of temperature.

  • Size effect is weakly dependent on the applied rate of deformation but is dependent on temperature.

Abstract

Multiscale discrete dislocation plasticity (MDDP) simulations are carried out to investigate the mechanical response and microstructure evolution of BCC iron micropillars under combined high temperature and strain rate deformation. The simulations are conducted at sizes ranging between 0.25 μm and 2 μm under an applied deformation rate ranging between 103s−1 and 107s−1 and subjected to different temperatures. MDDP based constitutive equation interrelating the size effect exponent to strain rate and temperature is also proposed indicating that the exponent is relatively sensitive to temperature and at a lesser degree to strain rate.

Detailed investigation of the microstructure shows that self-multiplication of dislocations is responsible for the strengthening mechanism in BCC iron micropillars. At low temperatures and small sizes, screw dislocations have a weak effect on plasticity for a certain period of time but subsequently control the self-multiplication process. At larger sizes, the motion of screw dislocations is responsible for plasticity at low temperatures. Due to the large volume size, screw dislocations are entangled inside the sample leading to a self-multiplication of dislocations via cross slip and other dislocation–dislocation interactions. At higher temperatures and for all sample sizes, mixed dislocations control plasticity via the multiplication of a complex network of dislocations. MDDP generated results are in good agreement with previous experimental studies on BCC metals.

Introduction

Over the last decade, unraveling the mechanical properties and deformation mechanisms of micro and nano pillars has been a major concern for scientists and researchers. Unlike bulk materials, the yield strength of micro and submicron pillars increases as the sample size decreases. However, this “smaller is stronger” trend differs in an essential way between materials. At the micro and nano scales, the size effect is quantified by an empirical power-law that relates the yield strength to the size of the pillar such that:σy=kdnwhere σy is the yield strength, d is the sample size, k and n are material constants determined experimentally. The power dependence expressed in Eq. (1) is highly correlated with the crystal structure of the material. Experimental (Greer et al., 2008; Greer and Nix, 2006; Jennings et al., 2011; Kim et al., 2010; Rogne and Thaulow, 2015; Schneider et al., 2013, 2009) and simulation studies (Greer et al., 2008; Kim et al., 2009) show that, unlike FCC metals where size effect is significant with the exponent n ranging between 0.6 and 1 (Frick et al., 2008; Greer and Nix, 2006; Kim and Greer, 2009; Volkert and Lilleodden, 2006), size effect in BCC metals is less pronounced at temperatures below the critical temperature with an exponent ranging between 0.2 and 0.5. The critical temperature is the temperature above which the movement of dislocations becomes athermal, i.e., screw and edge dislocations have the same mobility and move by phonon drag. As the temperature exceeds the critical temperature, size effect is more pronounced and the behavior is similar to that of FCC metals (Greer et al., 2008; Rogne and Thaulow, 2015; Schneider et al., 2013, 2009; Torrents et al., 2016). This behavior is attributed to the high lattice friction of screw dislocations which is size independent but temperature dependent. Lattice friction is the resolved shear stress required to move a dislocation at a temperature greater than 0 K. In BCC metals, the complex core structure of screw dislocations (Domain and Monnet, 2005; Gilbert et al., 2011) is responsible for the high lattice friction as well as the difference in mobility between screw and edge dislocations. Below the critical temperature, screw dislocations have high lattice friction, thus, they move by thermal activation at a mobility less than their edge counterparts.

Experimental studies on different types of BCC pillars emphasize the dependency of the size exponent on the critical temperature. Schneider et al. (2009) performed a series of compression experiments at room temperature on four BCC pillars and found that the size exponent is inversely proportional to the respective critical temperature of the micropillars. The higher the critical temperature, the less pronounced the size effect. Kim et al. (2010) performed uniaxial compression and tension experiments and found that the yield strength and flow stresses for the metals that have the lowest critical temperature are strongly size-dependent. Besides the high dependency on critical temperature, Torrents et al. (2016) showed that size effect is also affected by dislocation density. Atomistic simulations on Tungsten, performed by Cui et al. (2016) revealed that even at temperatures below the critical temperature, the mobility of screw dislocations becomes dominated by phonon drag as the sample size decreases below 200 nm. Hence, mixed dislocations dominate plasticity as screw dislocations are unable to be stored in the limited sample volume.

Apart from the effect of the critical temperature, the power exponent is strain rate affected (Huang et al., 2015; Schneider et al., 2013; Wei et al., 2004). Schneider et al. (2013) investigated the effect of strain rate sensitivity on submicron Molybdenum pillars and found that these micropillars resemble their bulk counterparts. Huang et al. (2015) reported that the strain rate sensitivity decreases as the sample size is decreased and becomes negligible as the pillar is reduced below 200 nm. They attributed this behavior to the competing effect between lattice friction and bow out stresses. As the pillar size is reduced to a few nanometers, bow out stresses dominate the deformation mechanisms because few dislocations are present, leading to a decrease in strain rate sensitivity.

Even though the size exponent of BCC pillars is similar to FCC pillars at temperatures greater than the critical temperature, the deformation mechanisms in FCC and BCC pillars are not the same. Extensive studies on FCC pillars proposed two models to account for the size effect at the microstructure level. One model is the dislocation starvation model followed by surface nucleation. Greer and Nix (2006) were the first to elaborate the dislocation starvation model, where they observed that the dislocation density is less than the initial dislocation density when the sample is plastically deformed. This was attributed to the annihilation of dislocations at the free surfaces. After annihilation and to ensure plasticity, dislocations are generated and multiplied via surface nucleation (Greer et al., 2011; Jennings et al., 2011; Zhu and Li, 2010). The other model is the single arm source (SAS) model (Parthasarathy et al., 2007; Rao et al., 2008). In this model, dislocations are multiplied when single arm sources are truncated at the surface. On the contrary to FCC structure, these models are not found in BCC even at high temperatures where BCC and FCC size effects are similar. In BCC metals, dislocation multiplication is detected during plastic deformation (Weinberger and Cai, 2008). At temperatures below the critical temperature, screw dislocations are sessile and don't annihilate at the free surface, therefore they have enough time to interact between each other, cross slip and therefore multiply. At temperatures greater than the critical temperature, previous atomistic simulations (Greer et al., 2008; Weinberger and Cai, 2008) revealed that dislocations form segments inside the pillars interacting with each other's and with the free surface leading to dislocations multiplication.

Despite all the studies that have been developed in an effort to understand the deformation mechanisms and mechanical properties of the size effect in BCC metals, there are still some unresolved issues related to the effect of high temperature/ high strain rate on the size exponent and the dynamic evolution of dislocations in BCC iron. In order to get a comprehensive analysis of the combining effects of high temperature and high strain rate on size effect in iron single crystals, multiscale dislocation dynamics plasticity (MDDP) simulations have been utilized. The simulations are conducted to mimic loading conditions at temperatures ranging between 300 K and 900 K and strain rates ranging between 103s−1 and 107s−1 at different pillar sizes ranging between 0.25 μm and 2 μm under monotonic compressive loading conditions.

This paper is organized as follows. In Section 2, the MDDP approach, the simulation set up and the data used for the simulations are described. In Section 3, the results are discussed and a dislocation dynamics constitutive equation, based on MDDP, is established correlating the exponent with temperature and strain rate under high rate of deformation over a wide range of temperatures.

Section snippets

MDDP

MDDP (Zbib et al., 2005; Kattoura and Shehadeh, 2014, Shehadeh, 2012, Zbib and Diaz de la Rubia, 2002) is a multiscale simulation model that merges the three dimensional discrete dislocation dynamics (DDD) with finite element (FE) analysis resulting in a hybrid elasto-viscoplastic model.

At the macro level, the material is assumed to obey the basic laws of continuum mechanics where the linear momentum balance and the energy balance equations are applied:divS=ρv˙pρCvT=0.9Sε˙pwhere S is the Cauchy

Plastic deformation response and microstructure evolution

As presented in Section 2.1, it is expected that the activities and interactions of dislocations are influenced by the imposed strain rate and temperature due to their effects on the dislocation velocity, mobility and lattice friction. Additionally, in small volumes as in micropillars, the limited space sets of samples of 0.25 μm and 1 μm size subjected to a same strain rate of 105 s  1 and a deformation temperature of 300 K and 600 K. Fig. 2a-f depicts snapshots of the dislocations

Conclusion

We have carried out MDDP simulations to investigate the combined effect of size, high temperature and high strain rates on the mechanical response and microstructure evolution of BCC iron. The simulation results show good agreement with previous experimental studies performed on BCC metals. From the above results and discussion, the following can be deduced:

  • For smaller sizes and low temperatures, mixed dislocations control plasticity at the onset of yielding and up to 2% strain. At high strain,

Declaration of Competing Interest

None.

Acknowledgements

Support from the Research Board at the American University of Beirut is greatly acknowledged.

References (40)

Cited by (15)

  • On the investigation of size effect in bcc α-iron under high strain rate and high temperature: Multiscale dislocations dynamics simulations

    2023, European Journal of Mechanics, A/Solids
    Citation Excerpt :

    Simulation results are discussed and analyzed in Section 3. Previous MDDP simulations (El Ters and Shehadeh, 2019, 2020) investigated the mechanical response and microstructure evolution of bulk BCC iron as well as BCC iron micropillars subjected to high strain rate and high temperature. Here we carried out MDDP simulations to mimic the mechanical response of micro-polycrystalline BCC iron.

  • Assessment of the time-dependent behavior of dislocation multiplication under shock loading

    2022, International Journal of Plasticity
    Citation Excerpt :

    However, the collective behavior of dislocation assemblies is far from well understood, especially that at high strain rate, because present diagnostic techniques still cannot observe the collective behaviors of dislocation assemblies with high time resolution (Milathianaki et al., 2013; Sliwa et al., 2018; Turneaure et al., 2018) and first-principles simulation techniques are incapable of describing the collective behaviors of dislocation assemblies at larger time or spatial scales (Zepeda-Ruiz et al., 2017, 2021). Molecular dynamics (MD) (Gilbert et al., 2011; Queyreau et al., 2011; Branicio et al., 2013; Li et al., 2015; Prieto-Depedro et al., 2015; Xu et al., 2015; Xiong et al., 2015; Wang et al., 2014, 2015, 2017; Zepeda-Ruiz et al., 2017, 2021; Krasnikov and Mayer, 2015, 2017; Mayer et al., 2021; Zhu et al., 2015, 2017; Kuksin and Yanilkin, 2013; Yanilkin et al., 2014; Xiong et al., 2020; Fan et al., 2021; Li et al., 2021) and dislocation dynamics (Huang and Li, 2015; Huang et al., 2020; Leung et al., 2016; Zhang and Ngan, 2018; Agarwal et al., 2020; Shehadeh and Zbib, 2016; Ters and Shehadeh, 2019, 2020; Monavari and Zaiser, 2018; Gurrutxaga-Lerma et al., 2015a; Cui et al., 2015, 2016a, 2016b, 2018, 2019) are the most effective numerical methods deployed to investigate dislocation evolution. Given that the spatial and time scales accessible by MD remain limited, dislocation dynamics, including discrete dislocation dynamics (DDD) and continuum dislocation dynamics (Monavari and Zaiser, 2018), are more preferred to study the collective behavior of dislocation assemblies.

View all citing articles on Scopus
View full text