Modeling adiabatic shear failure from energy considerations

The numerical simulation of dynamic structural failure by localized shear is quite complex in terms of constitutive models and choice of adequate failure criteria, along with a pronounced mesh-sensitivity. As a result, the existing numerical procedures are usually quite sophisticated, so that their application for design purposes is still limited. This study is based on the implementation of a simple energy-based criterion, which was developed on experimental considerations (Rittel et al., 2006), and uses a minimal number of adjustable parameters. According to this criterion, a material point starts to fail when the total strain energy density reaches a critical value. Thereafter, the strength of the element decreases gradually to zero to mimic the actual structural behavior. The criterion was embedded into commercial finite element software and tested by simulating numerically four typical high-rate experiments. The first is the dynamic torsion test of a tubular specimen. The second concerns the failure mode transition in mode II fracture of an edge crack in plain strain. The last two involve dynamic shear localization under high rate compression of a cylindrical and a shear compression specimen. A very good adequation was found both qualitatively and quantitatively. Qualitatively, in terms of failure path selection, and quantitatively, in terms of local strains, temperatures and critical impact velocity. The proposed approach is enticing from an engineering perspective aimed at predicting the onset and propagation of dynamic shear localization in actual structures.     

Atomistic potentials based energy flux integral criterion for dynamic adiabatic shear banding

The energy flux integral criterion based on atomistic potentials within the framework of hyperelasticity–plasticity is proposed for dynamic adiabatic shear banding (ASB). System Helmholtz energy decomposition reveals that the dynamic influence on the integral path dependence is originated from the volumetric strain energy and partial deviatoric strain energy, and the plastic influence only from the rest part of deviatoric strain energy. The concept of critical shear banding energy is suggested for describing the initiation of ASB, which consists of the dynamic recrystallization (DRX) threshold energy and the thermal softening energy. The criterion directly relates energy flux to the basic physical processes that induce shear instability such as dislocation nucleations and multiplications, without introducing ad-hoc parameters in empirical constitutive models. It reduces to the classical path independent J-integral for quasi-static loading and elastic solids. The atomistic-to-continuum multiscale coupling method is used to simulate the initiation of ASB. Atomic configurations indicate that DRX induced microstructural softening may be essential to the dynamic shear localization and hence the initiation of ASB.     

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 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.     

A numerical methodology for investigating the formation of adiabatic shear bands

In this paper, we present a numerical approach for analyzing thermo-visco-plastic deformation in one dimension. The method, which is accurate to second order, is based on integration along the characteristic lines. It is able to simulate fully localized plastic flow with high resolution and good efficiency. We apply this numerical scheme to the analysis of shear localization, emphasizing the interactions between a single shear band and its surroundings and among the members of a periodic array of shear bands. It is found that a shear band may grow intermittently due to interactions with other bands. The developed method is specifically adequate for analyzing the self-organized multiple adiabatic formation process, which will be discussed in the follow-up paper.     

边坡稳定的剪切带计算

为了解决边坡稳定分析中剪切带有限元网格的依赖性问题,采用梯度塑性理论,从本构关系中引入特征长度入手,建立计算模型。提出了一种8节点缩减积分的梯度塑性单元,并采用梯度塑性理论推导了Drucker-Prager屈服准则的软化模型的有限元格式,在ABAQUS中进行了二次开发,嵌入了本文提出的8节点单元和本构模型,并用ABAQUS软件进行了边坡剪切带的计算。计算结果表明,本文提出的方法消除了经典有限元计算的网格依赖性问题,可以得到与单元剖分无关的稳定的剪切带宽度。本文所提出的方法可适用于其他场合的剪切带计算。     

On evolution of thermo-plastic shear band

The present paper briefly reviews analytical studies of the evolution of thermoplastic shear band, i.e. emergence from uniform deformation, post-instability growth and late stage behaviour. The case studied is the simple shear of temperature and rate-dependent materials with heat transfer. Uniform mode exists before a critical state, if no heat flows out of testpiece. Upon reaching the critical state, bifurcation appears as a result of disturbances, which leads to instability and the formation of narrow shear band. Initially, the band, due to temperature disturbance, can shrink with increasing temperature and strain rate owing to unsteady flow. Then heat conduction dominates and causes the shear band to expand. The postmortem appearance of thermo-plastic shear band manifests itself as balance of plastic work rate and heat diffusion. Melting may also take place within the band.     

