Finite versus small strain discrete dislocation analysis of cantilever bending of single crystals

Plastic size effects in single crystals are investi-gated by using finite strain and small strain discrete dislo-cation plasticity to analyse the response of cantilever beam specimens. Crystals with both one and two active slip sys-tems are analysed, as well as specimens with different beam aspect ratios. Over the range of specimen sizes analysed here, the bending stress versus applied tip displacement response has a strong hardening plastic component. This hardening rate increases with decreasing specimen size. The hardening rates are slightly lower when the finite strain discrete disloca-tion plasticity (DDP) formulation is employed as curving of the slip planes is accounted for in the finite strain formulation. This relaxes the back-stresses in the dislocation pile-ups and thereby reduces the hardening rate. Our calculations show that in line with the pure bending case, the bending stress in cantilever bending displays a plastic size dependence. How-ever, unlike pure bending, the bending flow strength of the larger aspect ratio cantilever beams is appreciably smaller. This is attributed to the fact that for the same applied bend-ing stress, longer beams have lower shear forces acting upon them and this results in a lower density of statistically stored dislocations.     

Investigation on the plastic work-heat conversion coefficient of 7075-T651 aluminum alloy during an impact process based on infrared temperature measurement technology

The plastic work-heat conversion coefficient is one key parameter for studying the work-heat conversion under dynamic deformation of materials. To explore this coefficient of 7075-T651 aluminum alloy under dynamic compression, dynamic compression experiments using the Hopkinson bar under four groups of strain rates were conducted, and the temperature signals were measured by constructing a transient infrared temperature measurement system. According to stress versus strain data as well as the corresponding temperature data obtained through the experiments, the influences of the strain and the strain rate on the coefficient of plastic work converted to heat were analyzed.The experimental results show that the coefficient of plastic work converted to heat of 7075-T651 aluminum alloy is not a constant at the range of 0.85–1 and is closely related to the strain and the strain rate. The change of internal structure of material under high strain rate reduces its energy storage capacity, and makes almost all plastic work convert into heat.     

FRACTURE ANALYSIS OF A FUNCTIONALLY GRADED STRIP UNDER PLANE DEFORMATION

In this paper the plane elasticity problem for a functionally graded strip containing a crack is considered. It is assumed that the reciprocal of the shear modulus is a linear function of the thickness-coordinate, while the Possion's ratio keeps constant. By utilizing the Fourier transformation technique and the transfer matrix method, the mixed boundary problem is reduced to a system of singular integral equations that are solved numerically. The influences of the geometric parameters and the graded parameter on the stress intensity factors and the strain energy release rate are investigated. The numerical results show that the graded parameters, the thickness of the strip and the crack size have significant effects on the stress intensity factors and the strain energy release rate.     

Combined grain size, strain rate and loading condition effects on mechanical behavior of nanocrystalline Cu under high strain rates

Molecular dynamics simulations of nanocrys-talline Cu with average grain sizes of 3.1 nm, 6.2 nm, 12.4 nm and 18.6 nm under uniaxial strain and stress tension at strain rates of 10 8 s 1 , 10 9 s 1 and 10 10 s 1 are performed to study the combined grain size, strain rate and loading condition effects on mechanical properties. It is found that the strength of nanocrystalline Cu increases as grain size increases regardless of loading condition. Both the strength and ductility of nanocrystalline Cu increase with strain rate except that there is no monotonic relation between the strength and strain rate for specimens under uniaxial strain loading. Moreover, the strength and ductility of specimens under uniaxial strain loading are lower than those under uniaxial stress loading. The nucleation of voids at grain boundaries and their subsequent growth characterize the failure of specimens under uniaxial strain loading, while grain boundary sliding and necking dominate the failure of specimens under uniaxial stress loading. The rate dependent strength is mainly caused by the dynamic wave effect that limits dislocation motion, while combined twinning and slipping mechanism makes the material more ductile at higher strain rates.     

Deformation of plastically compressible hardening-softening-hardening solids

