Shear banding analysis of geomaterials by strain gradient enhanced damage model

Shear band localization is investigated by a strain-gradient-enhanced damage model for quasi-brittle geomaterials. This model introduces the strain gradients and their higher-order conjugate stresses into the framework of continuum damage mechanics. The influence of the strain gradients on the constitutive behaviour is taken into account through a generalized damage evolutionary law. A weak-form variational principle is employed to address the additional boundary conditions introduced by the incorporation of the strain gradients and the conjugate higher-order stresses. Damage localization under simple shear condition is analytically investigated by using the theory of discontinuous bifurcation and the concept of the second-order characteristic surface. Analytical solutions for the distributions of strain rates and strain gradient rates, as well as the band width of localised damage are found. Numerical analysis demonstrates the shear band width is proportionally related to the internal length scale through a coefficient function of Poisson’s ratio and a parameter representing the shape of uniaxial stress–strain curve. It is also shown that the obtained distributions of strains and strain gradients are well in accordance with the underlying assumptions for the second-order discontinuous shear band boundary and the weak discontinuous bifurcation theory.     

Micromechanics based second gradient continuum theory for shear band modeling in cohesive granular materials following damage elasticity

Gradient theories, as a regularized continuum mechanics approach, have found wide applications for modeling strain localization failure process. This paper presents a second gradient stress–strain damage elasticity theory based upon the method of virtual power. The theory considers the strain gradient and its conjugated double stresses. Instead of introducing an intrinsic material length scale into the constitutive law in an ad hoc fashion, a microstructural granular mechanics approach is applied to derive the higher-order constitutive coefficients such that the internal length scale parameter reflects the natural granularity of the underlying material microstructure. The derivations of the required damage constitutive relationships, the strong form governing equations as well as its weak form for the second gradient model are described. The recently popularized Element-Free Galerkin (EFG) method is then employed to discretize the weak form equilibrium equation for accommodating the resultant higher-order continuity requirements and further handling the mesh sensitivity problem. Numerical examples for shear band simulations show that the proposed second gradient continuum model can produce stable, accurate as well as mesh-size independent solutions without a priori assumption of the shear band path.     

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

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

Cosserat连续体模型中本构参数对应变局部化模拟结果影响的数值分析

基于所发展的压力相关弹塑性Cosserat连续体模型及相应的数值方法,以一维剪切层及二维平板压缩问题为例,数值分析了Cosserat连续体模型中的本构参数Cosserat剪模、软化模量及内部长度参数对应变局部化数值模拟结果的影响.结果表明在一定取值范围内,Cosserat剪模对数值模拟结果几乎没有影响,并给出了具体数值计算时的取值范围;软化模量绝对值越大,后破坏段的荷载-位移曲线越陡,计算得到的剪切带宽度越窄;内部长度参数越大,后破坏段的荷载-位移曲线越平缓,计算得到的剪切带越宽.     

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.     

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.     

Small scale effects on buckling and postbuckling behaviors of axially loaded FGM nanoshells based on nonlocal strain gradient elasticity theory

By means of a comprehensive theory of elasticity, namely, a nonlocal strain gradient continuum theory, size-dependent nonlinear axial instability characteristics of cylindrical nanoshells made of functionally graded material (FGM) are examined. To take small scale effects into consideration in a more accurate way, a nonlocal stress field parameter and an internal length scale parameter are incorporated simultaneously into an exponential shear deformation shell theory. The variation of material properties associated with FGM nanoshells is supposed along the shell thickness, and it is modeled based on the Mori-Tanaka homogenization scheme. With a boundary layer theory of shell buckling and a perturbation-based solving process, the nonlocal strain gradient load-deflection and load-shortening stability paths are derived explicitly. It is observed that the strain gradient size effect causes to the increases of both the critical axial buckling load and the width of snap-through phenomenon related to the postbuckling regime, while the nonlocal size dependency leads to the decreases of them. Moreover, the influence of the nonlocal type of small scale effect on the axial instability characteristics of FGM nanoshells is more than that of the strain gradient one.     

