离散位错动力学算法及其在材料塑性行为模拟中的应用  

Discrete dislocation dynamics algorithms and their application in modeling of plastic behaviors of crystalline materials

作　　者：黄敏生[1,2] 黄嵩 梁爽[1,2] 李振环 Minsheng Huang;Song Huang;Shuang Liang;Zhenhuan Li(Department o f Mechanics, Huazhong University of Science & Technology, Wuhan 430074, China;Hubei Key Laboratory of Engineering Structural Analysis and Safety Assessment, Huazhong University of Science & Technology, Wuhan 430074,China)

机构地区：[1]华中科技大学力学系,武汉430074 [2]华中科技大学工程结构分析与安全评定湖北省重点实验室,武汉430074

出　　处：《科学通报》2019年第18期1864-1877,共14页Chinese Science Bulletin

基　　金：国家自然科学基金(11572136,11632007);国家优秀青年科学基金(11822204)资助

摘　　要：晶体材料的塑性变形由位错的运动演化而引起.离散位错动力学(discrete dislocation dynamics, DDD)通过直接模拟大量位错的演化而研究材料的塑性变形,因此能够揭示材料微结构-位错微结构-塑性力学行为之间内在的物理关联,并能够自然而然地捕捉塑性变形微米/亚微米特征尺度下本征的尺度效应.它所能模拟的尺度介于微观分子动力学模拟和宏观有限元模拟之间,在多尺度算法中起到承上启下的作用.本文首先系统地发展、完善和丰富了离散位错动力学-有限元(finite element method, FEM)叠加算法、DDD-FEM直接耦合算法(discrete-continuous method, DCM)以及离散位错动力学-扩展有限元(extended finite element method, XFEM)耦合算法等框架体系.在此基础上,利用这些方法对单晶镍基高温合金的塑性变形机理、晶体材料的断裂和损伤变形行为以及塑性行为的微尺度和微结构效应3个方面开展了系统的研究.所得模拟结果指导了基于微结构和位错机制的单晶镍基高温合金晶体塑性本构模型的建立,丰富和加深了人们对材料强化、循环塑性、断裂、损伤、尺度效应和微结构效应的认识.此外,离散位错动力学可进一步应用于诸如高温、高压、高应变率、化学腐蚀环境、高辐照等极端条件下晶体材料塑性行为的研究,是材料力学行为多尺度模拟研究中的重要一环.The plastic behaviors of crystalline materials are induced by the evolution of a series of dislocations. Discrete dislocation dynamics(DDD) can be used to model plastic deformation by taking into account the nucleation, annihilation, motion of dislocations, and formation of jogs, junctions, networks, and other dislocation structures directly. As a result, it can be used to investigate the physical correlation between material microstructures, dislocation structures, and plastic responses. It can also model the intrinsic size effect of plastic behaviors at the micron/submicron scales. Since its modeling scale is much larger than that of the microscopic molecular dynamics but is smaller than that of the macroscopic finite element method(FEM), DDD plays a connecting link role in multiscale modeling.In the present manuscript, three categories of DDD schemes, i.e., the DDD-FEM superposition algorithm, the DDD-FEM direct coupling algorithm, and DDD-extended FEM(XFEM) coupling algorithm, are further developed and optimized. In the DDD-FEM superposition algorithm, the consideration of dislocation climb is specifically included to take into account the high-temperature effect. In the DDD-FEM direct coupling algorithm, the so-called virtual dislocations are introduced to distribute the Eigen plastic strain by dislocation evolution to the FEM module more appropriately, and a special scheme is developed to calculate the Peach-Koehler force on the dislocation efficiently. To model the problems with strong or weak interfaces(such as cracks, voids, particles, and grain boundaries), the scheme coupling the DDD and XFEM, which is not sensitive to the finite element mesh, is further developed carefully. Each of these three DDD schemes has its own special characteristics, and so they are suitable for solving different plastic problems.Based on these three kinds of DDD algorithms, the plastic deformation mechanisms of nickel-based superalloy, the damage and fracture behaviors of crystalline materials, and the size effect and microst

关 键 词：离散位错动力学 有限元 扩展有限元 镍基高温合金 微尺度效应 微结构相关性 

分 类 号：TG132.3[一般工业技术—材料科学与工程]