Constitutive model for plastic deformation of nanocrystalline materials with shear band

A phase mixture model was used to study the plastic deformation behaviors in hardening stage of nanocrystalline materials. The material was considered as a composite of grain interior phase and grain boundary (GB) phase. The constitutive equations of the two phases were determined in term of their main deformation mechanisms. In softening stage, a shear band deformation mechanism was presented and the corresponding constitutive relation was established. Numerical simulations have shown that the predications fit well with experimental data. The investigation using the finite-element method (FEM) provided a direct insight into quantifying shear localization effect in nanocrystalline materials.     

On criteria for dynamic adiabatic shear band propagation

In this work, we postulate the physical criterion for dynamic shear band propagation, and based on this assumption, we implement a numerical algorithm and a computation criterion to simulate initiation and propagation of dynamic adiabatic shear bands (ASBs). The physical criterion is based on the hypothesis that material inside the shear band region undergoes a dynamic recrystallization process during deformation under high temperature and high strain-rate conditions. In addition to providing a new perspective to the physics of the adiabatic shearbanding process and identifying material properties that play a crucial role in defining the material's susceptibility to ASBs, the proposed criterion is instrumental in numerical simulations of the propagation of ASBs when multi-physics models are adopted to describe and predict the complex constitutive behavior of ASBs in ductile materials. Systematic and large scale meshfree simulations have been conducted to test and validate the proposed criterion by examining the formation, propagation, and post-bifurcation behaviors of ASBs in two materials, 4340 steel and OFHC copper. The effects of heat conduction, in particular the length scale introduced by heat conduction, are also studied. The results of the numerical simulations are compared with experimental observations and a close agreement is found for various characteristic features of ASBs, such as the shear band width, speed of propagation, and maximum temperature.     

Shear band formation in thermo-viscoplastic solids under plane strain compression

The plane strain compression of a rectangular block is numerically investigated for the study of dynamic shear band development in thermo-elasto-viscoplastic materials from an internal inhomogeneity. As expected, it plays an important role in triggering the onset of shear, localization as well as thermal softening. And the competition between the strain, strain-rate hardening and thermal softening exists throughout the process. It is found that shear band develops at a 45-degree angle to the compression axis. In the light of given patterns of deformation and temperature, shear band evolution accelerated by thermal softening is retarded by the inertial effects. Interestingly, a similar temperature band is also formed along the trajectory of the localized deformation band. The calculations also show the energy evolution during the coupled thermo-mechanical process of shear band propagation. Finally, the mesh effect is discussed in terms of the numerical results from two different meshes. The project is supported by the National Natural Sciences Foundation of China.     

Multiresolution continuum modeling of micro-void assisted dynamic adiabatic shear band propagation

A thermal-mechanical multiresolution continuum theory is applied within a finite element framework to model the initiation and propagation of dynamic shear bands in a steel alloy. The shear instability and subsequent stress collapse, which are responsible for dynamic adiabatic shear band propagation, are captured by including the effects of shear driven microvoid damage in a single constitutive model. The shear band width during propagation is controlled via a combination of thermal conductance and an embedded evolving length scale parameter present in the multiresolution continuum formulation. In particular, as the material reaches a shear instability and begins to soften, the dominant length scale parameter (and hence shear band width) transitions from the alloy grain size to the spacing between micro-voids. Emphasis is placed on modeling stress collapse due to micro-void damage while simultaneously capturing the appropriate scale of inhomogeneous deformation. The goal is to assist in the microscale optimization of alloys which are susceptible to shear band failure.     

The growth of unstable thermoplastic shear with application to steady-wave shock compression in solids*

The catastrophic growth of unstable thermoplastic shear following the transition from homogeneous deformation to heterogeneous localized deformation through distributed shear banding is studied through approximate analytic and computational methods. The calculations provide expressions for shear band widths, spacing, catastrophic growth times and the rate of stress communication between shear bands. The optimum shear band width and spacing are found to be consistent with a minimum work principle. The model predicts that the product of the energy dissipated and the localization time in the shear localization process is invariant with respect to changes in the driving strain rate. Such behavior has been noted in the steady-wave shock compression of a number of solids. The calculations are applied to heterogeneous shear localization observed in the shock compression of aluminum.     

表面加工塑性层对金属柱壳剪切带自组织单旋起始的影响