不同加载状态下TA2钛合金绝热剪切破坏响应特性

一般认为绝热剪切现象在宏观上表现为材料动态本构失稳,即热软化大于应变硬化.本文采用帽型受迫剪切试样研究TA2钛合金的动态力学特性和本构失稳过程.首先对剪切区加载应力状态进行理论和数值分析,通过合理设计帽型试样,剪切区变形可近似按剪切状态处理;结合二维数字图像相关法(two-dimensional digital image correlation,DIC-2D)直接测试试样剪切区应变演化,给出帽型受迫剪切实验的等效应力-应变响应曲线.进一步,利用Hopkinson压杆对TA2钛合金开展动态压缩及帽型剪切对比试验研究,比较压缩、剪切试验得到的等效应力-应变曲线,采用"冻结"试样方法分析试样中绝热剪切局域化演化过程,探讨不同加载状态下TA2钛合金的绝热剪切破坏现象及其动态力学响应特性.实验结果表明,在塑性变形初始阶段,动态压缩及剪切加载下的等效应力-应变曲线符合较好,但随塑性损伤发展及绝热剪切带形成,两者出现分离,表明损伤及绝热剪切演化过程与应力状态相关.剪切试样实验得到的本构"软化"特性能够反映绝热剪切带起始、破坏演化过程的力学响应特性,而在动态压缩实验中,即使试样中已出现双锥形的绝热剪切带及局部裂纹分布,其表观等效应力-应变曲线并不出现软化特征,动态压缩实验无法得到关于绝热剪切起始、发展以及破坏的本构软化响应特性.     

Multiresolution modeling of ductile reinforced brittle composites

Tungsten carbide-cobalt (WC-Co) is an important ductile reinforced brittle composite used in a range of important applications. The relationship between microstructure and mechanical properties of WC-Co is truly multiscale; micromechanical processes interact at different scales, resulting in permanent plastic deformation, damage accumulation and final failure of the composite. The goal of the current paper is to develop a continuum-based model, which captures the progressively finer scales of strain localization observed in WC-Co composites during plastic deformation and failure. This is achieved via a set of multiresolution governing equations; a microstress is introduced at each scale of strain localization, which represents the resistance to inhomogeneous strain localization at that scale. The extra constitutive models associated with these microstresses can be elucidated from the average response of separate computational cell models of a representative microstructure. The final multiresolution continuum model is capable of capturing the important length scales of deformation during the plastic stage of deformation without resorting to modeling microstructural scale features directly. The result is a more realistic continuum model; in particular the fracture toughness prediction is more physical when these length scales are incorporated compared to a conventional continuum approach.     

Evolving internal length scales in plastic strain localization for granular materials

A numerical model in the Cosserat continuum for strain localization phenomena in granular materials is developed and proposed in this paper. The model assumes a constant internal length scale that is used to describe the shear band thickness. However, it is observed that the internal length scales need to change to accommodate the possible change in the contact surface between the particles, damage of the particles or/and any change in the local void ratio within the domain, which will change the shear band thickness. The mathematical formulations used in the present numerical model were equipped with evolution equations for the length scales through the Micropolar theory, those formulations are proposed and discussed in this paper. The evolution equations of the internal length scales describe any possible change in the contact surface between the particles, damage of the particles if exists and/or any change in the local void ratio within the domain. Hence, the strain localization described by the enhanced model with evolving internal length scales is more accurate and closer to the real solution. The solution for the shear bands thickness shows more accurate correlation with the experimental results and less dependency on the mesh size when such evolution equations are used. Moreover, the shear band thickness and inclination evolve during the deformation process.     

Failure criteria with unilateral conditions for simulation of plate perforation

