含冷却孔镍基合金次级取向效应的应变梯度晶体塑性有限元研究

单晶镍基合金具有优异的耐高温、高强、高韧等性能, 这些力学性能受制造过程引入的次级取向和冷却孔的影响. 已有研究大多关注单孔薄板的变形机理和力学性能, 而工程中应用的往往是多孔薄板, 当前亟需阐明多孔的塑性滑移带变形机理、次级取向效应以及冷却孔引起的应变梯度效应. 文章采用基于位错机制的非局部晶体塑性本构模型对含冷却孔镍基单晶薄板的单拉变形进行了数值模拟. 此模型基于塑性滑移梯度与几何必需位错的关系引入了位错流动项, 因此可有效刻画非均匀变形过程中的应变梯度效应. 为了全面揭示含孔镍基薄板的次级取向效应, 系统研究了[100]和[110]取向(两种次级取向)下镍基薄板的单拉变形行为, 并重点探究了在两种次级取向下冷却孔数量对薄板塑性行为的影响. 此外, 还分析了镍基合金板变形过程中各个滑移系上分切应力变化、主导滑移系开动以及几何必需位错密度的演化过程, 并讨论了塑性滑移量及其分布特征对不同次级取向镍基合金板强度的影响. 研究表明, 单孔和多孔的[110]薄板抗拉强度均低于[100]薄板, 多孔薄板的塑性变形过程比单孔薄板更为复杂且受次级取向影响更大, 并且发生滑移梯度位置主要位于冷却孔附近以及塑性滑移带区域. 研究结果可为工程中镍基合金的设计和服役提供理论指导.      

Distribution based model for the grain boundaries in polycrystalline plasticity

This contribution focuses on the development of constitutive models for the grain boundary region between two crystals, relying on the dislocation based polycrystalline model documented in (Evers, L.P., Parks, D.M., Brekelmans, W.A.M., Geers, M.G.D., 2002. Crystal plasticity model with enhanced hardening by geometrically necessary dislocation accumulation. J. Mech. Phys. Solids 50, 2403–2424; Evers, L.P., Brekelmans, W.A.M., Geers, M.G.D., 2004a. Non-local crystal plasticity model with intrinsic SSD and GND effects. J. Mech. Phys. Solids 52, 2379–2401; Evers, L.P., Brekelmans, W.A.M., Geers, M.G.D., 2004b. Scale dependent crystal plasticity framework with dislocation density and grain boundary effects. Int. J. Solids Struct. 41, 5209–5230). The grain boundary is first viewed as a geometrical surface endowed with its own fields, which are treated here as distributions from a mathematical point of view. Regular and singular dislocation tensors are introduced, defining the grain equilibrium, both in the grain core and at the boundary of both grains. Balance equations for the grain core and grain boundary are derived, that involve the dislocation density distribution tensor, in both its regular and singular contributions. The driving force for the motion of the geometrically necessary dislocations is identified from the pull-back to the lattice configuration of the quasi-static balance of momentum, that reveals the duality between the stress and the curl of the elastic gradient. Criteria that govern the flow of mobile geometrically necessary dislocations (GNDs) through the grain boundary are next elaborated on these bases. Specifically, the sign of the projection of a lattice microtraction on the glide velocity defines a necessary condition for the transmission of incoming GNDs, thereby rendering the set of active slip systems for the glide of outgoing dislocations. Viewing the grain boundary as adjacent bands in each grain with a constant GND density in each, the driving force for the grain boundary slip is further expressed in terms of the GND densities and the differently oriented slip systems in each grain. A semi-analytical solution is developed in the case of symmetrical slip in a bicrystal under plane strain conditions. It is shown that the transmission of plastic slip occurs when the angle made by the slip direction relative to the grain boundary normal is less than a critical value, depending on the ratio of the GND densities and the orientation of the transmitted dislocations.     

An experimental study on grain deformation and interactions in an Al-0.5%Mg multicrystal