Motivated by a model of the response of vertically aligned carbon nanotube (VACNT) pillars in uniaxial compression, we consider the deformation of a class of compressible elastic-viscoplastic solids with a hardening-softening-hardening variation of flow strength with plastic strain. In previous work (Hutchens et al. 2011) a constitutive relation was presented and used to model the response of VACNT pillars in axisymmetric compression. Subsequently, it was found that due to a programming error the constitutive relation presented in the paper (Hutchens et al. 2011) was not the one actually implemented. In particular, the plastic flow rule actually used did not satisfy plastic normality. Here, we present the constitutive formulation actually implemented in the previous work (Hutchens et al. 2011). Dynamic, finite deformation, finite element calculations are carried out for uniaxial compression, uniaxial tension and for indentation of a "half-space" by a conical indenter tip. A sequential buckling-like deformation mode is found in com- pression when there is plastic non-normality and hardening-softening-hardening. The same material characterization gives rise to a Lüders band-like deformation mode in ten- sion. When there is a deformation mode with a sharp front along mesh boundaries, the overall stress-strain response contains high frequency oscillations that are a mesh artifact. The responses of non-softening solids are also analyzed and their overall stress-strain behavior and deformationmodes are compared with those of hardening-softening- hardening solids. We find that indentation with a sharp in- denter tip gives a qualitatively equivalent response for hardening and hardening-softening-hardening solids.     

On the homogenization of metal matrix composites using strain gradient plasticity

The homogenized response of metal matrix composites（MMC） is studied using strain gradient plasticity.The material model employed is a rate independent formulation of energetic strain gradient plasticity at the micro scale and conventional rate independent plasticity at the macro scale. Free energy inside the micro structure is included due to the elastic strains and plastic strain gradients. A unit cell containing a circular elastic fiber is analyzed under macroscopic simple shear in addition to transverse and longitudinal loading. The analyses are carried out under generalized plane strain condition. Micro-macro homogenization is performed observing the Hill-Mandel energy condition,and overall loading is considered such that the homogenized higher order terms vanish. The results highlight the intrinsic size-effects as well as the effect of fiber volume fraction on the overall response curves, plastic strain distributions and homogenized yield surfaces under different loading conditions. It is concluded that composites with smaller reinforcement size have larger initial yield surfaces and furthermore,they exhibit more kinematic hardening.     

CREEP BEHAVIOR OF Zr-1.5Nb-0.4Sn-0.1Fe-0.1Cu ALLOY

Creep behavior of the Zr-1.5Nb-0.4Sn-0.1Fe-0.1Cu alloy sheet is investigated from 300℃ to 400℃ in the stress range from 50 MPa to 180 MPa along the rolling direction. The measured strain rates range from 8.8 × 10^-10 s^-1 to 4.7 × 10^-7 s^-1. The activation energies are estimated to assess the creep deformation mechanisms in this alloy. The strain rate is slightly different at low stress, while it shows a distinct difference at high stresses. Stress exponents of this alloy increase with increasing applied stress at all testing temperatures. It is concluded that the creep deformation of the Zr-1.5Nb-0.4Sn-0.1Fe-0. 1 Cu alloy is controlled by the diffusion creep at low stress region and by the climbing of dislocations at high stress region.     

Bauschinger and size effects in thin-film plasticity due to defect-energy of geometrical necessary dislocations

The Bauschinger and size effects in the thinfilm plasticity theory arising from the defect-energy of geometrically necessary dislocations (GNDs) are analytically investigated in this paper. Firstly, this defect-energy is deduced based on the elastic interactions of coupling dislocations (or pile-ups) moving on the closed neighboring slip plane. This energy is a quadratic function of the GNDs density, and includes an elastic interaction coefficient and an energetic length scale L. By incorporating it into the work- conjugate strain gradient plasticity theory of Gurtin, an energetic stress associated with this defect energy is obtained, which just plays the role of back stress in the kinematic hardening model. Then this back-stress hardening model is used to investigate the Bauschinger and size effects in the tension problem of single crystal Al films with passivation layers. The tension stress in the film shows a reverse dependence on the film thickness h. By comparing it with discrete-dislocation simulation results, the length scale L is determined, which is just several slip plane spacing, and accords well with our physical interpretation for the defect- energy. The Bauschinger effect after unloading is analyzed by combining this back-stress hardening model with a friction model. The effects of film thickness and pre-strain on the reversed plastic strain after unloading are quantified and qualitatively compared with experiment results.     