Finite element analysis of plate perforation is performed using a thermoelastic–thermoviscoplastic constitutive model and two uncoupled damage criteria with unilateral conditions. The first criterion is based on continuum damage mechanics, while the second is the Cockcroft–Latham criterion. Fracture is assumed to occur at a critical damage, and crack propagation is modelled with element erosion. The identification of the material parameters is based on inverse modelling of a uniaxial tension test. The constitutive model and damage criteria are applied in explicit finite element analysis of the perforation of a 12 mm thick steel plate struck by a blunt projectile in the sub-ordinance velocity regime. The results show the importance of the unilateral conditions in obtaining accurate prediction of the adiabatic shear banding process. The two fracture criteria have comparable dependence on the stress triaxiality and the Lode parameter when the microcrack-closure parameter is set to zero in the criterion based on continuum damage mechanics, and give similar failure mode and ballistic limit in the plate perforation problem.     

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

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

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

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

饱和土中剪切带宽度的研究

主要对饱和土在简单剪切和耦合载荷下的剪切带宽进行了理论分析. 在模型中考虑了剪胀、孔压耗散、惯性及孔隙水和土颗粒之间阻力等因素的耦合效应,将有效应力考虑为应变、应变率和孔隙压力的函数. 通过对模型的简化分析得到简单剪切和耦合应力状态下的剪切带宽度的近似计算公式. 分析表明,当土体失稳后,会出现有限厚度的剪切带, 且其形成是一个后期过程,可以忽略惯性效应. 剪切带的厚度依赖于孔隙压力特性和土体骨架的剪胀特性,但是不依赖应力状态. 最后将分析与实验结果进行了对比,在一定程度上,两者是符合的.     

A rate-dependent damage model for brittle materials based on the dominant crack

A rate-dependent, continuum damage model is developed for brittle materials under dynamic loading. This model improves on the approach (ISOSCM) of [Addessio, F.L., Johnson, J.N., 1990. A constitutive model for the dynamic response of brittle materials. Journal of Applied Physics 67, 3275–3286] in several respects. (1) A new damage surface is found by applying the generalized Griffith instability criterion to the dominant crack (having the most unstable orientation), rather than by averaging the instability condition over all crack orientations as done previously. The new surface removes a discontinuity in the damage surface in ISOSCM when the pressure changes sign. (2) The strain due to crack opening is more consistent with crack mechanics, with only the tensile principal stresses contributing to the crack opening strain. This is achieved by incorporating a projection operator in the equation for the crack opening strain. One consequence of incorporating the projection operator is a prediction of shear dilatancy, which is not accounted for in ISOSCM. (3) The evolution of damage, which is based on the energy-release rate for the dominant crack, has a physical basis, whereas in the previous approach the damage growth rate was assumed to be an exponential function of the distance from the stress state to the damage surface without specific physical justification.An implicit algorithm has been developed so that a larger time step can be used than with the explicit algorithm used in ISOSCM. The numerical results of a silicon carbide (SiC) ceramic under several loading paths (hydrostatic tension/compression, uniaxial strain, uniaxial stress, and shear) and strain rates are presented to illustrate the main features of the model.     

土的工程力学性质的颗粒流模拟

基于颗粒流理论,引入不同的颗粒接触连接本构模型,分别建立了砂土和粘性土的颗粒流模型．通过颗粒流数值模型试验,对砂土和粘性土的室内平面应变试验及其剪切带形成和发展进行了数值模拟,分别对比了不同围压下颗粒流试样与室内试验的应力应变关系曲线,基本再现了砂土和粘性土试样应力．应变关系．通过砂土和粘性土PFC试样剪切带模拟表明,当围压较小时试样内部颗粒位移量小而且分布范围较广,当围压增大时,试样内部颗粒位移量也增大,而且发生较大位移颗粒的分布范围趋于集中,同时随着围压的增大试样内部形成明显的剪切带．无论砂土还是粘性土的PFC试样,随着围压的增加剪切带的形状趋于集中,而且剪切带宽度在减小．在围压很小时,试样内形成大的破坏区域,在围压较大时出现明显的线破坏区．这些规律基本与室内试验结果相似。     

Micro-macro scale instability in 2D regular granular assemblies