在外爆加载金属柱壳高速坍塌过程中, 发生塑性变形失稳形成的剪切带具有高度的自组织特征, 甚至出现剪切带排列的单旋现象—剪切带在顺时针和逆时针两个方向呈现一个方向占优的现象. 柱壳在坍塌时, 最大剪切应力位于柱壳内表面, 剪切带的形核及扩展行为受内表面材料介观状态的影响显著. 本文通过选材和控制柱壳加工工艺, 获得了内表面具有不同厚度塑性层的20钢柱壳, 采用厚壁圆筒实验技术, 研究了表面加工塑性层对金属柱壳绝热剪切带自组织单旋现象起始的影响规律及其物理机制. 研究结果表明, 金属柱壳内表面加工塑性层显著改变了试样剪切带的起始条件, 沿顺时针或逆时针方向排列的剪切带形核数量在总剪切带数量中所占比例取决于表面加工塑性层的厚度和晶粒取向, 具有单一晶粒拉伸方向的较厚塑性层样品更容易形成单向螺旋剪切带结构. 在相同变形条件下, 随着塑性层厚度增加, 剪切带平均形核速率和扩展速率增大, 剪切带平均间距减小. 结果可为理解金属柱壳在高速塌陷过程中绝热剪切带占优取向现象提供有价值的参考.      

金属材料平面应变拉伸变形局部化有限元分析

本文对于涉及韧性金属大变形中颈缩与剪切带断裂一类高度非线性变形局部化问题进行了弹塑性有限元数值模拟。采用改进的Ｊ２形变理论微分形式公式与交叉三角形四边形单元有限元网格，详细研究了应变硬化指数及初始表面不均匀特性的平面应变拉伸颈缩和剪切带形成的综合影响，给出此类问题的断裂机制图。     

Macroscopic constitutive relations for elastomers reinforced with short aligned fibers: Instabilities and post-bifurcation response

This paper is concerned with the characterization of the macroscopic response and possible development of instabilities in a certain class of anisotropic composite materials consisting of random distributions of aligned rigid fibers of elliptical cross section in a soft elastomeric matrix, which are subjected to general plane strain loading conditions. For this purpose, use is made of an estimate for the stored-energy function that was derived by Lopez-Pamies and Ponte Castañeda (2006b) for this class of reinforced elastomers by means of the second-order linear comparison homogenization method. This homogenization estimate has been shown to lose strong ellipticity by the development of shear localization bands, when the composite is loaded in compression along the (in-plane) long axes of the fibers. The instability is produced by the sudden, collective rotation of a band of fibers to partially release the high stresses that develop in the elastomer matrix when the composite is compressed along the stiff, long-fiber direction. Consistent with the mode of the impending instability, a lower-energy, post-bifurcation solution is constructed where “striped domain” microstructures consisting of layers with alternating fiber orientations develop in the composite. The volume fractions of the layers and the fiber orientations within the layers adjust themselves to satisfy equilibrium and compatibility across the layers, while remaining compatible with the imposed overall deformation. Mathematically, this construction is shown to correspond to the rank-one convex envelope of the original estimate for the energy, and is further shown to be polyconvex and therefore quasiconvex. Thus, it corresponds to the “relaxation” of the stored-energy function of the composite, and can in turn be viewed as a stress-driven “phase transition,” where the symmetry of the fiber microstructures changes from nematic to smectic.     

The formation of multiple adiabatic shear bands

In a previous paper, Zhou et al. [2006. A numerical methodology for investigating adiabatic shear band formation. J. Mech. Phys. Solids, 54, 904-926] developed a numerical method for analyzing one-dimensional deformation of thermoviscoplastic materials. The method uses a second order algorithm for integration along characteristic lines, and computes the plastic flow after complete localization with high resolution and efficiency. We apply this numerical scheme to analyze localization in a thermoviscoplastic material where multiple shear bands are allowed to form at random locations in a large specimen. As a shear band develops, it unloads neighboring regions and interacts with other bands. Beginning with a random distribution of imperfections, which might be imagined as arising qualitatively from the microstructure, we obtain the average spacing of shear bands through calculations and compare our results with previously existing theoretical estimates. It is found that the spacing between nucleating shear bands follows the perturbation theory due to Wright and Ockendon [1996. A scaling law for the effect of inertia on the formation of adiabatic shear bands. Int. J. Plasticity 12, 927-934], whereas the spacing between mature shear bands is closer to that predicted by the momentum diffusion theory of Grady and Kipp [1987. The growth of unstable thermoplastic shear with application to steady-wave shock compression in solids. J. Mech. Phys. Solids 35, 95-119]. Scaling laws for the dependence of band spacing on material parameters differ in many respects from either theory.     

Adiabatic shear localization in the impact of edge-notched specimens