Heterogeneous plastic deformation behavior of a coarse-grained Al-0.5%Mg multicrystal was investigated experimentally at the individual grain level. A flat uniaxial tensile specimen consisting of a single layer of millimeter-sized grains was deformed quasi-statically up to an axial strain of 15% at room temperature. The initial local crystallographic orientations of the grains and their evolutions after 5, 12, and 15% plastic strains were measured by electron backscattered diffraction pattern analysis in a scanning electron microscope. The local in-plane plastic strains and rigid body rotations of the grains were measured by correlation of digital optical video images of the specimen surface acquired during the tensile test. It is found that both intergranular and intragranular plastic deformation fields in the aluminum multicrystal specimen under uniaxial tension are highly heterogeneous. Single or double sets of slip-plane traces were predominantly observed on the electro-polished surfaces of the millimeter-sized grains after deformation. The active slip systems associated with these observed slip-plane traces were identified based on the grain orientation after deformation, the Schmid factor, and grain interactions in terms of the slip-plane trace morphology at grain boundaries. It is found that the aluminum multicrystal obeys neither the Sachs nor the Taylor polycrystal deformation models but deforms heterogeneously to favor easy slip transmission and accommodation among the grains.     

Computer simulation of tensile deformation of aluminium bicrystalline specimen by means of elasto-plastic finite element method of crystal

Based on the experimental data of strain-stress relation of single crystal and the complexity of polycrystal deformation, a calculating model is proposed and the calculating program is designed. Finally, the tensile deformation of the bicrystal aluminium specimen is simulated by means of the crystal elasto-plastic finite element method. The strain-stress curve of tension, the details of the slip deformation and the stress distribution are obtained. So the effect of the region near grain boundary is shown obviously.     

梯度纳米结构金属力学性能、变形机理和多尺度计算研究进展

近年来，梯度纳米结构金属因其优越的力学性能和独特的塑性变形机理受到广泛关注，已成为材料与力学学科的热点和前沿。本文首先介绍梯度纳米结构金属的强度、塑性、加工硬化和抗疲劳等核心力学性能，以及晶粒长大、塑性应变梯度和几何必需位错等塑性变形机理及其力学研究。其次介绍梯度纳米结构金属的多尺度计算与模拟研究。最后讨论梯度纳米结构金属研究领域存在的挑战。     

Non-equilibrium grain boundary structure and inelastic deformation using atomistic simulations

Grain boundary influence on material properties becomes increasingly significant as grain size is reduced to the nanoscale. Nanostructured materials produced by severe plastic deformation techniques often contain a higher percentage of high-angle grain boundaries in a non-equilibrium or energetically metastable state. Differences in the mechanical behavior and observed deformation mechanisms are common due to deviations in grain boundary structure. Fundamental interfacial attributes such as atomic mobility and energy are affected due to a higher non-equilibrium state, which in turn affects deformation response. In this research, atomistic simulations employing a biased Monte Carlo method are used to approximate representative non-equilibrium bicrystalline grain boundaries based on an embedded atom method potential, leveraging the concept of excess free volume. An advantage of this approach is that non-equilibrium boundaries can be instantiated without the need of simulating numerous defect/grain boundary interactions. Differences in grain boundary structure and deformation response are investigated as a function of non-equilibrium state using Molecular Dynamics. A detailed comparison between copper and aluminum bicrystals is provided with regard to boundary strength, observed deformation mechanisms, and stress-assisted free volume evolution during both tensile and shear simulations.     

Strain hardening in bent copper foils

A series of systematic tensile and microbend tests were conducted on copper foil specimens with different thicknesses. The specimens were made of a copper foil having almost unidirectional crystal orientations that was considered to be nearly single-crystal. In order to investigate the effects of slip system interactions, two different crystal orientations relative to the tensile direction were considered in the tests: one is close to coplanar double-slip orientation, and the other is close to the ideal cube orientation (the tensile direction nearly coincides to [0 0 1]) that yields multi-planar multi-slip deformation. We extended the microbend test method to include the reversal of bending, and we attempted to divide the total amount of strain-hardening into isotropic and kinematic hardening components. In the tensile tests, no systematic tendency of size dependence was observed. In the microbend tests, size-dependent kinematic hardening behavior was observed for both the crystal orientations, while size dependence of isotropic hardening was observed only for the multi-planar multi-slip case. We introduce an extended crystal plasticity model that accounts for the effects of the geometrically necessary dislocations (GNDs), which correspond to the spatial gradients of crystallographic slips. Through numerical simulations performed using the model, the origin of the size-dependent behavior observed in the microbend tests is discussed.     