Prediction of incipient shear band trajectories in a thick wall cylinder explosion test

A simple but precise and physical mechanism-based mathematical expression is proposed to predict shear band trajectories in thick wall cylinders subject to external surrounding pressure. The expression is based on the Coulomb-Mohr fracture criterion and can be applied to various compression-sensitive materials, especially ceramics. The predicted result closely matches the experimental observations, which makes this method quite useful in testing material behavior. This expression also permits the extraction of the parameter μ in the Coulomb-Mohr criterion from experimental observations. Furthermore, no pre-assumptions or after-test measurement are necessary in order to carry out the prediction. The only two values needed to conduct the prediction are the initial inner radius and the friction coefficient μ. A comparison between the newly proposed model and existing theory is made to reveal their relations and demonstrate the effectiveness of the newly derived mathematical expression.     

Adiabatic shear band localization of inelastic single crystals in symmetric double-slip process

Summary The main objective of the present paper is the development of a viscoplastic regularization procedure valid for an adiabatic dynamic process for multi-slips of single crystals. The next objective is to focus attention on the investigation of instability criteria, and particularly on shear band localization conditions.To achieve this aim, an analysis of acceleration waves is given, and advantage is taken of the notion of the instantaneous adiabatic acoustic tensor. If zero is an eigenvalue of the acoustic tensor, then the associated discontinuity does not propagate, and one speaks of a stationary discontinuity. This situation is referred to as the strain localization condition, and corresponds to a loss of hyperbolicity of the dynamical equations. It has been proved that for an, adiabatic process of rate-dependent (elastic-viscoplastic) crystal, the wave speed of discontinuity surface always remains real and different from zero. It means that for this case the initial-value problem is well-posed. However, for an adiabatic process of rate-independent(elastic-plastic) crystal, the wave speed of discontinuity surface can be equal zero. Then the necessary condition for a localized plastic deformation along the shear band to be formed is as follows: the determinant of the instantaneous adiabatic acoustic tensor is equal to zero. This condition for localization is equivalent to that obtained by using the standard bifurcation method. Based on this idea, the conditions for adiabatic shear band localization of plastic deformation have been investigated for single crystals. Particular attention has been focused on the discussion of the influence of thermal expansion, thermal plastic, softening and spatial covariance effects on shear band localization criteria for a planar model of an f.c.c. crystal undergoing symmetric primary-conjugate double slip. The results obtained have been compared with available experimental observations.Finally, it is noteworthy that the viscoplasticity regularization procedure can be used in the developing of an unconditionally stable numerical integration algorithm for simulation of adiabatic inelastic flow processes in ductile single crystals, cf. [21].The paper has been prepared within research programme sponsored by the Committee of Scientific Research under Grant 3 P404 031 07.     

Deformation localization and shear band fracture in strong anisotropy sheet tension

The tensile deformation localization and the shear band fracture behaviors of sheet metals with strong anisotropy are numerically simulated by using Updating Lagrange finite element method, Quasi-flow plastic constitutive theory^[1] and B-L planar anisotropy yield criterion^[2]. Simulated results are compared with experimental ones. Very good consistence is obtained between numerical and experimental results. The relationship between the anisotropy coefficientR and the shear band angle θ is found. The project supported by the National Natural Science Foundation of China and the Excellent Youth Teacher Foundation of the State Education Commission of China     

Near-tip out of plane displacement fields for dynamic crack propagation in functionally graded materials

Asymptotic expansion for the out of plane displacement field around a crack propagating along the gradient in a functionally graded material is developed. The irregular behavior of one of the terms in the expansion at low crack speeds is further examined and a remedial solution, which is well behaved at low crack speeds, is proposed. The developed out of plane displacement field is used to estimate stress intensity factor from quasi-static finite element solution. The results indicate that inclusion of the proposed nonhomogeneity specific terms gives estimates of stress intensity factor, which are consistent with existing analytical predictions.