Impact experiments are performed on edgenotched specimens in the two-dimensional punch geometry. Materials tested include 18Ni(350) maraging steel; S7 tool steel; 4340, 300M, HP 9-4-20 and D-6ac ultra high-strength steels; and Ti6Al4V alloy. These materials have shown a high susceptibility to dynamic shear failure in previous studies. Impact velocity ranged from 25 m/s to 45 m/s, and shear bands were found to form at the notch tip and at the die corner on the back side of the specimen for all materials tested. Metallurgical analysis confirms the existence of adiabatic shear bands followed by a crack propagating through the fully developed shear band. High-speed photography was used to observe the initiation of adiabatic shear bands shortly after impact. Laser-etched lines on the specimen surfaces allowed the determination of the time of impact and the initiation time of shear failure. The elapsed time between the two was used to estimate the stress intensity factor at the time of shear band initiation. Comparisons of shear band initiation stress intensity factors at the notch tip and die corner are made. It is seen that the shear bands initiate at approximately the same stress intensity factor at both the notch tip and die corner. Finite element simulations support the use of a square root singularity for the stress in the plate near the corners of a deformable punch or die.     

再结晶组织对TA2纯钛绝热剪切行为的影响

使用二辊轧机对TA2工业纯钛进行多道次大应变冷轧处理,制备了冷轧总变形量为70%的TA2纯钛板。通过对冷轧TA2纯钛板进行500℃加热、不同保温时间的退火处理,获得了具有不同再结晶组织的钛板。基于帽形试样和限位环变形控制技术,在分离式霍普金森压杆装置上对不同再结晶组织的试样进行动态冲击冻结实验,结合光学显微镜和扫描电子显微镜表征试样冲击前后微观组织的变化,研究了再结晶组织对TA2纯钛绝热剪切行为的影响。结果表明,随着退火保温时间的延长,试样再结晶晶粒占比逐渐增大,晶粒分布由分散向局部聚集转变;在相同应变和应变率下,在所有试样中都观察到了绝热剪切带,再结晶晶粒占比高的试样更易诱发绝热剪切带中裂纹形核扩展。对比变形前后试样再结晶组织和几何必需位错变化,结合剪切区整体温升分析发现,再结晶晶粒作为材料软化点能够诱发剪切带的形成,而剪切带发展后期产生的绝热温升会促进剪切带内材料发生二次再结晶,提高剪切带内材料的韧性,延缓剪切裂纹的形成。     

A new constitutive model for shear banding instability in metallic glass

Inelastic deformation of metallic glass is through shear banding, characterized by significantly localized deformation and emerged expeditiously under certain stress state. This study establishes a new constitutive model addressing the physical origin of the shear banding. In the modeling, the atomic structural change and the free volume generation are embodied by the plastic shear strain and the associated dilatation. The rugged free energy landscape is adopted to naturally reflect the rate-independent flow stress and flow serrations. Based on this, the conditions for the onset of shear banding instability are established, which enables the explicit calculation of the shear band inclination angle and its extension speed. The study concludes that shear band angle is significantly influenced by the diltancy factor and pressure sensitivity, that a shear band does not increase its thickness once emanated from a deformation unit, that the spreading speed of a shear band is intersonic, and that more shear bands, which lead to higher ductility, can be induced by high strain rates or by the introduction of a second material phase. The analysis also demonstrates that the ductility of metallic glass depends on the sample geometry and/or the stress state.     

表面粗糙度对TC4钛合金柱壳剪切带形成的影响

剪切带是材料在高应变率加载条件下特有的变形和损伤形式之一,关于影响金属材料中剪切带形成的敏感性因素及其机理的研究,一直是科学研究和工程设计中关注的重点问题. 在柱壳高速坍塌过程中,剪切带优先在内表面形核, 其形核及扩展行为受内表面介观状态的影响显著.本文采用爆轰加载厚壁圆筒坍塌实验技术,结合材料表面处理技术、微结构表征技术和剪切带理论模型分析,研究了内表面粗糙度变化对TC4钛合金柱壳剪切带形成影响的细观动力学规律.结果表明, 在爆炸加载形成的高应变率条件下,表面粗糙度对TC4钛合金柱壳中剪切带形成具有明显影响. 在相同的变形条件下,随着试样内表面粗糙度的增大, 剪切带数量、长度和形核速率均增大;表面粗糙度越大, 部分剪切带扩展速率越快, 剪切带长度差异越大,剪切带的屏蔽效应增强. 分析表明,实验获得的剪切带间距与W-O模型和M模型预测结果基本吻合,具体数值受试样内表面粗糙度影响, 随着表面粗糙度的增大,实验结果逐渐小于预测数值.      

Post-critical plastic deformation in incrementally nonlinear materials