Size and strain rate effects in tensile strength of penta-twinned Ag nanowires

Penta-twinned Ag nanowires(pt-AgNWs) have recently attracted much attention due to their interesting mechanical and physical properties. Here we perform largescale atomistic simulations to investigate the influence of sample size and strain rate on the tensile strength of pt-AgNWs. The simulation results show an apparent size effect in that the nanowire strength(defined as the critical stress for dislocation nucleation) increases with decreasing wire diameter. To account for such size effect, a theoretical model involving the interaction between an emerging dislocation and the twin boundary has been developed for the surface nucleation of dislocations. It is shown that the model predictions are in quantitative agreement with the results from atomistic simulations and previous experimental studies in the literatures. The simulations also reveal that nanowire strength is strain-rate dependent, which predicts an activation volume for dislocation nucleation in the range of 1–10b~3,where b is the magnitude of the Burgers vector for a full dislocation.     

基于微纳米压入法提取无铅焊料合金弹塑性力学参数的研究

本文通过微纳米压入法结合数值模拟研究了无铅焊料合金SnAg3.5 的弹塑性力学性能，分别采用圆柱形压头及两种不同锥角压头对无铅焊料合金进行压入测试：基于圆柱形压头测试过程中接触面积恒定的特点得到了无铅焊料的弹性模量，进一步采用塑性应变梯度理论对两种锥角压头的测试结果予以修正并通过数值模拟反分析得到相应的特征应力值，同时基于压入特征塑性应变与压头锥角的关系式得到两种不同锥角压头下的特征应变值，在此基础上经求解方程组得到焊料合金的初始屈服应力与应变强化因子，进而得到了焊料合金的幂强化弹塑性本构关系．该方法剔除了压入尺度效应的影响并保证了所得本构关系的唯一性，给出了一种通过原位压入测试表征金属材料弹塑性力学性能的有效方法．     

A study on robust indentation techniques to evaluate elastic–plastic properties of metals

Spherical indentation is studied based on numerical analysis and experiment, to develop robust testing techniques to evaluate isotropic elastic–plastic material properties of metals. The representative stress and plastic strain concept is critically investigated via finite element analysis, and some conditions for the representative values are suggested. The representative values should also be a function of material properties, not only indenter angle for sharp indenter and indentation depth for spherical indenter. The pros and cons of shallow and deep spherical indentation techniques are also discussed. For an indentation depth of 20% of an indenter diameter, the relationships between normalized indentation parameters and load–depth data are characterized, and then numerical algorithm to estimate material elastic–plastic curve is presented. From the indentation load–depth curve, the new approach provides stress–strain curve and the values of elastic modulus, yield strength, and strain-hardening exponent with an average error of less than 5%. The method is confirmed to be valid for various elastic properties of indenter. Experimental validation of the approach then is performed by using developed micro-indentation system. For the material severely disobeying power law hardening, a modified method to reduce errors of predicted material properties is contrived. It is found that our method is robust enough to get ideal power law properties, and applicable to input of more complex physics.     

DETERMINATION OF CREEP PARAMETERS FROM INDENTATION CREEP EXPERIMENTS

