Effect of porosity on the properties of strain localization in porous media under undrained conditions

Within the general framework of mixture theory and by introducing the fictitious “fluid phase” as a mixture of a liquid and a gas, the conditions for localization of deformation into a shear band in the incremental response of partially saturated and fully saturated elastic–plastic porous media under undrained conditions are derived. The effect of porosity is included in the derivation. The explicit analytical expressions of the direction of shear band initiation and the corresponding hardening modulus of the porous media for the plane strain case are deduced, and a parametric analysis is made of the influence of the porosity on the properties of strain localization based on Mohr–Coulomb yield criterion. It is found that the dependence of the shear banding properties of partially saturated porous media on the porosity is related to the stress states and Poisson's ratio. However, the properties of the strain localization for the fully saturated porous media are almost independent of Poisson's ratio. Finally, on the basis of Mohr–Coulomb yield criterion, some solutions of the shear banding orientation for water-saturated granular materials are obtained, which are proved to be in good agreement with the experimental results reported by other researchers.     

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

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

Prediction of the initial thickness of shear band localization based on a reduced strain gradient theory

Plastic flow localization in ductile materials subjected to pure shear loading and uniaxial tension is investigated respectively in this paper using a reduced strain gradient theory, which consists of the couple-stress (CS) strain gradient theory proposed by Fleck and Hutchinson (1993) and the strain gradient hardening (softening) law (C–W) proposed by Chen and Wang (2000). Unlike the classical plasticity framework, the initial thickness of the shear band and the strain rate distribution in both cases are predicted analytically using a bifurcation analysis. It shows that the strain rate is obviously non-uniform inside the shear band and reaches a maximum at the center of the shear band. The initial thickness of the shear band depends on not only the material intrinsic length l_cs but also the material constants, such as the yield strength, ultimate tension strength, the linear hardening and softening shear moduli. Specially, in the uniaxial tension case, the most possible tilt angle of shear band localization is consistent qualitatively with the existing experimental observations. The results in this paper should be useful for engineers to predict the details of material failures due to plastic flow localization.     

A finite element approach to study cavitation instabilities in non-linear elastic solids under general loading conditions

This paper proposes an effective numerical method to study cavitation instabilities in non-linear elastic solids. The basic idea is to examine—by means of a 3D finite element model—the mechanical response under affine boundary conditions of a block of non-linear elastic material that contains a single infinitesimal defect at its center. The occurrence of cavitation is identified as the event when the initially small defect suddenly grows to a much larger size in response to sufficiently large applied loads. While the method is valid more generally, the emphasis here is on solids that are isotropic and defects that are vacuous and initially spherical in shape. As a first application, the proposed approach is utilized to compute the entire onset-of-cavitation surfaces (namely, the set of all critical Cauchy stress states at which cavitation ensues) for a variety of incompressible materials with different convexity properties and growth conditions. For strictly polyconvex materials, it is found that cavitation occurs only for stress states where the three principal Cauchy stresses are tensile and that the required hydrostatic stress component at cavitation increases with increasing shear components. For a class of materials that are not polyconvex, on the other hand and rather counterintuitively, the hydrostatic stress component at cavitation is found to decrease for a range of increasing shear components. The theoretical and practical implications of these results are discussed.     

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.     

分子动力学模拟金属玻璃板拉伸时剪切带的形成和演化过程

金属玻璃在低温高应力条件下容易生成剪切带而导致结构的破坏,大大限制了它的推广应用。本文采用分子动力学模拟研究了三种Cu64Zr36（不带缺口、一侧带缺口、两侧带缺口）金属玻璃板试样在拉伸过程中剪切带的形成和演化过程及其力学性能。结果表明：不带缺口金属玻璃板试样在低温高应力的拉伸过程中会自发出现局部剪切转变区,发生剪切局部化,继续拉伸会在与加载轴大约成45°方向上形成剪切带。剪切带的形成与剪切转变区的分布和局部化有关,带缺口比不带缺口的试样会更早出现应变局部化,即在较低的拉伸应变下便形成剪切带,其拉伸强度也相应较低。相同条件下,一侧带缺口与两侧带缺口的试样在拉伸强度上几乎相同,但两侧带缺口试样的应变局部化程度稍低,主要是两侧缺口处均出现了剪切转变区,导致其分布和局部化不够集中,这也是形成主剪切带和次剪切带的主要原因。以上结果为进一步从原子尺度理解金属玻璃剪切带的形成和演化特征提供了重要的信息。     

