Acta Polytechnica CTU Proceedings doi:10.14311/APP.2018.15.0069 Acta Polytechnica CTU Proceedings 15:69–73, 2018 © Czech Technical University in Prague, 2018 available online at http://ojs.cvut.cz/ojs/index.php/app MODELING OF MONOCRYSTALLINE MAGNESIUM MICROBEAM BENDING Jiří Němečeka, ∗, Jan Maňákb, Jiří Němečeka a Czech Technical University in Prague, Faculty of Civil Engineering, Department of Mechanics, Thákurova 7, 166 29 Prague 6, Czech Republic b Institute of Physics of ASCR, Na Slovance 1999/2, 182 21 Prague 8, Czech Republic ∗ corresponding author: jiri.nemecek.1@fsv.cvut.cz (Ph.D. candidate) Abstract. This paper presents a numerical simulation of a micro-scale experiment on a magnesium alloy. Micro cantilever beams were fabricated using Focused Ion Beam technology in a single crystal of Mg. The cantilever beams have dimensions in the order of a few micrometers and a pentagonal cross section. Nanoindenter was used for cantilever beam bending and load-displacement curve was received. Cantilevers with two different crystallographic orientations were chosen for the experiment. Three dimensional numerical FE model with elastoplastic behavior respecting crystal anisotropy was used to fit experimental load displacement curves. Strengths and deformation energy were evaluated from the models for each cantilever. Keywords: Magnesium, microbeams, crystal anisotropy, FE model. 1. Introduction Magnesium crystallizes in a hexagonal close packed structure (HPC) with various deformation mecha- nisms such as basal plane slipping {0001} <1120>, prismatic slip {1010} <1120> and pyramidal slip {1011} <1123> or twinning {1012} <1011> [1]. The deformations start from a defect or on boundaries between different crystalline grains. In the first sce- nario, slipping happens when the crystal planes slide along each other and in the second, twinning is a concentrated shift of atomic positions [2]. Magnesium alloys fulfill the demand for a lightweight construction material with good recycling potential. Nowadays magnesium is commonly used in automotive, train and aircraft industry for its high tensile strength, low strength to weight ratio and ex- cellent corrosion resistance. However, magnesium also has negative properties such as tendency for crack- ing and brittle failure which is probably given due to intercrystalline failure and local transcrystalline fracture at twin zone [3, 4]. Nanoindentation is a powerful tool which uses a small diamond tip pressed into the material to measure mechanical properties such as modulus of elasticity, hardness and creep from load-displacement curves [5]. Nanoindentation usually operates at microscale and is a perfect tool for measuring properties of an individ- ual crystal. Besides classical nanoindentaion testing there are possibilities to use nanoindenter tip for bend- ing microscale samples. Micrometer samples can be fabricated using Focused Ion Beam (FIB) technology. FIB instrument is very similar to a scanning electron microscope (SEM) but instead of electrons it uses Gal- lium ions to hit the material and cause its milling [6]. FIB energy is typically between 10 and 50 kV, with beam currents varying between 1 pA and 10 nA [7] and is capable of creating samples with varied ge- ometry, mostly pillars or triangular and pentagonal micro beams. Micro-beam bending is advantageous for evaluating tensile strength and fracture energy. In this work, micro cantilevers beams were tested in specific crystal directions where twinning occurs. During the testing, load-displacement curves were measured. Also, a simple numerical model with good reproducibility of the curves gained on different crys- talline orientations has been developed. Figure 1. Crystal orientation with respect to loading direction. 2. Experimental part A sample from pure magnesium was prepared. The sample was cut by a precise diamond saw and finely polished to achieve average roughness in the range of a few nanometers. EBSD (Electron backscatter diffraction) was performed on the polished surface to find its crystal orientation. Micro beams were fabricated on the identified places by the FEI Quanta 3D FEG dual beam instrument combining SEM and FIB. All samples were micro cantilever beams with length of about 22 µm and have a pentagonal cross section with a width of about 3 µm and variable height 69 http://dx.doi.org/10.14311/APP.2018.15.0069 http://ojs.cvut.cz/ojs/index.php/app J. Němeček, J. Maňák, J. Němeček Acta Polytechnica CTU Proceedings Figure 2. SEM images of fabricated micro beam with C crystal orientation: side view (top) and front view (bottom). from 3.5 to 5 µm. Two types of cantilevers with two prismatic orientations, labelled C and D, were created, see Figure 1. SEM images of the fabricated cantilevers are shown in Figures 2 and 3. Circle marks were created on each cantilever approx- imately 2 µm away from the free end. At this mark, nanoindentation was performed with a Berkovich tip to bend the cantilever. Standard displacement con- trolled tests were performed with a constant increment of displacement up to a maximal value of 3000 nm. In total, four cantilevers were bent for each orientation, however, the curves had the same trend. In this paper were picked only the typical curves. Typical experi- mental load-displacement curves measured are shown in Figure 5. In both cases firstly, elastic deformation characterized by reversible response of the lattice to the displacement. Secondly, plastic deformation of the crystal, caused mostly by twinning for the cases Figure 3. SEM images of fabricated micro beam with D crystal orientation: side view (top) and front view (bottom). of magnesium crystal orientations, takes place. After reaching a critical strain, twinning causes a release of energy that is demonstrated as a radical drop of load in Figure 5. 3. Numerical modelling A numerical model was created in ANSYS Workbench software 17.2. A three dimensional model was pro- posed due to the irregular cross section of the can- tilever micro beam and due to previous experience of authors with micro beam modelling [10]. Material constants for elastic anisotropy of pure Mg were taken from the literature [8, 9] and are given in Table 1. For the hexagonal crystal, five independent elastic constants are needed and a sixth is calculated from 70 vol. 15/2018 Modeling of Monocrystalline Magnesium Microbeam Bending Figure 4. Directional deformation (Z axis) at maximal prescribed displacement 3 µm of the cantilever labelled as C with shown mesh. C11 [GPa] C12 [GPa] C13 [GPa] C33 [GPa] C44 [GPa] C66 [GPa] Basal 58.58 25.02 20.79 61.11 16.58 16.78 Prismatic 61.11 20.79 25.02 58.58 16.58 20.16 Table 1. Elastic constant of pure Mg for basal and prismatic crystal orientation [8, 9]. 0 1000 2000 3000 0 100 200 300 400 500 Displacement [nm] L oa d [µ N ] C - Experimental D - Experimental Figure 5. Typical experimental load-displacement curves of two different crystal orientations measured during cantilever bending. Plastic strain [µm µm−1] Stress [MPa] C D 0 392 260 0.07 470 410 0.010 495 420 0.015 530 425 0.020 545 430 Table 2. Plasticity data for the two Prismatic orien- tations of crystal C and D. the equation C66 = 1 2 (C11 − C12). (1) The constants are given for basal crystal orientation only and need to be transformed for prismatic orien- tation. Plasticity of the beams was modelled with multilin- ear isotropic hardening. The Yield stress was directly calculated from the bending experiment, where the values for cantilevers C and D are 392 MPa and 260 MPa respectively. The data of plastic strain and stress are shown in Table 2. The cantilever is filled with 3-D SOLID186 elements which are defined by 20 nodes and have three degrees of freedom per node, translation in the nodal x, y, and z directions. This el- ement supports von Mises plasticity and is suitable for large deflection [11]. The boundary conditions were set similar to the real cantilever with fixed support on one end and a displacement increment of up to 3000 nm. The deformed shape of the cantilever model is shown in Figure 4. The deformation process of cantilever bending can be characterized with deformation energy. A part of the total deformation energy that occurs before the first twinning was calculated from the stress-strain curves as follows G = ∫ �crit 0 σd�, (2) where σ is the axial stress and �crit is the critical strain when the first twinning happens. 71 J. Němeček, J. Maňák, J. Němeček Acta Polytechnica CTU Proceedings 0 1000 2000 3000 0 100 200 300 400 500 Twinning Displacement [nm] L oa d [µ N ] C - Experimental C - Fit Figure 6. Experimental and calculated load- displacement curves of C prismatic crystal orientation with marked the first twinning. 4. Results and discussion Calculated load-displacement curves can be seen in Figures 6 and 7. In both cases fitting curves exhibit the same behavior as in the experiment on the elas- tic part that matches perfectly. The curves are also with good agreement after Yield point until the first twining which is observed on the experimental curves when a radical drop of load occurs. Unfortunately this drop brings very difficult task for modelling and it was not considered in this work. According to re- search on austenitic high managanese steels showing twinning-included plasticity by C.Y. Sun et al. [12], it would be necessary to employ crystal plasticity in- corporating the twinning mechanism, where eighteen material constants are needed to be obtained experi- mentally. The first twinning, labelled as wcrit, occurs at displacement of 1.2 µm for C and 2 µm for D orien- tation. The shape of the fitting curves after twinning has the same trend and would match well with the experimental curve if no twinning had happened with only one exception. The slope of the un-loading curve for D orientation is slightly different. The axial stress in the X axis direction is shown in Figure 9. The strain was calculated as the equivalent total strain composed from elastic strain and equiv- alent plastic stain computed of von Mises criterion. Stress-strain curves calculated at the top fixed end of the cantilever are shown in Figure 8. However, this diagram is again partly correct till the wcrit value; the data after this point are not relevant. Thus, crit- ical strain (�crit) and stress were deducted from the model. The normal stress till the first drop of load is 619.9 MPa with a strain of 0.0323 for the C can- tilever and 592.9 MPa with value of strain 0.0548. The stress-strain curves made it possible to compare the deformation behavior of micro-beams with dif- ferent orientations. Based on this comparison, the 0 1000 2000 3000 0 100 200 300 400 500 Twinning Displacement [nm] L oa d [µ N ] D - Experimental D - Fit Figure 7. Experimental and calculated load- displacement curves of D prismatic crystal orientation with marked the first twinning. 0.000 0.032 0.055 0.080 0.120 0 100 200 300 400 500 592 620 700 800 900 Strain [-] St re ss [M P a] C D Figure 8. Calculated stress-strain curves of C and D prismatic crystal orientations with labelled critical stresses and strains. Yield point is higher for C orientation while the stored energy before the first twinning is higher for D. The energies calculated according to Equation 2 for C and D are 16.09 Jm−3 and 28.08 Jm−3 , respec- tively. That makes D orientation more ductile. 72 vol. 15/2018 Modeling of Monocrystalline Magnesium Microbeam Bending Figure 9. Deformed shape and axial stress (X axis) at critical displacement of 2 µm with uninformed wireframe of D cantilever beam. 5. Conclusions This paper shows numerical modelling of micro can- tilever beams made from pure magnesium with two prismatic orientations of the crystals. The elastic part, modelled with five independent anisotropic constants, gave a perfect match to the experiments. Inelastic deformations were captured with isotropic plasticity with multilinear hardening law until the first twinning. The twinning mechanism was not captured and left for future work. It was found that the fitting curves have the same trend as experimental curves if the twinning was removed, which suggests that the model is correct and only the twinning mechanism needs to be incorporated using e.g.crystal plasticity by Sun et al. [12]. Critical displacements wcrit, corresponding to the critical strains and stresses at the time of the first twinning were evaluated. It was found that the critical stress for C orientation is 30 MPa higher compared to D. However, the critical strain is higher for D. Thus, the deformation energy stored in the material before the first twinning is higher for D. List of symbols wcrit Displacement at first twinning [nm] �crit Strain at first twinning [–] σ Axial stress (X axis) [MPa] G Deformation energy [J m−3] C66 Elastic constant FIB Focused ion beam SEM Scanning electron microscope EBSD Electron backscatter diffraction Acknowledgements Financial support of the Czech Science Foundation (projects 17-18652S) is gratefully acknowledged. References [1] K. Fukuda, Y. Koyanagi, M. Tsushida, et al. Activation stress for slip systems of pure magnesium single crystals in pure shear test. Materials transactions 58(4):587–591, 2017. [2] N. K. Bourne. Unexpected twins. Physics 9:19, 2016. [3] F. Kaiser, K. Kainer. Magnesium alloys and technology. John Wiley & Sons, 2003. [4] T. Homma, N. Kunito, S. Kamado. Fabrication of extraordinary high-strength magnesium alloy by hot extrusion. Scripta Materialia 61(6):644–647, 2009. [5] A. C. Fischer-Cripps. Introduction to contact mechanics, vol. 87. Springer, 2000. [6] F. A. Stevie, L. Giannuzzi, B. Prenitzer. The focused ion beam instrument. In Introduction to Focused Ion Beams, pp. 1–12. Springer, 2005. [7] S. Reyntjens, R. Puers. A review of focused ion beam applications in microsystem technology. Journal of micromechanics and microengineering 11(4):287, 2001. [8] R. Hearmon. The elastic constants of anisotropic materials. Reviews of Modern Physics 18(3):409, 1946. [9] M. A. Kumar, I. J. Beyerlein, R. J. McCabe, C. N. Tome. Grain neighbour effects on twin transmission in hexagonal close-packed materials. Nature communications 7:13826, 2016. [10] J. Nemecek, J. Nemecek, J. Manak. Fracture properties of cement studied by nanoindentation and fib In Press 2018. [11] ANSYS, Inc. Ansys reference manual version 17.2, 2017. https://www.ansys.com/. [12] C. Sun, N. Guo, M. Fu, S. Wang. Modeling of slip, twinning and transformation induced plastic deformation for twip steel based on crystal plasticity. International journal of plasticity 76:186–212, 2016. 73 https://www.ansys.com/ Acta Polytechnica CTU Proceedings 15:69–73, 2018 1 Introduction 2 Experimental part 3 Numerical modelling 4 Results and discussion 5 Conclusions List of symbols Acknowledgements References