ＩｎｔｒｏｄｕｃｔｉｏｎＣｒｅｅｐｉｓｔｈｅｔｉｍｅｄｅｐｅｎｄｅｎｔｐｌａｓｔｉｃｄｅｆｏｒｍａｔｉｏｎｏｆｍａｔｅｒｉａｌｓｗｈｉｃｈｉｓｃｈａｒａｃｔｅｒｉｚｅｄｂｙａｓｔｒｏｎｇｄｅｐｅｎｄｅｎｃｅｏｆｔｈｅｃｒｅｅｐｒａｔｅｆｒｏｍｓｔｒｅｓｓσｏｎｃｏｎｓｔａｎｔｔｅｍｐｅｒａｔｕｒｅＴ .Ｔｈｅｓｅｃｏｎｄａｒｙｃｒｅｅｐｒａｔｅ εｄｅｐｅｎｄｓｏｎｓｔｒｅｓｓσａｓ ε=Ｃｃ·σｎＣ, ( 1 )ｗｈｅｒｅＣｃｉｓａｃｏｎｓｔａｎｔ,ｎＣｉｓｔｈｅｃｒｅｅｐｓｔｒｅｓｓｅｘｐｏｎｅｎｔ.Ｔｈ…     

An expanding cavity model incorporating strain-hardening and indentation size effects

An expanding cavity model (ECM) for determining indentation hardness of elastic strain-hardening plastic materials is developed. The derivation is based on a strain gradient plasticity solution for an internally pressurized thick-walled spherical shell of an elastic power-law hardening material. Closed-form formulas are provided for both conical and spherical indentations. The indentation radius enters these formulas with its own dimensional identity, unlike that in classical plasticity based ECMs where indentation geometrical parameters appear only in non-dimensional forms. As a result, the newly developed ECM can capture the indentation size effect. The formulas explicitly show that indentation hardness depends on Young’s modulus, yield stress, strain-hardening exponent and strain gradient coefficient of the indented material as well as on the geometry of the indenter. The new model reduces to existing classical plasticity based ECMs (including Johnson’s ECM for elastic–perfectly plastic materials) when the strain gradient effect is not considered. The numerical results obtained using the newly developed model reveal that the hardness is indeed indentation size dependent when the indentation radius is very small: the smaller the indentation, the larger the hardness. Also, the indentation hardness is seen to increase with the Young’s modulus and strain-hardening level of the indented material for both conical and spherical indentations. The strain-hardening effect on the hardness is observed to be significant for materials having strong strain-hardening characteristics. In addition, it is found that the indentation hardness increases with decreasing cone angle of the conical indenter or decreasing radius of the spherical indenter. These trends agree with existing experimental observations and model predictions.     

Mechanics of indentation of plastically graded materials—I: Analysis

The introduction of controlled gradients in plastic properties is known to influence the resistance to damage and cracking at contact surfaces in many tribological applications. In order to assess potentially beneficial effects of plastic property gradients in tribological applications, it is essential first to develop a comprehensive and quantitative understanding of the effects of yield strength and strain hardening exponent on contact deformation under the most fundamental contact condition: normal indentation. To date, however, systematic and quantitative studies of plasticity gradient effects on indentation response have not been completed. A comprehensive parametric study of the mechanics of normal indentation of plastically graded materials was therefore undertaken in this work by recourse to finite element method (FEM) computations. On the basis of a large number of computational simulations, a general methodology for assessing instrumented indentation response of plastically graded materials is formulated so that quantitative interpretations of depth-sensing indentation experiments could be performed. The specific case of linear variation in yield strength with depth below the indented surface is explored in detail. Universal dimensionless functions are extracted from FEM simulations so as to predict the indentation load versus depth of penetration curves for a wide variety of plastically graded engineering metals and alloys for interpretation of, and comparisons with, experimental results. Furthermore, the effect of plasticity gradient on the residual indentation pile-up profile is systematically studied. The computations reveal that pile-up of the graded alloy around the indenter, for indentation with increasing yield strength beneath the surface, is noticeably higher than that for the two homogeneous reference materials that constitute the bounding conditions for the graded material. Pile-up is also found to be an increasing function of yield strength gradient and a decreasing function of frictional coefficient. The stress and plastic strain distributions under the indenter tip with and without plasticity gradient are also examined to rationalize the predicted trends. In Part II of this paper, we compare the predictions of depth-sensing indentation and pile-up response with experiments on a specially made, graded model Ni-W alloy with controlled gradients in nanocrystalline grain size.     

考虑压头曲率半径和应变梯度的微压痕分析