分子动力学模拟金属玻璃板拉伸时剪切带的形成和演化过程

金属玻璃在低温高应力条件下容易生成剪切带而导致结构的破坏,大大限制了它的推广应用。本文采用分子动力学模拟研究了三种Cu64Zr36（不带缺口、一侧带缺口、两侧带缺口）金属玻璃板试样在拉伸过程中剪切带的形成和演化过程及其力学性能。结果表明：不带缺口金属玻璃板试样在低温高应力的拉伸过程中会自发出现局部剪切转变区,发生剪切局部化,继续拉伸会在与加载轴大约成45°方向上形成剪切带。剪切带的形成与剪切转变区的分布和局部化有关,带缺口比不带缺口的试样会更早出现应变局部化,即在较低的拉伸应变下便形成剪切带,其拉伸强度也相应较低。相同条件下,一侧带缺口与两侧带缺口的试样在拉伸强度上几乎相同,但两侧带缺口试样的应变局部化程度稍低,主要是两侧缺口处均出现了剪切转变区,导致其分布和局部化不够集中,这也是形成主剪切带和次剪切带的主要原因。以上结果为进一步从原子尺度理解金属玻璃剪切带的形成和演化特征提供了重要的信息。     

基于晶体塑性理论研究铝材料高压高应变率下的强度特性

为了了解金属材料在极端加载下复杂动态响应过程中的多种机制和效应,重点针对Al材料在高压、高应变率加载下的塑性变形机制,在经典晶体塑性模型的基础上,对其中的非线性弹性、位错动力学和硬化形式进行改进,建立适用于高压、高应变率加载下的热弹-黏塑性晶体塑性模型。该模型可以较好地描述单晶铝和多晶铝材料屈服强度随压力的变化过程,相比宏观模型,用该模型还获得了多晶Al材料在冲击加载下的织构演化规律,揭示了织构择优取向行为和压力的关系。     

The void-size effect on plastic flow localization in the Gurson model

Recent studies have shown that the size of microvoids has a significant effect on the void growth rate. The purpose of this paper is to explore whether the void size effect can influence the plastic flow localization in ductile materials. We have used the extended Gurson‘s dilatational plasticity theory, which accounts for the void size effect, to study the plastic flow localization in porous solids with long cylindrical voids. The localization model of Rice is adopted, in which the material inside the band may display a different response from that outside the band at the incipient plastic flow localization. The present study shows that it has little effect on the shear band angle.     

帽型试样动态绝热剪切破坏演化分析

利用分离式霍普金森压杆加载Ta2钛合金扁平闭合帽形受迫剪切试样,结合数字图像相关法和“冻结”试样的微观金相观察,研究剪切区剪切应变的演化、绝热剪切带形成条件等。结果显示：受迫剪切试样在动态加载过程,剪切区剪切应变不断集中,形成绝热剪切带,裂纹沿绝热剪切带发展;随加载率提高,绝热剪切起始临界应变减小;进一步利用数字图像相关法DIC场应变分析及金相微观观测对比,利用卸载回复特性对绝热剪切带起始临界条件进行了讨论,计算的绝热剪切带起始时温升仅为86℃。材料软化可能不是绝热剪切带起始的控制条件,相反是由于绝热剪切带形成造成的应变高度集中发展导致温度急剧升高。     

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

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

Simulation of Cavitation in Water Saturated Porous Media Considering Effects of Dissolved Air