Crystal plasticity model with enhanced hardening by geometrically necessary dislocation accumulation

A strain gradient dependent crystal plasticity approach is used to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. Material points are considered as aggregates of grains, subdivided into several fictitious grain fractions: a single crystal volume element stands for the grain interior whereas grain boundaries are represented by bi-crystal volume elements, each having the crystallographic lattice orientations of its adjacent crystals. A relaxed Taylor-like interaction law is used for the transition from the local to the global scale. It is relaxed with respect to the bi-crystals, providing compatibility and stress equilibrium at their internal interface. During loading, the bi-crystal boundaries deform dissimilar to the associated grain interior. Arising from this heterogeneity, a geometrically necessary dislocation (GND) density can be computed, which is required to restore compatibility of the crystallographic lattice. This effect provides a physically based method to account for the additional hardening as introduced by the GNDs, the magnitude of which is related to the grain size. Hence, a scale-dependent response is obtained, for which the numerical simulations predict a mechanical behaviour corresponding to the Hall-Petch effect. Compared to a full-scale finite element model reported in the literature, the present polycrystalline crystal plasticity model is of equal quality yet much more efficient from a computational point of view for simulating uniaxial tension experiments with various grain sizes.     

Predicting the Hall–Petch effect in fcc metals using non-local crystal plasticity

Many conventional continuum approaches to solid mechanics do not address the size sensitivity of deformation to microstructural features like grain boundaries, and are therefore unable to capture much of the experimentally observed behavior of polycrystal deformation. We propose a non-local crystal plasticity model, in which the geometrically necessary dislocation (GND) density is calculated using a non-local integral approach. The model is based on augmented F^eF^p kinematics, which account for the initial microstructure (primarily grain boundaries) present in the polycrystal. With the augmented kinematics, the initial GND and the evolving GND state are determined in a consistent manner. The expanded kinematics and the non-local crystal plasticity model are used to simulate the tensile behavior in copper polycrystals with different grain sizes ranging from 14 μm to 244 μm. The simulation results show a grain size dependence on the polycrystal’s yield strength, which are in good agreement with the experimental data.     

Grain-boundary sliding and separation in polycrystalline metals: application to nanocrystalline fcc metals

In order to model the effects of grain boundaries in polycrystalline materials we have coupled a crystal-plasticity model for the grain interiors with a new elastic-plastic grain-boundary interface model which accounts for both reversible elastic, as well irreversible inelastic sliding-separation deformations at the grain boundaries prior to failure. We have used this new computational capability to study the deformation and fracture response of nanocrystalline nickel. The results from the simulations reflect the macroscopic experimentally observed tensile stress-strain curves, and the dominant microstructural fracture mechanisms in this material. The macroscopically observed nonlinearity in the stress-strain response is mainly due to the inelastic response of the grain boundaries. Plastic deformation in the interior of the grains prior to the formation of grain-boundary cracks was rarely observed. The stress concentrations at the tips of the distributed grain-boundary cracks, and at grain-boundary triple junctions, cause a limited amount of plastic deformation in the high-strength grain interiors. The competition of grain-boundary deformation with that in the grain interiors determines the observed macroscopic stress-strain response, and the overall ductility. In nanocrystalline nickel, the high-yield strength of the grain interiors and relatively weaker grain-boundary interfaces account for the low ductility of this material in tension.     

Dislocation nucleation from bicrystal interfaces with dissociated structure