Instability and stress–strain behavior were investigated for 2D regular assemblies of cylindrical particles. Biaxial shear experiments were performed on three sets of assemblies with regular, albeit increasingly defective structures. These experiments revealed unique instability behavior of these assemblies. Continuum models for the assemblies were then constructed using the granular micromechanics approach. In this approach, the constitutive equations governing the behavior of inter-particle contacts are written in local or microscopic level. The behavior of the RVE is then retrieved by using either kinematic constraint or least squares (static constraint) along with the principle of virtual work to equate the work done by microscopic force–displacement conjugates to that of the macroscopic stress and strain tensor conjugates. The ability of the two continuum approaches to describe the measured stress–strain behavior was evaluated. The continuum models and the local constitutive laws were used to perform instability analyses. The onset of instability and orientation of shear band was found to be well predicted by the instability analyses with the continuum models. Further, macro-scale instability was found to correlate with the instability of inter-particle contacts, although with some variations for the two modeling approaches.     

Shear bands in dense metallic granular materials

In this paper, the phenomenon of strain localization, i.e., shear banding, in densely distributed metallic assemblies has been studied. A discrete element methodology for analyzing metallic granular materials has been put forward. In this numerical model, elastoplastic contact, as well as friction, rolling resistance and cohesion between spheres, are explicitly taken into account. The calculations reveal that the shear banding mechanism in dense assemblies can be thought as an instability triggered by initial imperfections. Within a band, the motion, deformation and rearrangement of spheres soften the resistance of the aggregate, as these mechanisms create additional geometric imperfections. Additionally, the simulations showed that the shear-band width does not change conclusively with the friction, rolling resistance and plasticity parameters. However, cohesive strength, even in small amounts, drastically increased the shear-band width. Finally, the shear-band thickness and inclination angles are strongly dependent on the degree of initial imperfection. Whereas for a perfect assembly the shear band inclinations were consistently around 60°, more heterogeneous assemblies lead to shear band angles closer to the continuum mechanics solution, which is 45°. This was found to be in agreement with recent experimental observations.     

The initiation and growth of macroscopic martensite band in nano-grained NiTi microtube under tension

The superelastic behavior of polycrystalline nano-grained NiTi shape memory alloy microtube under uniaxial tension is studied in this paper. The nominal stress–strain curve of the microtube during superelastic deformation is recorded. Both direct surface observation and observation by using a special surface coating show that the deformation of the tube is via the nucleation and propagation of macroscopic stress-induced martensite band. It is also found that the martensite nucleates in the form of a spiral lens-shaped narrow band that inclines at about 33^o to the plane of cross section of tube when the stress reaches the peak of stress–strain curve. The spiral band grew via gradual increase in both width and length of the band and finally merged into a single cylindrical band. The subsequent deformation of the tube is realized by the growth of this cylindrical martensite band. Several other deformation features of the tube are also observed and the results are discussed and compared with the theoretical analysis in this paper.     

An analysis of key characteristics of the Frank-Read source process in FCC metals

A well-known intragranular dislocation source, the Frank-Read (FR) source plays an important role in size-dependent dislocation multiplication in crystalline materials. Despite a number of studies in this topic, a systematic investigation of multiple aspects of the FR source in different materials is lacking. In this paper, we employ large scale quasistatic concurrent atomistic-continuum (CAC) simulations to model an edge dislocation bowing out from an FR source in Cu, Ni, and Al. First, a number of quantities that are important for the FR source process are quantified in the coarse-grained domain. Then, two key characteristics of the FR source, including the critical shear stress and critical dislocation configuration, are investigated. In all crystalline materials, the critical stresses and the aspect ratio of the dislocation half-loop height to the FR source length scale well with respect to the FR source length. In Al, the critical stress calculated by CAC simulations for a given FR source length agrees reasonably well with a continuum model that explicitly includes the dislocation core energy. Nevertheless, the predictions of the isotropic elastic theory do not accurately capture the FR source responses in Cu and Ni, which have a relatively large stacking fault width and elastic anisotropy. Our results highlight the significance of directly simulating the FR source activities using fully 3D models and shed light on developing more accurate continuum models.