The formation of multiple macroscopic shear bands is investigated as a mechanism of advanced plastic flow of polycrystalline metals. The overall deformation pattern and material characteristics are determined beyond the critical instant of ellipticity loss, without the need of introducing an internal length scale. This novel approach to the modelling of post-critical plastic deformation is based on the concept of a representative nonuniform solution in a homogeneous material. The indeterminacy of a post-critical representative solution is removed by eliminating unstable solution paths with the help of the energy criterion of path instability. It is shown that the use of micromechanically based, incrementally nonlinear corner theories of time-independent plasticity leads then to gradual concentration of post-critical plastic deformation. The volume fraction occupied by shear bands is found to have initially a well-defined, finite value insensitive to the mesh size in finite element calculations. Further deformation depends qualitatively on details of the constitutive law. In certain cases, the volume fraction of active bands decreases rapidly to zero, leading to material instability of dynamic type. However, for physically hardening materials with the yield-vertex effect, the localization volume typically remains finite over a considerable deformation range. At later stages of the plane strain simulation, differently aligned secondary bands are formed in a series of bifurcations.     

Numerical modeling of the thermomechanical behavior of steels with allowance for the propagation of Luders bands

A thermomechanical model based on physical representations of the motion of dislocation continuum and a model for the initiation and propagation of plastic shear are proposed to describe slow flows of the type of Luders bands. Two-dimensional calculations of Luders band propagation are performed for HSLA-65 steel samples under compression at various strain rates and temperatures. The calculation results are in good agreement with experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 5, pp. 146–155, September–October, 2007.     

Simulation of crack propagation in fiber-reinforced bulk metallic glasses

Based on a phase-field model for deformation in bulk metallic glasses (BMGs), shear band formation and crack propagation in the fiber-reinforced BMG are investigated. Ideal unbroken fibers embedded in the BMG matrix are found to significantly influence the shear banding and crack propagation in the matrix. The crack propagation affected by fibers’ length and orientation is quantitatively characterized and is described by micromechanics models for composite materials. Furthermore, fractures in some practical fiber-reinforced BMG composites such as tungsten-reinforced Zr-based BMG are simulated. The relation between the enhanced fracture toughness and the mechanical properties of fiber reinforcements is determined. Different fracture modes of BMG-matrix composites are identified from the systematic simulation studies, which are found to be consistent with experiments. The simulation results suggest that the phase-field modeling approach could be a useful tool to assist the fabrication and design of BMG composites with high fracture toughness and ductility.     

A NOVEL APPROACH TO TESTING THE DYNAMIC SHEAR RESPONSE OF Ti-6Al-4V

Modifications were made on the traditional split Hopkinson pressure bar (SHPB) system to conduct dynamic shear tests. The shear response of Ti-6Al-4V was acquired at a shear strain rate of 10 4 s 1 by using this modified apparatus. The geometry as well as the clamping mode of the double-notch specimen was optimized by commercial FEM software ABAQUS, and the feasibility of the experiment set-up was validated. A shear stress calibration coefficient of τ = 1.03 and a shear strain calibration coefficient of Γ = 0.50 were obtained.We have employed high-speed photography to record the deformation process, especially the initiation and propagation of adiabatic shear band (ASB), during the dynamic shear test. The frames show that the time duration from ASB initiation to its completion is less than 2 μs, from which we can estimate that the propagation speed of ASB within Ti-6Al-4V is more than 1250 m/s under such loading conditions. The temperature rise within ASB is also estimated to be T 2 ≈ 1460℃ based on energy balance. Such high temperature has led to softening of the material within the ASBs, and has intensified the shear localization and finally resulted in fracture of the material.     

A failure criterion based on material instability

Scaling laws for adiabatic shear bands are used to parameterize a model that is suitable for introducing shear damage within engineering calculations. One-dimensional solutions to the governing equations for a single shear band provide laws that connect the driving deformation, the imperfections, and the physical characteristics of the material to the process of stress collapse [International Journal of Plasticity 8 (1992) 583, Mechanics of Materials 17 (1994) 215]. The current model uses homogeneous material response and the scaling laws to anticipate the correct timing beyond the maximum stress at which stress collapse should occur. The model is implemented into a finite element code for wave propagation and used in the analysis of boundary value problems that are dominated by shear failure. Finally, implications of the model for simulations of material failure are discussed.     

高速冲击表面处理对金属材料力学性能和组织结构的影响

高速冲击表面处理过程中的应变率对金属材料的宏观力学性能和微观组织结构都具有重要影响。根据当前应变率效应的研究成果,从宏观与微观相结合的角度出发,综述了高速冲击表面处理过程中应变率对金属材料强度和塑性的影响规律,并重点阐述了不同应变率下金属材料内部微观组织结构的演变规律,主要包括晶粒结构、绝热剪切带、相变、位错组态和析出相以及变形孪晶等。此外,还分析了组织结构随应变率的演化和微观变形机制的转变对材料力学性能的强化和弱化机理。最后,对高速冲击表面处理梯度组织的变形特点进行了总结。提出了不同组织结构对材料性能影响的综合效应模型,以期为应变率效应的深入研究奠定基础。