Atomistic calculations are used to model the nucleation of partial dislocations during a tensile deformation from bicrystal interfaces with dissociated structure. Interfaces with this type of structure occur primarily in materials with low intrinsic stacking fault energies. In this work, the initial structure of each bicrystal interface is refined using energy minimization techniques. Molecular dynamics simulations are then used to study the deformation of each interface in uniaxial tension perpendicular to the boundary plane at a constant strain rate. Analysis focuses on the evolution of the dissociated interface structure prior to the dislocation nucleation event and the resulting structure of the boundary after the emission of partial dislocations from the interface. Dislocation nucleation occurs predominantly at the dissociated interface structural unit, while the spacing between interface features is identified as an important length scale that affects the failure mode. The evolution of the dissociated interface structure and the nucleation of partial dislocations are found to be similar to results obtained in a previous atomistic study of the stress dependence of a lock formation containing a stair-rod dislocation.     

Numerical study of the plastic behaviour in tension of welds in high strength steels

The influence of the mismatch between material properties and constraint on the plastic deformation behaviour of the heat affected zone of welds in high strength steels is investigated in this study, using finite element simulations. An elastoplastic implicit three-dimensional finite element code (EPIM3D) was used in the analysis. The paper presents the mechanical model of the code and the methodology used for the numerical simulation of the tensile test of welded joints. Numerical results of the tensile test of welded samples with different hypothetical widths for the Heat Affected Zone and various material mismatch levels are shown. The analysis concerns the overall strength and ductility of the joint and in relation to the plastic behaviour of the heat affected zone. The influence of the yield stress, tensile strength and constraint on the stress and plastic strain distribution in the soft heat affected zone is also discussed.     

Mechanism-based strain gradient crystal plasticity—I. Theory

We have been developing the theory of mechanism-based strain gradient plasticity (MSG) to model size-dependent plastic deformation at micron and submicron length scales. The core idea has been to incorporate the concept of geometrically necessary dislocations into the continuum plastic constitutive laws via the Taylor hardening relation. Here we extend this effort to develop a mechanism-based strain gradient theory of crystal plasticity. In this theory, an effective density of geometrically necessary dislocations for a specific slip plane is introduced via a continuum analog of the Peach-Koehler force in dislocation theory and is incorporated into the plastic constitutive laws via the Taylor relation.     

CrMnFeCoNi纳米晶温度相关的拉伸行为研究

摘要：高熵合金是一种由多种主元元素组成的新型合金。实验研究表明等原子比CrMnFeCoNi高熵合金在低温下具有比室温更高的拉伸强度和断裂韧性。本文针对这一现象，利用分子动力学模拟对平均晶粒尺寸为6 nm的CrMnFeCoNi纳米晶在300、200和77 K下分别进行拉伸模拟。模拟研究揭示了纳米尺度CrMnFeCoNi高熵合金力学行为的温度效应和强韧机理。微结构演化分析表明：低温下，塑性变形阶段，滑移系开动的较少，位错滑移所受的阻力越大，屈服强度和抗拉强度越大；模型破坏时，孔洞缺陷形核较慢，更多孔洞缺陷演化成断口，更多的断口分摊拉伸应变，使得高熵合金纳米晶的低温韧性更好。     

双相不锈钢循环微变形行为的原位试验研究

对双相不锈钢开展拉伸和应力循环控制下的原位试验研究,选择特征区域对其局部微变形行为的产生、分布和演化进行了研究和讨论．结果表明：双相不锈钢在较低的载荷作用下内部出现明显的应变梯度,随着拉伸载荷或循环周次的增长,其局部微应变不断增加,最大应变在奥氏体相、小晶粒、狭窄晶粒和不平整的晶界区不断汇集,同一时刻的局部最大微应变远大于宏观应变;循环载荷作用下双相不锈钢在微观尺度产生了局部微棘轮变形效应,微棘轮应变随着循环周次的增加不断累积,其增长速度逐渐降低并趋于饱和;随着循环周次的增加,最大微棘轮变形区逐渐扩展、联通形成大塑性变形带,晶粒走向呈45°方向时,变形带贯穿两相向垂直加载方向发展,而晶粒呈树根状走向时,由于受到两相交互和铁素体的阻碍作用,微棘轮变形带主要集中在奥氏体中沿着奥氏体晶粒走向发展．     

Discrete dislocation dynamics simulations to interpret plasticity size and surface effects in freestanding FCC thin films

