Numerical study of metal matrix nanocomposites using discrete dislocation approach

NUMERICAL STUDY OF METAL MATRIX NANOCOMPOSITES USING DISCRETE DISLOCATION APPROACH ELLIOT LAW NATIONAL UNIVERSITY OF SINGAPORE 2011 NUMERICAL STUDY OF METAL MATRIX NANOCOMPOSITES USING DISCRETE DISLOCATION APPROACH ELLIOT LAW (B. Eng. (Civil) (Hons), NUS) A THESIS SUBMITTED FOR THE DEGREE OF DOCTOR OF PHILOSOPHY DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING NATIONAL UNIVERSITY OF SINGAPORE 2011 Acknowledgements First and foremost, praise be to God for His goodness and faithfulness in bringing me through this entire study. I thank God for blessing me with wisdom and sustaining me with good health as I worked on this project. I would also like express my deepest thanks to my parents and family members for their unconditional love and support throughout all these years. I would like to express my utmost gratitude to my supervisors, Dr. Pang Sze Dai and Prof. Quek Ser Tong, for their advice, guidance and counsel since my undergraduate years. I have learnt so much from them and it has been such a great pleasure and joy to be able to work with them. I am also very grateful for their care and support; they have inspired me to likewise toward others and to give my best in all my undertakings. I would also like to thank the National University of Singapore (NUS) for supporting me with the Research Scholarship for the entire duration of my study. I would also like to acknowledge my fellow students whom I have had the opportunity to work with throughout the duration of my studies. A special note of appreciation to Mr. Tran Diep Phuoc Thao, Ms. Liu Lihui, Ms. Matilda Loh and Mr. Too Jun Lin for their friendship, encouragement and support; the joy of working with them has made my work more meaningful as well as given me the impetus and motivation to complete this pursuit. i Last but not least, I would like to express my thanks to the staff at the Structural Engineering Laboratory, especially Mdm. Annie Tan (who has since transferred to the Engineering Design and Innovation Centre) and Mr. Ang Beng Oon, for their support in this project. I would also like to acknowledge Dr. Sharon Nai Mui Ling from the Singapore Institute of Manufacturing Technology (SIMTech) for her tremendous assistance in the experimental work conducted for this study. Special thanks also to Assoc. Prof. Manoj Gupta and Dr. Khin Sandar Tun from the Department of Mechanical Engineering at NUS as well as Mr. Lam Kim Song from the Fabrication Support Centre of the same department for their help in the fabrication and machining of the nanocomposite specimens as well as sharing of experimental results on tensile properties and microstructural characteristics. Give thanks to the Lord, for He is good; His love endures forever. (Ps. 107:1) ii Table of Contents i Acknowledgements iii Table of Contents viii Summary List of Tables xi List of Figures xii List of Symbols xxii 1.0 Introduction 1.1 Metal matrix composites – overview 1.2 Numerical studies on metal matrix composites 1.2.1 Metal matrix composite systems 1.2.2 Influence of various microstructural features on mechanical properties of metal matrix composites 1.2.3 Microstructural modelling using representative volume element approach 1.2.4 Boundary conditions for representative volume element 1.3 Metal matrix nanocomposites – overview 1.4 Size effects on mechanical properties of metal matrix composites 13 1.5 Methods of computer simulations for dislocations 15 1.5.1 Atomistic methods 16 1.5.2 Continuum methods 17 1.5.3 Extended finite element method 19 iii 1.5.4 2.0 iv 20 1.6 Numerical simulations of metal matrix nanocomposites 21 1.7 Objective 25 1.8 Scope and limitations 25 1.9 Organization of the study 26 Theoretical framework for composite material model using discrete dislocation method 27 2.1 Dislocations – basic concepts and characteristics 27 2.2 Discrete dislocation formulation 32 2.2.1 Instantaneous state of dislocated body 32 2.2.2 Forces between dislocations 35 2.2.3 Constitutive relations for motion of dislocations 36 2.2.4 Constitutive relations for creation and annihilation of dislocations 39 Analytical solutions for dislocation fields in two-dimensional space 43 2.3.1 Special case of two-dimensional space with horizontal slip planes 45 2.4 Implementation of discrete dislocation formulation for twodimensional unit cell analyses 50 2.5 Computational scheme 54 2.3 3.0 Multi-scale modelling and coupled atomisticcontinuum methods Numerical implementation issues 59 3.1 Computational time-step and efficiency 59 3.1.1 59 Evaluation of boundary nodal forces due to dislocation stress field 3.1.2 3.2 3.3 4.0 Effect of cut-off distance on calculation of dislocation fields 63 3.1.3 Evaluation of dislocation glide force 65 3.1.4 Time-step 67 3.1.5 Tracking of dislocation events and processes 69 Calibration of material parameters for dislocation processes 72 3.2.1 Density and strength of impurities or obstacles 74 3.2.2 Density of dislocation sources 80 3.2.3 Nucleation strength of dislocation sources 82 3.2.4 Calibration procedure 84 Size of representative volume element 90 3.3.1 Density of dislocation sources 94 3.3.2 Strength and density of impurities 98 3.3.3 Change in mean overall response with number of realizations 100 3.3.4 Overall response of composite material 104 Effects of microstructural features and constituent material properties on mechanical response of metal matrix nanocomposites 111 4.1 Inclusion volume fraction 111 4.2 Inclusion size 115 4.3 Inclusion aspect ratio and orientation 118 4.4 Arrangement of inclusions 124 4.5 Material properties of constituent phases 132 4.6 Nature of dislocation pile-ups in metallic nanocomposites 137 v 5.0 Simulation of damage in metal matrix nanocomposites 141 5.1 Damage of inclusions 144 5.1.1 Effects of fracture strength, volume fraction and size of inclusions 145 5.1.2 156 5.2 6.0 vi Matrix damage 160 5.2.1 Effect of void formation in a pure metallic matrix 162 5.2.2 Effect of void formation in a metallic nanocomposite 166 5.2.3 Effect of inclusion volume fraction 170 5.2.4 174 Effect of inclusion arrangement Experimental verification 181 6.1 Experimental method 181 6.1.1 189 Fixture assembly and installation 6.1.2 Testing procedure 191 6.2 Materials processing and specimen fabrication 193 6.3 Experimental results 196 6.3.1 196 6.4 7.0 Effect of inclusion arrangement Hardness and tensile properties 6.3.2 Stress-strain response under shear 196 Comparison with numerical results 201 Conclusions and future work 205 7.1 Conclusions 205 7.2 Recommendations for future work 208 7.2.1 Effect of interfacial zone 208 7.2.2 Consideration of crystallographic details 213 References 217 Appendix A: Modifications to standard Iosipescu shear test fixture for testing of smaller specimens 231 Appendix B: List of publications 235 vii Summary A metal matrix composite (MMC) is a composite material with reinforcement phase which is dispersed within a continuous metallic host to improve the thermomechanical properties of the host metal. Recent experiments show that reducing the reinforcement size to the nanoscale dramatically increases the mechanical strength of MMCs. While extensive numerical studies on the mechanical properties of conventional MMCs have been conducted, only a handful of such studies exist for metal matrix nanocomposites (MMNCs). Numerical simulations are useful for performing virtual experiments on MMNCs to explore effects which are currently difficult to investigate experimentally and to analyse the underlying processes that govern the mechanical response of these materials. Hence, the objective of this study is to investigate the mechanical properties of MMNCs using numerical simulations in order to determine the best combinations of constituent material properties, compositions and microstructure for optimum mechanical performance of these materials. Two-dimensional discrete dislocation analysis implemented using the multi-inclusion representative volume element (RVE) approach is adopted for the simulations. A calibration procedure is developed to determine the suitable values for various parameters which describe dislocation processes in a pure metallic matrix. Suitable RVE sizes required for the modelling of MMNCs are also established; it is found that a statistically representative RVE should have approximately 70 to 80 inclusions. viii Fleck, N. A., Muller, G. M., Ashby, M. F. and Hutchinson, J. W. 1994. Strain gradient plasticity: theory and experiment. Acta Metallurgica et Materialia 42(2): 475-487. Fleck, N. A. and Hutchinson, J. W. 2001. A reformulation of strain gradient plasticity. Journal of the Mechanics and Physics of Solids 49: 2245-2271. Freund, L.B. 1994. The mechanics of dislocations in strained-layer semiconductor materials. Advances in Applied Mechanics 30: 1-66. G’Sell, C., Boni, S. and Shrivastava, S. 1983. Application of the plane simple shear test for determination of the plastic behaviour of solid polymers at large strains. Journal of Materials Science 18: 903-918. Ganesh, V. V. and Chawla, N. 2005. Effect of particle orientation anisotropy on the tensile behavior of metal matrix composites: experiments and microstructure-based simulation. Materials Science and Engineering A 391: 342-353. Gao, H., Huang, Y., Nix, W. D. and Hutchinson, J. W. 1999. Mechanism-based strain gradient plasticity - I. Theory. Journal of the Mechanics and Physics of Solids 47: 1239-1263. Gao, H. and Huang, Y. 2003. Geometrically necessary dislocation and sizedependent plasticity. Scripta Materialia 48: 113-118. Ghosh, S. and Moorthy, S. 1998. Particle fracture simulation in non-uniform microstructures of metal-matrix composites. Acta Materialia 46(3): 965-982. Gitman, I.M., Gitman, M.B. and Askes, H. 2006. Quantification of stochastically stable representative volumes for random heterogeneous materials. Archive of Applied Mechanics 75: 79-92. Gitman, I.M., Askes, H. and Sluys, L.J. 2007. Representative volume: existence and size determination. Engineering Fracture Mechanics 74(16): 2518-2534. Gnosh, S., Bai, J. and Raghavan, P. 2007. Concurrent multi-level model for damage evolution in microstructurally debonding composites. Mechanics of Materials 39: 241-266. Gracie, R., Ventura, G. and Belytschko, T. 2007. A new finite element method for dislocations based on interior discontinuities. International Journal for Numerical Methods in Engineering 69: 423-441. Gracie, R., Oswald, J. and Belytschko, T. 2008. On a new extended finite element method for dislocations: core enrichment and nonlinear formulation. Journal of the Mechanics and Physics of Solids 56: 200-214. Greengard, L. and Rokhlin, V. 1987. A fast algorithm for particle simulations. Journal of Computational Physics 73: 325-348. 221 Groh, S., Devincre, B., Kubin, L.P., Roos, A., Feyel, F. and Chaboche, J.L. 2005. Size effects in metal matrix composites. Materials Science and Engineering A 400401: 279-282. Gulluoglu, A.N., Srolovitz, D.J., LeSar, R. and Lomdahl, P.S. 1989. Dislocation distributions in two dimensions. Scripta Metallurgica 23(8): 1347-1352. Gupta, M. and Wong, W. L. E. 2005. Enhancing overall mechanical performance of metallic materials using two-directional microwave assisted rapid sintering. Scripta Materialia 52: 479-483. Han, W., Eckschlager, A. and Bohm, H. J. 2001. The effects of three-dimensional multi-particle arrangements on the mechanical behavior and damage initiation of particle-reinforced MMCs. Composites Science and Technology 61: 1581-1590. Harik, V. M. and Salas, M. D., eds. 2003. Trends in nanoscale mechanics: analysis of nanostructured materials and multi-scale modeling. ICASE/LaRC Interdisciplinary Series in Science and Engineering. Dordrecht, The Netherlands: Kluwer Academic Publishers. Hasnaoui, A., Derlet, P. M. and Van Swygenhoven, H. 2004. Interaction between dislocations and grain boundaries under an indenter - a molecular dynamics simulation. Acta Materialia 52(8): 2251-2258. Hassan, S. F. and Gupta, M. 2005. Development of high performance magnesium nano-composites using nano-Al2O3 as reinforcement. Materials Science and Engineering A 392: 163-168. Hassan, S. F. and Gupta, M. 2006a. Effect of length scale of Al2O3 particulates on microstructural and tensile properties of elemental Mg. Materials Science and Engineering A 425: 22-27. Hassan, S. F. and Gupta, M. 2006b. Effect of different types of nano-size oxide particulates on microstructural and mechanical properties of elemental Mg. Journal of Materials Science 41: 2229-2236. Hbaieb, K., Wang, Q. X., Chia, Y. H. J. and Cotterell, B. 2007. Modelling stiffness of polymer/clay nanocomposites. Polymer 48: 901-909. Heness, G. L., Ben-Nissan, B., Gan, L. H. and Mai, Y.W. 1999. Development of a finite element micromodel for metal matrix composites. Computational Materials Science 13: 259-269. Hsu, C. J., Chang, C. Y., Kao, P. W., Ho, N. J. and Chang, C. P. 2006. Al-Al3Ti nanocomposites produced in situ by friction stir processing. Acta Materialia 54(19): 5241-5249. Huang, Y., Gao, H., Nix, W. D. and Hutchinson, J. W. 2000. Mechanism-based strain gradient plasticity - II. Analysis. Journal of the Mechanics and Physics of Solids 48: 99-128. 222 Huang, Y., Qu, S., Hwang, K. C., Li, M. and Gao, H. 2004. A conventional theory of mechanism-based strain gradient plasticity. International Journal of Plasticity 20: 753-782. Hull, D. and Bacon, D.J. 1984. Introduction to dislocations. 3rd ed. Oxford: Pergamom Press. Humphreys, J. and Bate, P. 2003. Gradient plasticity and deformation structures around inclusions. Scripta Materialia 48: 173-178. Iosipescu, N. 1967. New accurate procedure for single shear testing of metals. Journal of Materials 2(3): 537-566. Jiang, M., Ostoja-Starzewski, M. and Jasiuk, I. 2001. Scale-dependent bounds on effective elastoplastic response of random composites. Journal of the Mechanics and Physics of Solids 49: 655-673. Jonsson, A. 2003. Discrete dislocation dynamics by an O(N) algorithm. Computational Materials Science 27: 271-288. Kalpakjian, S. and Schmid, S.R. 2006. Manufacturing engineering and technology. 5th ed. Upper Saddle River, NJ: Prentice Hall. Kang, Y. C. and Chan, S. L. I. 2004. Tensile properties of nanometric Al2O3 particulate-reinforced aluminum matrix composites. Materials Chemistry and Physics 85: 438-443. Kanit, T., Forest, S., Galliet, I., Mounoury, V. and Jeulin, D. 2003. Determination of the size of the representative volume element for random composites: statistical and numerical approach. International Journal of Solids and Structures 40: 3647-3679. Kenesei, P., Borbély, A. and Biermann, H. 2004. Microstructure based threedimensional finite element modelling of particulate reinforced metal-matrix composites. Materials Science and Engineering A 387-389: 852-856. Khisaeva, Z.F. and Ostoja-Starzewski, M. 2006. On the size of RVE in finite elasticity of random composites. Journal of Elasticity 85: 153-173. Knap, J. and Ortiz, M. 2001. An analysis of the quasicontinuum method. Journal of the Mechanics and Physics of Solids 49: 1899-1923. Kohler, C., Kizler, P. and Schmauder, S. 2005. Atomistic simulation of the pinning of edge dislocations in Ni by Ni3Al precipitates. Materials Science and Engineering A 400-401: 481-484. Kubin, L.P., Canova, G., Condat, M., Devincre, B., Pontikis, V. and Brechet, Y. 1992. Dislocation microstructures and plastic flow: a 3D simulation. Solid State Phenomena 23-24: 455-472. 223 Kuhn, H.A. 2000. Shear, torsion and multiaxial testing. In ASM handbook 8, mechanical testing and evaluation, 185-194. ASM International. Kumar, R., Nicola, L. and Van der Giessen, E. 2009. Density of grain boundaries and plasticity size effects: a discrete dislocation dynamics study. Materials Science and Engineering A 527: 7-15. Lee, H. K. and Pyo, S. H. 2008. Multi-level modeling of effective elastic behavior and progressive weakened interface in particulate composites. Composites Science and Technology 68: 387-397. Li, M., Chu, W. Y., Qian, C. F., Gao, K. W. and Qiao, L. J. 2003. Molecular dynamics simulation of dislocation intersections in aluminum. Materials Science and Engineering A 363(1-2): 234-241. Liu, H. B., Canizal, G., Jimenez, S., Espinosa-Medina, M. A. and Ascensio, J. A. 2003. Molecular dynamics simulation on edge dislocation in the bulk and nanoparticles of iron. Computational Materials Science 27(3): 333-341. Liu, Y., Van der Giessen, E. and Needleman, A. 2007. An analysis of dislocation nucleation near a free surface. International Journal of Solids and Structures 44(6): 1719-1732. Liu, X. Y. and Biner, S. B. 2008. Molecular dynamics simulations of the interactions between screw dislocations and self-interstitial clusters in body-centered cubic Fe. Scripta Materialia 59(1): 51-54. Llorca, J., Needleman, A., Suresh, S. 1991. An analysis of the effects of matrix void growth on deformation and ductility in metal-ceramic compsoites. Acta Metallurgical et Materialia 39(10): 2317-2335. Llorca, J. and González, C. 1998. Microstructural factors controlling the strength and ductility of particle-reinforced metal-matrix composites. Journal of the Mechanics and Physics of Solids 46(1): 1-18. Lubarda, V.A., Blume, J.A and Needleman, A. 1993. An analysis of equilibrium dislocation distributions. Acta Metallurgica et Materialia 41(2): 625-642. Luo, J. J. and Daniel, I. M. 2003. Characterization and modeling of mechanical behavior of polymer/clay nanocomposites. Composites Science and Technology 63: 1607-1616. Luo, D., Wang, W. X. and Takao, Y. 2007. Effects of the distribution and geometry of carbon nanotubes on the macroscopic stiffness and microscopic stresses of nanocomposites. Composites Science and Technology 67(14): 2947-2958. Ma, Q. and Clarke, D. R. 1995. Size dependent hardness of silver single crystals. Journal of Materials Research 10: 853-863. 224 Ma, Z.Y., Li, Y.L., Liang, Y., Zheng, F., Bi, J. and Tjong, S.C. 1996. Nanometric Si3N4 particulate-reinforced aluminum composite. Materials Science and Engineering A 219: 229-231. Ma, Z.Y., Tjong, S.C., Li, Y.L. and Liang, Y. 1997. High temperature creep behavior of nanometric Si3N4 particulate reinforced aluminum composite. Materials Science and Engineering A 225: 125-134. Martin, L.P., Dadon, D. and Rosen, M. 1996. Evaluation of ultrasonically determined elasticity-porosity relations in zinc oxide. Journal of the American Ceramic Society 79(5): 1281-1289. Mazahery, A., Abdizadeh, H. and Baharvandi, H.R. 2009. Development of highperformance A356/nano-Al2O3 composites. Materials Science and Engineering A 518: 61-64. Miller, R., Tadmor, E. B., Phillips, R. and Ortiz, M. 1998. Quasicontinuum simulation of fracture at the atomic scale. Modelling and Simulation in Materials Science and Engineering 6: 607-638. Mishnaevsky, L. Jr., Dong, M., Honle, S. and Schmauder, S. 1999. Computational mesomechanics of particle-reinforced composites. Computational Materials Science 16: 133-143. Mishnaevsky, L.L. Jr. 2004. Three-dimensional numerical testing of microstructures of particle reinforced composites. Acta Materialia 52: 4177-4188. Mishnaevsky, L. Jr., Derrien, K. and Baptiste, D. 2004. Effect of microstructure of particle reinforced composites on the damage evolution: probabilistic and numerical analysis. Composites Science and Technology 64: 1805-1818. Mishnaevsky, L. L. Jr. 2006. Functionally gradient metal matrix composites: Numerical analysis of the microstructure-strength relationships. Composites Science and Technology 66: 1873-1887. Moes, N., Dolbow, J. and Belytschko, T. 1999. A finite element method for crack growth without remeshing. International Journal for Numerical Methods in Engineering 46: 131-150. Mondal, D. P., Ramakrishnan, N. and Das, S. 2006. FEM modeling of the interface and its effect on the elasto-plastic behavior of metal matrix composites. Materials Science and Engineering A 433: 286-290. Mughrabi, H. 2004. On the current understanding of strain gradient plasticity. Materials Science and Engineering A 387-389: 209-213. Muskhelishvili, N.I. 1953. Some basic problems of the mathematical theory of elasticity: fundamental equations, plane theory of elasticity, torsion and bending. Groningen, Holland: P. Noordhoff Ltd. 225 Nabarro, F.R.N. 1952. The mathematical theory of stationary dislocations. Advances in Physics 1(3): 269-394. Needleman, A. and Van der Giessen, E. 2001. Discrete dislocation and continuum descriptions of plastic flow. Materials Science and Engineering A 309-310: 1-13. Needleman, A., Van der Giessen, E. and Deshpande, V.S. 2006. Statistical aspects of discrete dislocation plasticity. Scripta Materialia 54: 729-733. Nicholas, T. 1971. Strain-rate and strain-rate-history effects in several metals in torsion. Experimental Mechanics 11(8): 370-374. Nicola, L., Xiang, Y., Vlassak, J.J., Van der Giessen, E. and Needleman, A. 2006. Plastic deformation of freestanding thin films: experiments and modeling. Journal of the Mechanics and Physics of Solids 54: 2089-2110. Odegard, G. M., Clancy, T. C. and Gates, T. S. 2005. Modeling of the mechanical properties of nanoparticle/polymer composites. Polymer 46: 553-562. Odegard, G. M., Gates, T. S., Wise, K. E., Park, C. and Siochi, E. J. 2003. Constitutive modeling of nanotube-reinforced polymer composites. Composites Science and Technology 63: 1671-1687. Pelissou, C., Baccou, J., Monerie, Y. and Perales, F. 2009. Determination of the size of the representative volume element for random quasi-brittle composites. International Journal of Solids and Structures 46: 2842-2855. Poole, W. J., Ashby, M. F. and Fleck, N. A. 1996. Micro-hardness of annealed and work-hardened copper polycrystals. Scripta Materialia 34: 559-564. Qiao, R. and Brinson, L.C. 2009. Simulation of interphase percolation and gradients in polymer nanocomposites. Composites Science and Technology 69: 491-499. Qu, S., Siegmund, T., Huang, Y., Wu, P. D., Zhang, F. and Hwang, K. C. 2005. A study of particle size effect and interface fracture in aluminum alloy composite via an extended conventional theory of mechanism-based strain-gradient plasticity. Composites Science and Technology 65: 1244-1253. Ramasubramaniam, A., Ariza, M.P. and Ortiz, M. 2007. A discrete mechanics approach to dislocation dynamics in BCC crystals. Journal of the Mechanics and Physics of Solids 55: 615-647. Rintoul, M.D. and Torquato, S. 1997. Reconstruction of the structure of dispersions. Journal of Colloid and Interface Science 186(2): 467-476. Rosenberger, M. R., Forlerer, E. and Schvezov, C. E. 2007. Modeling the microindentation of metal matrix composites. Materials Science and Engineering A 463: 275-283. 226 Segurado, J. and Llorca, J. 2002. A numerical approximation to the elastic properties of sphere-reinforced composites. Journal of the Mechanics and Physics of Solids 50(10): 2107-2121. Segurado, J., González, C. and Llorca, J. 2003. A numerical investigation of the effect of particle clustering on the mechanical properties of composites. Acta Materialia 51: 2355-2369. Segurado, J. and Llorca, J. 2005. A computational micromechanics study of the effect of interface decohesion on the mechanical behavior of composites. Acta Materialia 53(18): 4931-4942. Selmi, A., Friebel, C., Doghri, I. and Hassis, H. 2007. Prediction of the elastic properties of single walled carbon nanotube reinforced polymers: a comparative study of several micromechanical models. Composites Science and Technology 67: 20712084. Shati, F. K., Esat, I. I. and Bahai, H. 2001. FEA modelling of visco-plastic behavior of metal matrix composites. Finite Elements in Analysis and Design 37: 263-272. Shilkrot, L. E., Miller, R. E. and Curtin, W. A. 2004. Multiscale plasticity modeling: couple atomistics and discrete dislocation mechanics. Journal of the Mechanics and Physics of Solids 52: 755-787. Shin, C.S., Fivel, M.C., Verdier, M and Kwon, S.C. 2006. Numerical methods to improve the computing efficiency of discrete dislocation dynamics simulations. Journal of Computational Physics 215(2): 417-429. Shu, J. Y. and Fleck, N. A. 1999. Strain gradient crystal plasticity: size-dependent deformation of bicrystals. Journal of the Mechanics and Physics of Solids 47(2): 297324. Shu, J.Y., Fleck, N.A., Van der Giessen, E. and Needleman, A. 2001. Boundary layers in constrained plastic flow: comparison of nonlocal and discrete dislocation plasticity. Journal of the Mechanics and Physics of Solids 49: 1361-1395. Song, M. and Xiao, D. 2008. Modeling the fracture toughness and tensile ductility of SiCp/Al metal matrix composites. Materials Science and Engineering A 474: 371375. Soppa, E., Schmauder, S., Fischer, G., Brollo, J. and Weber, U. 2003. Deformation and damage in Al/Al2O3. Computational Materials Science 28: 574-586. Stewart, D. and Cheong, K. S. 2008. Molecular dynamics simulations of dislocations and nanocrystals. Current Applied Physics 8(3-4): 494-497. Stolken, J. S. and Evans, A. G. 1998. A microbend test method for measuring the plasticity length scale. Acta Materialia 46: 5109-5115. 227 Stroeven, M., Askes, H. and Sluys, L.J. 2004. Numerical determination of representative volumes for granular materials. Computer Methods in Applied Mechanics and Engineering 193: 3221-3238. Su, X. F., Chen, H. R., Kennedy, D. and Williams, F. W. 1999. Effects of interphase strength on the damage modes and mechanical behaviour of metal-matrix composites. Composites: Part A 30: 257-266. Sun, L.Z., Liu, H.T. and Ju, J.W. 2003. Effect of particle cracking on elastoplastic behaviour of metal matrix composites. International Journal for Numerical Methods in Engineering 56: 2183-2198. Tadmor, E. B., Phillips, R. and Ortiz, M. 1996. Mixed atomistic and continuum models of deformation in solids. Langmuir 12(6): 4529-4534. Taya, M. and Arsenault, R.J. 1989. Metal matrix composites – thermomechanical behavior. New York, NY: Pergamom Press. Tian, L. and Rajapakse, R.K.N.D. 2007. Finite element modelling of nanoscale inhomogeneities in an elastic matrix. Computational Materials Science 41: 44-53. Tjong, S. C. 2007. Novel nanoparticle-reinforced metal matrix composites with enhanced mechanical properties. Advanced Engineering Materials 9(8): 639-652. Tun, K.S., Jayaramanavar, P., Nguyen, Q.B., Chan, J., Kwok, R. and Gupta, M. 2011. Investigation into tensile and compressive responses of Mg-ZnO composites. Materials Science and Technology (accepted for publication). Van der Giessen, E. and Needleman, A. 1995. Discrete dislocation plasticity: a simple planar model. Modelling and Simulation in Materials Science and Engineering 3: 689-735. Van der Giessen, E. and Needleman, A. 2003. GNDs in nonlocal plasticity theories: lessons from discrete dislocation simulations. Scripta Materialia 48: 127-132. Ventura, G., Moran, B. and Belytschko, T. 2005. Dislocations by partition of unity. International Journal for Numerical Methods in Engineering 62: 1463-1487. Voyiadjis, G. Z., Aifantis, E. C. and Weber, G. 2003. Constitutive modeling of plasticity in nanostructured materials. In Trends in nanoscale mechanics: analysis of nanostructured materials and multi-scale modeling, eds. Harik, V.M. and Salas, M. D., 123-146. Dordrecht, The Netherlands: Kluwer Academic Publishers. Walrath, D.E. and Adams, D.F. 1983. The Iosipescu shear test as applied to composite materials. Experimental Mechanics 23(1): 105-110. Walrath, D.E. and Adams, D.F. 1985. Iosipescu shear properties of graphite fabric/epoxy composite laminates. Report no. UMWE-DR-501-103-1, Department of Mechanical Engineering, University of Wyoming, NASA Grant no. NAG-1-272. Also available as NASA-CR-176316. 228 Wang, H.Y. and LeSar, R. 1995. O(N) algorithm for dislocation dynamics. Philosophical Magazine A 71(1): 149-164. Wang, Y. U., Jin, Y. M., Cuitino, A. M. and Khachaturyan, A. G. 2001. Nanoscale phase field microelasticity theory of dislocations: model and 3D simulations. Acta Materialia 49(10): 1847-1857. Ward, D. K., Curtin, W. A. and Qi, Y. 2006. Mechanical behavior of aluminumsilicon nanocomposites: A molecular dynamics study. Acta Materialia 54: 44414451. Weibull, W. 1951. A statistical distribution function of wide applicability. Journal of Applied Mechanics 18: 293-297. Wong, W. L. E., Karthik, S. and Gupta, M. 2005. Development of high performance Mg-Al2O3 composites containing Al2O3 in submicron length scale using microwave assisted rapid sintering. Materials Science and Technology 21(9): 1063-1070. Wong, W. L. E. and Gupta, M. 2006. Simultaneously improving strength and ductility of magnesium using nano-size SiC particulates and microwaves. Advanced Engineering Materials 8(8): 735-740. Wong, W. L. E. and Gupta, M. 2007a. Improving overall mechanical performance of magnesium using nano-alumina reinforcement and energy efficient microwave assisted processing route. Advanced Engineering Materials 9(10): 902-909. Wong, W. L. E. and Gupta, M. 2007b. Development of Mg/Cu nanocomposites using microwave assisted rapid sintering. Composites Science and Technology 67: 1541-1552. Xia, Z. H. and Curtin, W. A. 2001. Multiscale modeling of damage and failure in aluminum-matrix composites. Composites Science and Technology 61: 2247-2257. Xia, Z., Curtin, W. A. and Peters, P. W. M. 2001. Multiscale modeling of failure in metal matrix composites. Acta Materialia 49: 273-287. Xue, Z., Huang, Y. and Li, M. 2002. Particle size effect in metallic materials: a study by the theory of mechanism-based strain gradient plasticity. Acta Materialia 50: 149160. Yamakov, V., Wolf, D., Phillpot, S. R. and Gleiter, H. 2002. Deformation twinning in nanocrystalline Al by molecular-dynamics simulation. Acta Materialia 50(20): 5005-5020. Yamakov, V., Wolf, D., Phillpot, S. R. and Gleiter, H. 2003. Dislocation-dislocation and dislocation-twin reactions in nanocrystalline Al by molecular dynamics simulation. Acta Materialia 51(14): 4135-4147. 229 Yan, Y. W., Geng, L. and Li, A. B. 2007. Experimental and numerical studies of the effect of particle size on the deformation behavior of the metal matrix composites. Materials Science and Engineering A 448: 315-325. Yang, J., Cady, C., Hu, M.S., Zok, F., Mehrabian, R. and Evans, A.G. 1990. Effects of damage on the flow strength and ductility of a ductile Al alloy reinforced with SiC particulates. Acta Metallurgical et Materialia 38(12): 2613-2619. Yang, Y, Lan, J. and Li, X. 2004. Study on bulk aluminum matrix nano-composite fabricated by ultrasonic dispersion of nano-sized SiC particles in molten aluminum alloy. Materials Science and Engineering A 380: 378-383. Yao, Z., Wang, J.S., Liu, G.R., Cheng, M. 2004. Improved neighbor list algorithm in molecular simulations using cell decomposition and data sorting method. Computer Physics Communications 161: 27-35. Zhang, Z. and Chen, D. L. 2006. Consideration of Orowan strengthening effect in particulate-reinforced metal matrix nanocomposites: a model for predicting their yield strength. Scripta Materialia 54: 1321-1326. Zhang, J. T., Liu, L. S., Zhai, P. C., Fu, Z. Y. and Zhang, Q. J. 2007a. The prediction of the dynamic responses of ceramic particle reinforced MMCs by using multiparticle computational micro-mechanical model. Composites Science and Technology 67: 2775-2785. Zhang, W. X., Li, L. X. and Wang, T. J. 2007b. Interphase effect on the strengthening behavior of particle-reinforced metal matrix composites. Computational Materials Science 41: 145-155. Zhang, Z. and Chen, D. L. 2008. Contribution of Orowan strengthening effect in particulate-reinforced metal matrix nanocomposites. Materials Science and Engineering A 483-484: 148-152. 230 Appendix A: Modifications to standard Iosipescu shear test fixture for testing of smaller specimens The standard specimen dimensions for the Iosipescu shear test are shown in Figure 6.2 in which the total width is 20 mm while the nominal width of the gauge region is 12 mm. However, the total width of the Mg-ZnO nanocomposite specimen used in this study is restricted to only mm, with the nominal width between the notches reduced proportionally to mm. Hence, additional modifications are made to the test fixture in order to accommodate a smaller specimen without resorting to using a scaled-down fixture which is much more expensive than the standard-size fixture. As shown in Figure 6.7, additional metal pieces are required to reduce the distance between the contact points adjacent to the notches (i.e. reduce the length b shown in Figure 6.4) to prevent premature failure of the smaller specimen due to high bending stresses, as well as to fill in the gap in the grip regions due to the reduced specimen width. Calculations to determine the dimensions of these additional metal pieces are shown here. As shown in Figure 6.4, the maximum bending moment in the specimen is given by M max  Pb (A.1) Since the cross-section of the specimen is rectangular, the elastic section modulus S of the cross-section in the grip region about the (horizontal) bending axis is given by S hd1 (A.2) 231 where h is the thickness of the specimen while d1 is the total width as shown in Figure 6.2. Assuming that the shear stress τ is constant along the cross-section of the specimen in the notched region, the relationship between τ and the applied load P is given by P   hw (A.3) where w is the width of the specimen across the notched region. Therefore, if the grip region remains elastic throughout the duration of the test, the relationship between the maximum bending stress σmax in the grip region and the shear stress τ in the notched region can be determined as follows: M max 3bw  2 S d1  max  (A.4) To prevent premature failure of the grip region, the maximum bending stress σmax must be less than the yield strength σy of the specimen. Assuming that shear strength y y and substituting into Equation (A.4), the following relationship between b and the specimen dimensions can be obtained: bw d1  3  bw d1  0.577 (A.5) For the standard test fixture used in this study, b = 0.5 in ≈ 12.7 mm. Hence, for the standard-size specimen, bw d1 232  12.7  12  0.381  0.577 20 (A.6) Consequently, the standard-size specimen will not fail prematurely due to bending in the grip region with a factor of safety of 1.5. For the smaller specimen, d1 and w are limited to mm and mm respectively. Therefore, the corresponding limit for b is given by b 0.577d1 0.577    5.2 mm w (A.7) Hence, a value of b = 0.15 in ≈ 3.8 mm is suggested for the smaller specimen to prevent bending failure in the grip region. The corresponding factor of safety is approximately 1.35. Based on the suggested value of b and the dimensions of the smaller specimen, the dimensions of the additional metal pieces are determined as shown in Figure A.1. These additional metal pieces are inserted into the test fixture as shown in Figure 6.7 when testing of the smaller specimen is conducted. 0.500 1.625 0.125 0.550 0.550 1.800 0.250 0.075 (Units: inches) Figure A.1 Dimensions for additional metal pieces for testing of smaller specimens. 233 This page is intentionally left blank. 234 Appendix B: List of publications Law, E., Pang, S.D. and Quek, S.T. 2009. Numerical study of metal matrix nanocomposites using discrete dislocation approach. 22nd KKCNN Symposium on Civil Engineering, Chiang Mai, Thailand, October 31 – November 2. Quek, S.T., Law, E. and Pang, S.D. 2010. Mechanical response of metal matrix nanocomposites using discrete dislocation approach. 18th International Conference on Composites or Nano Engineering (ICCE-18), Anchorage, Alaska, USA, July 4-10. Law, E., Pang, S.D. and Quek, S.T. 2010. Effect of thermal residual stress on the mechanical response of metal matrix composites with nanosized reinforcement particles. 23rd KKCNN Symposium on Civil Engineering, Taipei, Taiwan, November 13-15. Law, E., Pang, S.D. and Quek, S.T. 2010. Effect of RVE size in the discrete dislocation simulation of the mechanical response of metal matrix composites with nanosized reinforcement particles. 7th Asian-Australasian Conference on Composite Materials (ACCM-7), Taipei, Taiwan, November 15-18. Law, E., Pang, S.D. and Quek, S.T. 2011. Discrete dislocation analysis of the mechanical response of silicon carbide reinforced aluminum nanocomposites. Composites Part B: Engineering 42(1): 92-98. Law, E., Pang, S.D. and Quek, S.T. 2011. Numerical analysis of the effect of particle arrangement on mechanical behavior and particle damage in metal matrix nanocomposites. 18th International Conference on Composite Materials (ICCM18), Jeju Island, Korea, August 21-26. Law, E., Pang, S.D. and Quek, S.T. 2012. Effects of particle arrangement and particle damage on the mechanical response of metal matrix nanocomposites: a numerical analysis. Acta Materialia 60(1): 8-21. 235 This page is intentionally left blank. 236 [...]... Distribution of dislocations within metallic matrix with τobs = 0.60 GPa and ρnuc of (a) 40 μm-2, and (a) 160 μm-2 82 xiii Figure 3.13 Mean overall response of metallic matrix for different values of τ*nuc with τobs of (a) 0.15 GPa, and (b) 1.20 GPa 83 Figure 3.14 Overall response of metallic matrix with different deformation stages indicated 87 Figure 3.15 Distribution of dislocations within metallic matrix. .. and τobs of (a) 0.15 GPa, and (b) 0.60 GPa 76 Figure 3.9 Mean overall response of metallic matrix for different values of τobs with ρobs of (a) 80 μm-2, and (b) 160 μm-2 78 Figure 3.10 Distribution of dislocations within metallic matrix with τobs = 0.60 GPa and ρobs of (a) 80 μm-2, and (a) 320 μm-2 80 Figure 3.11 Mean overall response of metallic matrix for different values of ρnuc with τobs of (a) 0.15... element in which a dislocation is located 71 Figure 3.6 Overall response of aluminum matrix for different realizations of random dislocation source and impurity distributions 75 Figure 3.7 Mean overall response of metallic matrix for different values of ρobs with τobs of (a) 0.15 GPa, and (b) 1.20 GPa 75 Figure 3.8 Deformation of RVE and distribution of dislocations within metallic matrix with ρobs =... = 200 MPa 159 Figure 5.14 (a) Mean overall response of pure metallic matrix with different values of εfailure, with the corresponding (b) percentage of voids, and (c) density of active dislocations (i.e excluding annihilated dislocations) in the matrix 163 Figure 5.15 Distribution of dislocations in pure metallic matrix with failure strain εfailure of (a) infinity, (b) 0.020 and (c) 0.015 at γave =... Figure 3.17 Mean overall response of metallic matrix for different RVE sizes with τobs = 0.15 GPa and ρnuc of (a) 160 μm-2, and (b) 40 μm-2 95 Figure 3.18 Deformation of RVE and distribution of dislocations within metallic matrix for RVE sizes of (a) 0.5 μm × 0.5 μm and (b) 3 μm × 3 μm with τobs = 0.15 GPa and ρnuc = 160 μm-2 97 Figure 3.19 Mean overall response of metallic matrix for different RVE sizes... and ρnuc of (a) 40 μm-2, and (b) 160 μm-2 98 Figure 3.20 Mean overall response of metallic matrix for different RVE sizes with τobs = 0.15 GPa and ρobs of (a) 80 μm-2, and (b) 320 μm-2 99 Figure 3.21 Mean overall response of metallic matrix for different RVE sizes with τobs = 0.60 GPa and ρobs of (a) 80 μm-2, and (b) 320 μm-2 100 Figure 3.22 Overall response for various realizations of metallic matrix. .. size of 25 nm for different values of εfailure, with the corresponding (b) percentage of voids, and (c) density of active dislocations (i.e excluding annihilated dislocations) in the matrix 175 Figure 5.21 Distribution of dislocations in composite material with regular rectangular arrangement of 2 per cent inclusion volume fraction and inclusion size of 25 nm at γave = 1.05% for different values of εfailure... nucleation strength of dislocation sources τobs Strength of obstacles or impurities xxiii This page is intentionally left blank xxiv 1.0 Introduction This chapter gives an overview of metal matrix composites (MMCs) and metal matrix nanocomposites (MMNCs), the numerical studies which have been performed on these materials, the size effects observed in MMNCs, and the numerical methods for modelling dislocations... values, and (b) maximum number of subdivisions for adaptive quadrature used to evaluate boundary nodal forces due to dislocation stress field 62 Figure 3.3 Overall response for composite material with different cutoff distances rcut used in the evaluation of dislocation fields 64 Figure 3.4 Distribution of dislocations within matrix of composite material for cut-off distance rcut of (a) 100 b, and (b) 10000... 5.16 (a) Mean overall response of composite material with nonclustered random arrangements of 2 per cent inclusion volume fraction and inclusion size of 25 nm for different values of εfailure, with the corresponding (b) percentage of voids, and (c) density of active dislocations (excluding annihilated dislocations) in the matrix 167 xviii Figure 5.17 Distribution of dislocations in composite material . NUMERICAL STUDY OF METAL MATRIX NANOCOMPOSITES USING DISCRETE DISLOCATION APPROACH ELLIOT LAW NATIONAL UNIVERSITY OF SINGAPORE 2011 NUMERICAL STUDY OF METAL. of inclusions 124 4.5 Material properties of constituent phases 132 4.6 Nature of dislocation pile-ups in metallic nanocomposites 137 vi 5.0 Simulation of damage in metal matrix nanocomposites. response of metallic matrix for different values of ρ obs with τ obs of (a) 0.15 GPa, and (b) 1.20 GPa. 75 Figure 3.8 Deformation of RVE and distribution of dislocations within metallic matrix