A stepwise-equilibrium and adaptive molecular dynamics (MD) simulation scheme for investigating the fracture toughness of single crystals is proposed in this study. The critical fracture toughness is found by conducting MD simulations along with the gradually increasing external load. At each load step, an equilibrium state is obtained by relaxing the system from the initial state generated. This is done by adjusting the atomic position using an additional displacement of linear elastic solution corresponding to the current load increment. The load increment is adjusted at each step in an adaptive way in order to achieve high computational efficiency and accuracy. A nickel crystal having 14256 atoms is investigated using this technique. The critical stress intensity factor in the (1[UNKNOWN]0) plane is found to be 0.7436 MPa √m, while the fracture stress is 4.7776 GPa. The effects of vacancies on the critical stress intensity factors are also investigated.
Abraham, F.F.2001. “Crack Dynamics in Brittle Fracture: An Atomistic Study,”Nuclear Instruments and Methods in Physics Research, B180: 72–76 .
2.
Bachlechner, M.E. , Omeltchenko, A. , Nakano, A. , Kalia, R.K. , Vashishta, P. , Ebbsjo, I. and Madhukar, A.2000. “Dislocation Emission at the Silicon/Silicon Nitride Interface: A Million Atom Molecular Dynamics Simulation on Parallel Computers,”Phys. Rew. Lett., 84: 322–325 .
3.
Carlsson, A.E. and Thomson, R.1998. “Fracture Toughness of Materials: From Atomistics to Continuum Theory,”Solid State Phys., 51: 233–280 .
4.
Cleri, F. , Yip, S. , Wolf, D. and Phillpot, S.R.1997. “Atomic-scale Mechanism of Crack-tip Plasticity: Dislocation Nucleation and Crack-tip Shielding,”Phys. Rew. Lett., 79: 1309–1312 .
5.
Crawford, E.D. and Lesser, A.J.1999. “Brittle to Ductile: Fracture Toughness Mapping on Controlled Epoxy Networks,”Polym Eng. Sci., 39: 385–392 .
6.
Devincre, B. and Roberts, S.G.1996. “3-D Simulation of Dislocation-crack Interactions in B.C.C. Metals at the Mesoscopic Scale,”Acta Metallurgica et Materialia, 44: 2891–2900 .
7.
Dienes, G.J. and Paskin, A.1983. Latanision, R.M. and Pickens, J.R. (eds), Atomistics of Fracture, p. 671-671, Plenum .
8.
Farkas, D. , Duranduru, M. , Curtin, W.A. and Ribbens, C.2001. “Multiple-dislocation Emission from the Crack tip in the Ductile Fracture of Al,”Phil. Mag., A81: 1241–1255 .
9.
Hauch, J.A. , Holland, D. , Marder, M.P. and Swinney, H.L.1999. “Dynamic Fracture in Single Crystal Silicon,”Phys. Rev., B82: 3823–3826 .
10.
Heino, P. , Hakkinen, H. and Kaski, K.1998. “Molecular dynamics Study of Copper with Defects Under Strain,”Phys. Rev., B58: 641–652 .
11.
Hirosaki, N. , Ogata, S. and Kitagawa, H.1999. “Molecular Dynamics Study of the Stress Distribution Near a Crack tip in β-Silicon Nitride,”Mater. Sci. Res. Int., 5: 253–257 .
12.
Holland, D. and Marder, M.1999. “Cracks and Atoms,”Advanced Materials, 11: 793–806 .
13.
Kitamura, T. , Yashiro, K. and Ohtani, R.1997. “Atomic Simulation on Deformation and Fracture of Anao-single Crystal of Nickel in Tension,”JSME Int. J., 40: 430–435 .
14.
Knap, J. and Sieradzki, K.1999. “Crack Tip Dislocation Nucleation in FCC Solids,”Phys. Rew. Lett., 82: 1700–1703 .
15.
Li, M. and Selinger, R.L.B.2003. “Molecular Dynamics Simulations of Dislocation Instability in a Stress Gradient,”Phys. Rew., B67: 2003-2003 .
16.
Li, Z.J. , Chu, W.Y. , Li, J.X. , Gao, K.W. and Qiao, L.J.2002. “Threedimensional Molecular Dynamics Simulation of Hydrogen-enhanced Dislocation Emission and Crack Propagation,”Prog. Nat. Sci., 12: 717–720 .
17.
Lynden-Bell, R.M.1992. “The Fracture of Perfect Crystals Under Uniaxial Tension at High Temperatures,”J. Phys. Condensed Matter, 4: 2127–2138 .
18.
Lynden-Bell, R.M.1994. “Computer Simulations of Fracture at the Atomic Level,”Science, 263: 1704–1705 .
Meguid, S.A.1989. Engineering Fracture Mechanics, Elsevier Appl. Sci., New York .
21.
Miller, R. , Tadmor, E.B. , Phillips, R. and Ortiz, M.1998a. “Quasicontinuum Simulation of Fracture at the Atom Scale,”Modelling Simul. Mater. Sci. Eng., 6: 607–638 .
22.
Miller, R. , Ortiz, M. , Phillips, R. , Shenoy, V. and Tadmor, E.B.1998b. “Quasicontinuum Models of Fracture and Plasticity,”Eng. Fracture Mech., 61: 427–444 .
23.
Ramasubramaniam, A.Curtin, W.A. and Farkas, D.2002. “Fracture in Nanolamellar Materials: Continuum and Atomistic Models with Application to Titanium Aluminides,”Phil. Mag. A82: 2397–2417 .
24.
Robach, J.S. , Robertson, I.M. , Wirth, B.D. and Arsenlis, A.2003. “In-situ Transmission Electron Microscopy Observations and Molecular Dynamics Simulations of Dislocation-defect Interactions in Ion-irradiated Copper,”Phil. Mag., A83: 955–967 .
25.
Rifkin, J.2000. XMD – Molecular Dynamics Program, University of Connecticut .
26.
Schoeck, G.2003. “The Emission of Dislocations from Crack Tips – A Critical Assessment,”Mater. Sci. Eng., A356: 93–101 .
27.
Tsuji, T. and Noda, N.1991. “Investigations of Stress Singularity and Distraction Mechanism in the Neighborhood of a V-notch Suing Molecular Dynamics Method,”Trans. Japan Soc. Mech. Engineers, A57: 2936–2941 .
28.
Voter, A.F.1994. “The Embedded-atom Method,” In: Westbrook, J.H. and Fleischer, R.L. (eds), Intermetallic Compounds: Principles and Applications, Vol. 1, pp. 77–90, John Wiley & Sons Ltd.
29.
Wunderlich, W. and Awaji, H.2001. “Molecular Dynamics simulation of the Fracture Toughness of Sapphire,”Mater. Design, 22: 53–59 .
30.
Yanagida, N. and Watanabe, O.1996. “Molecular Dynamics Simulation of Effects of Lattice Orientation on Crack Propagation in α-iron when the Primary Slip Direction is in the Plane of Tensile Stress,”JSME Int. J., A39: 321–329 .
31.
Zhou, S.J. , Preston, D.L. , Lomdahl, P.S. and Beazley, D.M.1998a. “Large-scale Molecular Dynamics Simulations of Dislocation Intersection in Copper,”Science, 279: 1525–1527 .