Strong size effects have been experimentally observed when microstructural features approach the geometric dimensions of the sample. In this work experimental investigations and discrete dislocation analyses of plastic deformation in metallic thin films have been performed. Columnar grains representative of the film microstructure are here considered. Simulations are based on the assumptions that sources are scarcely available in geometrically confined systems and nucleation sites are mainly located at grain boundaries. Especially, we investigated the influence on the mesoscopic constitutive response of the two characteristic length scales, i.e., film thickness and grain size. The simulated plastic response qualitatively reproduces the experimentally observed size effects while the main deformation mechanisms appear to be in agreement with TEM analyses of tested samples. A new interpretation of size scale plasticity is here proposed based on the probability of activating grain boundary dislocation sources. Moreover, the key role of a parameter such as the grain aspect ratio is highlighted. Finally, the unloading behavior has been investigated and a strong size dependent Bauschinger effect has been found. An interpretation of these phenomena is proposed based on the analysis of the back stress distribution within the samples.     

Geometrically necessary dislocations in viscoplastic single crystals and bicrystals undergoing small deformations

In this study we develop a gradient theory of small-deformation single-crystal plasticity that accounts for geometrically necessary dislocations (GNDs). The resulting framework is used to discuss grain boundaries. The grains are allowed to slip along the interface, but growth phenomenona and phase transitions are neglected. The bulk theory is based on the introduction of a microforce balance for each slip system and includes a defect energy depending on a suitable measure of GNDs. The microforce balances are shown to be equivalent to nonlocal yield conditions for the individual slip systems, yield conditions that feature backstresses resulting from energy stored in dislocations. When applied to a grain boundary the theory leads to concomitant yield conditions: relative slip of the grains is activated when the shear stress reaches a suitable threshold; plastic slip in bulk at the grain boundary is activated only when the local density of GNDs reaches an assigned threshold. Consequently, in the initial stages of plastic deformation the grain boundary acts as a barrier to plastic slip, while in later stages the interface acts as a source or sink for dislocations. We obtain an exact solution for a simple problem in plane strain involving a semi-infinite compressed specimen that abuts a rigid material. We view this problem as an approximation to a situation involving a grain boundary between a grain with slip systems aligned for easy flow and a grain whose slip system alignment severely inhibits flow. The solution exhibits large slip gradients within a thin layer at the grain boundary.     

Three-dimensional local stress analysis on grain boundaries in polycrystalline material

In order to understand the initiation behavior of microstructurally small cracks in a stress corrosion cracking condition, it is important to know the tensile normal stress acting on the grain boundary (normal GB stress). The local stress in a polycrystalline body is enhanced by the inhomogeneity which stems from the shape and orientation of each grain. The stress in a three-dimensional polycrystalline body consisting of 100 grains with random orientation, under a remote uniform tensile stress condition, is evaluated by the finite element method. It was revealed that the local stress on the polycrystalline body is inhomogeneous under uniform applied stress and becomes large at those grain boundaries that are perpendicular to the load axis, though there is large fluctuation. It was also shown that the normal GB stress tends to be large near the triple points due to the deformation constraint caused by adjacent grains. Finally, the maximum stress on the surface of a large component caused by the inhomogeneity was evaluated by using Gumbel statistics.     

An interfacial energy incorporated couple stress crystal plasticity and the finite element simulation of grain subdivision

A couple stress crystal plasticity formulation that incorporates interfacial couple stress energy was proposed in terms of the virtual work-rate principle for finite element method. By applying the assumed constitutive models of couple stress at the grain boundary as well as the grain interior, finite element simulations were conducted for various crystal models, with different grain subdivision models to examine how plastic deformation work is affected by grain subdivision from the interfacial couple stress energy effect.Finite element simulation results showed that the amount of predicted plastic deformation work depends on grain subdivision, and that the amount of work can be minimized for a particular grain subdivision. We inferred from the simulation results that actual grain subdivision might correspond to the minimum amount of plastic deformation work and, if this correlation is validated, actual grain subdivision might be predicted based on the interfacial energy incorporated couple stress crystal plasticity.