在压头尖端曲率半径取100nm的前提下，采用Chen和Wang的应变梯度理论，对微压痕实验进行了系统的数值分析. 首先通过拟合载荷-位移实验曲线的后半段来确定材料的屈服应力和幂硬化指数值,然后用有限元方法数值模拟压痕实验，并将计算得到的整段载荷-位移曲线及硬度-位移曲线和实验结果进行了比较. 结果表明应变梯度理论所预测的计算结果和实验结果很好地符合,包括压痕深度在亚微米和微米范围内的整段曲线.     

The equivalent axisymmetric model for Berkovich indenters in power-law hardening materials

Nix and Gio [Nix, W.D., Gao, H.J., 1998. Indentation size effects in crystalline materials: a law for strain gradient plasticity. Journal of the Mechanics and Physics of Solids 46, 411–425] established an important relation between the micro-indentation hardness and indentation depth for axisymmetric indenters. For the Berkovich indenter, however, this relation requires an equivalent cone angle. Qin et al. [Qin, J., Huang, Y., Xiao, J., Hwang, K.C., 2009. The equivalence of axisymmetric indentation model for three-dimensional indentation hardness. Journal of Materials Research 24, 776–783] showed that the widely used equivalent cone angle from the criterion of equal base area leads to significant errors in micro-indentation, and proposed a new equivalence of equal cone angle for iridium. It is shown in this paper that this new equivalence holds for a wide range of plastic work hardening materials. In addition, the prior equal-base-area criterion does not hold because the Berkovich indenter gives much higher density of geometrically necessary dislocations than axisymmetric indenter. The equivalence of equal cone angle, however, does not hold for Vickers indenter.     

纳米压痕试验中压头尺寸效应的准连续介质法分析

采用准连续介质法模拟了单晶铝纳米压痕试验过程,分析了不同宽度的刚性矩形压头所引起的初始塑性变形特点,获得了载荷-压深、应变能-位移和硬度-压深曲线.从位错理论的角度分析了压头尺寸对纳米压痕测试结果的影响.研究发现:随着压头宽度的不断增大,压头下方位错形核所需要的载荷和压深程度增大,需要的应变能增加,应变能的变化速率递增,纳米硬度值减小,呈现出明显的尺寸效应.同时表明在一定的压入深度下,硬度与压头尺寸之间存在着一定的比例关系,不同尺寸压头获得的硬度值可以相互换算,但当矩形刚性压头宽度大于或等于120时这种尺寸效应消失.研究结果为纳米压痕实验过程中压头尺寸的选择提供了参考依据.     

纳米压痕过程的三维有限元数值试验研究

采用有限元方法模拟了纳米压痕仪的加、卸载过程，三维有限元模型考虑了纳米压痕仪的标准Berkovich压头．介绍了有限元模型的几何参数、边界条件、材料特性与加载方式，讨论了摩擦、滑动机制、试件模型的大小对计算结果的影响，进行了计算结果与标准试样实验结果的比较，证实了模拟的可靠性．在此基础上，重点研究了压头尖端曲率半径对纳米压痕实验数据的影响．对比分析了尖端曲率半径r＝0与r＝100nm两种压头的材料压痕载荷—位移曲线．结果表明，当压头尖端曲率半径r≠0时，基于经典的均匀连续介质力学本构理论、传统的实验手段与数据处理方法，压痕硬度值会随着压痕深度的减小而升高．     

蠕变压痕实验的计算机模拟

采用有限元的方法对双相材料的蠕变压痕实验进行了数值模拟，在有限元数值解的基础上，定义了相应于传统单轴蠕变实验的“等效应力”和“等效应变”，建立了蠕变压痕实验同传统单轴蠕变实验之间的关系，给出了确定薄膜蠕变应力指数和蠕变常数的方法；同时数值解的结果表明，实验中通过控制压痕深度不超过薄膜厚度的 ５％～１０％，忽略基体的硬化指数对确定薄膜性能的影响存在一定的误差，但基体的弹性模量对确定薄膜的蠕变性能影响不大。