An extension of a mathematical model for non-isothermal multiphase materials to consider the dissolution of air in liquid water and air mass sources during its desorption at lower water pressure is presented. The solid skeleton is assumed elasto-plastic; heat, water and air flows and water phase changes are taken into account. Physics of air dissolution and desaturation due to the air released from liquid water during cavitation in porous media are discussed. A numerical example where cavitation develops during shear band development in undrained water-saturated dense sands is solved with the developed model as discretized in space and time with the Finite Element Method.     

饱和多孔介质分析解的唯一性与应变局部化分岔

基于不连续性分岔基本理论导出了静态非渗流状态下弹塑性饱和多孔介质应变局部化发生时的临界硬/软化模量,利用二阶功正定性原理研究了两相问题分析解的唯一性问题,并给出了基于主轴空间下解的显式表达式。研究工作表明,在静态非渗流状态下,弹塑性饱和多孔介质分析的唯一性与应变局部化发生的临界条件除了在量值上与单相介质有着明显的不同外两者之间还有许多一致的特性,这些一致的特性对问题的分析是十分重要的。     

侧斜与负载对螺旋桨无空化和空化水动力性能的影响

为了定量分析空化初始航速的影响因素,首先分析侧斜和负载对螺旋桨无空化和有空化时性能的影响。以NSRDC4381无侧斜桨和4383100％侧斜桨为对象,采用改进Sauer空化模型和修正SST湍流模型,对空化崩溃性能、空化初生性能和无空化时正车、倒车以及紧急倒车敞水性能进行了计算与比较。结果表明,预报两个桨的敞水性能曲线和多空化数下的空化崩溃性能曲线均与实验值吻合较好。在中度负载区间（J=0．5-0．9）内,侧斜对正车和紧急倒车时敞水性能以及空化时推力和力矩崩溃性能均无明显影响,但会使倒车敞水性能显著下降。在重载和轻载条件下,侧斜均能明显改善空化初生性能。侧斜一定时,负载会直接影响尾流湍流速度脉动量和涡核集中区分布,影响轴面速度流管收缩程度,进而影响无空化和有空化条件下的推进性能。     

Shear band analysis and shear moduli calibration

Strain localization is a well known phenomenon, generally associated with plastic deformation and rupture in solids, especially in geomaterials. In this process, deformation is observed to concentrate in narrow zones called shear bands. This phenomenon has been studied extensively in the last 20 years by different researchers, experimentally, theoretically and numerically. A criterion for the onset of localization can be predicted solely on the basis of the constitutive law of the material, using the so-called shear band analysis. This criterion gives the critical orientation, and the critical stress state and strain for a given loading history. An important point, already stressed by Vardoulakis in 1980, is that in particular, out-of-axes shear moduli play a central role in the criterion. These are the moduli involved in the response to a deviatoric stress increment with principal axes oriented at 45° from total stress principal axes. Out-of-axes shear moduli are difficult parameters to calibrate; common tests, with fixed principal stress and strain directions, do not provide any information on these moduli, as long as they remain homogeneous. Still, real civil engineering and environmental problems are definitely not simple axisymmetric triaxial tests; practical modeling involves complex stress paths, and need complex parameters to be calibrated. Only special tests, like compression–torsion on hollow cylinder tests, or even more complex tests can be used for shear moduli calibration. However, shear band initiation in homogeneous, fixed-axes tests does activate out-of-axes shear. Hence, it is natural that shear band analysis makes shear moduli enter into the analysis.Then, a typical inverse analysis approach can be used here: experimental observation of strain localization in triaxial tests can be used together with a proper shear band analysis for the model considered, in order to determine out-of-axes shear moduli.This approach has been used for a stiff marl in the framework of a calibration study on a set of triaxial tests. The steps of the method are presented, and the bifurcation surface in the stress space is exhibited.     

Plastic flow in oriented glassy polymers

A manufactured product often possesses residual texture which was either incidentally or deliberately acquired during its processing history. This is particularly true for the case of polymers, where the ability to easily preferentially preorient the material in specific directions is exploited in order to obtain a higher strength product. Specific examples include synthetic fibers, and biaxially-oriented films and containers. The response of the preoriented/textured product to normal service life loading conditions will differ considerably from that of a product composed of isotropic material. This paper addresses the issue of the effects of texture on the deformation behavior of glassy polymers. Here, the physically-based constitutive model of Boyce, Parks, and Argon describing the rate, temperature, and pressure dependent inelastic deformation of initially isotropic glassy polymers is used to model the effects of preorientation (i.e., initial texture), via the use of appropriate initial conditions on internal state variables. The model is then utilized in an analysis of the effects of texture on the yield of glassy polymers and the shear localization which normally follows yielding in oriented polymers. These results are compared with trends found in experiments as reported in the literature. The effectiveness of the proposed model for the present application is also compared with earlier models of yielding of anisotropic materials such as Hill's criterion.     

The growth of shear bands in composite microstructures

Shear band formation in materials with inhomogeneous and composite microstructures is influenced by factors that usually do not come into play in monolithic materials. Experiments and calculations have shown that inhomogeneities in material properties enhance the localization of deformation. This investigation concerns the propagation of shear bands in a two-phase tungsten composite under the conditions of nominally pure shear deformation. Finite element calculations are carried out to delineate the effects of different grain–matrix morphologies. In the numerical models, the initiation of shear bands is triggered by a notch, simulating the effect of defects such as microcracks and microvoids in materials. Calculations demonstrate that phase morphology, particle size and the relative location of initiation site have significant influences on the development of localized deformation. The work and energy evolutions are tracked for each constituent phase in the microstructures. In addition, the exchange of thermal energy through heat flow between the phases is analyzed. The results show that a strong correlation exists between the course of shear band propagation and the thermomechanical coupling between microscopic phases.     

Multi-slip gradient formulation for modeling microstructure effects on shear bands in granular materials

This paper presents a higher order gradient multi-slip formulation to model the effect of inhomogeneous deformation in granular materials. The effects of heterogeneity and porosity anisotropy within the multi-slip formulation are taken into consideration through the modification of the mobilized friction. The mobilized friction is assumed to be a direct function of either the gradient of the porosity distribution or the fabric tensor. The formulation with two active slip planes was implemented into a finite element code and used to simulate biaxial shear tests on dry sand. The analysis quantifies most of the shear band characteristics observed by past experimentation. It is shown that the localization and shear band characteristics in granular materials are very much dependent on the initial fabric and slip system arrangement.     

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.     

Modeling shear localization along granular soil–structure interfaces using elasto-plastic Cosserat continuum

The current study presents finite element simulations of shear localization along the interface between cohesionless granular soil and bounding structure under large shearing movement. Micro-polar (Cosserat) continuum approach is applied in the framework of elasto-plasticity in order to overcome the numerical problems of localization modeling seen in the conventional continuum mechanics. The effects of different micro-polar kinematic boundary conditions, along the interface, on the evolution and location of shear band are shown by the numerical results. Furthermore, shear band thickness is also investigated for its dependence on the initial void ratio, vertical pressure and mean grain size. Here, the distribution and evolution of static and kinematic quantities are the main focuses regarding infinite layer of micro-polar material during plane shearing, especially with advanced large movement of bounding structure. The influence of such movement has not been investigated yet in the literature. Based on the results obtained from this study, shear localization appears parallel to the direction of shearing. It occurs either in the middle of granular layer or near boundaries, regarding the assumed micro-polar kinematic boundary conditions at the bottom and top surfaces of granular soil layer. Narrower shear band is observed in lower rotation resistance of soil particles along the interface. It is emphasized that the displacement magnitude of bounding structure has significant effect on the distribution and evolution of state variables and polar quantities in the granular soil layer. However, continuous displacement has no meaningful effect on the thickness of shear band. Here, smooth distributions of void ratio and shear stress components are obtained within the shear band, what the other previous numerical investigations did not receive. Despite indirect linking of Lade’s model to the critical state soil mechanics, state variables tend towards asymptotical stationary condition in large shear deformation.