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

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

土质边坡失稳的突变性分析

对土质边坡圆弧滑动稳定性问题，采用应力软化模型，由系统能量导出极限平衡方程，运用突变理论方法对土坡失稳进行了分析。分析表明，边坡参数变化为某一临界值时，土坡的滑动位移突然增大，发生突变失稳现象，由刚体极限平衡稳定理论无法体现应力软化土质边坡的这种失稳的突变性。     

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

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

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

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

基于偶应力理论剪切带问题的弹塑性有限元分析

对于软化材料的剪切带问题,传统弹塑性有限元分析遇到了困难,进入弹塑性阶段,计算结果对网格划分敏感,出现所谓的有限元网格依赖性问题,随着网格的细分,计算常常因不收敛导致失效. 用有限元软件ABAQUS计算了3个例题,证实了传统弹塑性有限元分析软化材料剪切带问题的局限性,同时证实对于无剪切带的厚壁筒问题不会出现上述问题. 进一步引入细观非局部化理论,对非局部理论含有的细观参数\ell进行了深入讨论,并采用可通过C0 -1分片检验的18参偶应力三角形单元,重新计算了3个例题,结果避免了上述问题,说明细观偶应力有限元尤其适用于分析剪切带问题.      

冲击载荷下钨合金圆台试件绝热剪切变形局部化的数值模

采用有限元计算编码ABAQUS模拟了钨合金圆台试件在冲击载荷下的变形和剪切局部化行为。计算采用二维轴对称应变条件下的绝热模型。钨合金的本构方程采用热粘塑性形式的Johnson Cook模型。为了得到不同尺度的变形信息 ,计算中用了两种网格 ;先用粗糙网格分析试件变形局部化的概貌 ;接着 ,用细密网格 (在变形局部化区域 ,网格尺寸达到 10 m)分析绝热剪切带的形成和发展。有限元模拟得到的绝热剪切带位置和方向与实验一致。计算结果表明 ,绝热剪切带的形成和发展与试件的应力状态密切相关。     

关于美国国家科学基金会主办的局部化塑性不稳定性和破坏准则讨论会的报告

<正> 在大塑性变形下,以狭窄剪切带形式表现的局部化不稳定性是必不可免的。这样的剪切带一旦形成,它们就要继续发展,以至带中的应变可以变得非常大。所以,剪切带是大塑性变形的一种重要的机理,同时它们又经常是最终韧性断裂的重要先兆和所在位置。自70年代初期以来,对剪切局部化现象已进行了广泛的研究,1989年11月14—15日,国家科学基金会(NSF)在华盛顿资助和组织了一个小型讨论会,以评估当前对“局部化塑性不稳定性      

关于美国国家科学基金会主办的局部化塑性不稳定性和破坏准则讨论会的报告

<正> 在大塑性变形下,以狭窄剪切带形式表现的局部化不稳定性是必不可免的。这样的剪切带一旦形成,它们就要继续发展,以至带中的应变可以变得非常大。所以,剪切带是大塑性变形的一种重要的机理,同时它们又经常是最终韧性断裂的重要先兆和所在位置。自70年代初期以来,对剪切局部化现象已进行了广泛的研究,1989年11月14—15日,国家科学基金会(NSF)在华盛顿资助和组织了一个小型讨论会,以评估当前对“局部化塑性不稳定性     

冲击载荷下钨合金圆台试件绝热剪切变形局部化的数值模拟

采用有限元计算编码ABAQUS模拟了钨合金圆台试件在冲击载荷下的变形和剪切局部化行为.计算采用二维轴对称应变条件下的绝热模型.钨合金的本构方程采用热粘塑性形式的Johnson-Cook模型.为了得到不同尺度的变形信息,计算中用了两种网格;先用粗糙网格分析试件变形局部化的概貌;接着,用细密网格(在变形局部化区域,网格尺寸达到10μm)分析绝热剪切带的形成和发展.有限元模拟得到的绝热剪切带位置和方向与实验一致.计算结果表明,绝热剪切带的形成和发展与试件的应力状态密切相关.     

稳定渗流条件下土坡稳定性分析的一种新方法

在提出随机角生成曲线滑动面方法的基础上,通过改变选定参数来搜索临界滑动面,对渗流作用下的不同土质边坡进行稳定性分析。 考虑渗流效应时,将土体视为三种土条划分情况,并进行受力分析,推导出简化Janbu法中土条的浮力和渗透力计算公式。同时,对实际渗流场做出合理简化,以便应用于计算机编程。算例对比分析表明: ①此方法计算出的最小安全系数与已有研究成果非常接近,得到的最危险滑动面也颇为相似,可证明本文方法的正确性; ②以非圆弧滑动面与圆弧滑动面计算出的最小安全系数相比,前者明显要小,建议在有渗流作用时应采用非圆弧滑动面; ③对两个不同水位下的均质土坡,以水位与坡高之比为参数对边坡状态进行分析时,结果表明随着地下水位的升高,边坡安全系数急剧下降,因而,渗流效应是边坡发生失稳的一个重要因素,工程实际中应特别注意边坡排水设施的设置。     

Strain gradient crystal plasticity effects on flow localization

In metal grains one of the most important failure mechanisms involves shear band localization. As the band width is small, the deformations are affected by material length scales. To study localization in single grains a rate-dependent crystal plasticity formulation for finite strains is presented for metals described by the reformulated Fleck–Hutchinson strain gradient plasticity theory. The theory is implemented numerically within a finite element framework using slip rate increments and displacement increments as state variables. The formulation reduces to the classical crystal plasticity theory in the absence of strain gradients. The model is used to study the effect of an internal material length scale on the localization of plastic flow in shear bands in a single crystal under plane strain tension. It is shown that the mesh sensitivity is removed when using the nonlocal material model considered. Furthermore, it is illustrated how different hardening functions affect the formation of shear bands.     

NONLINEAR UNIFIED STRENGTH MODEL OF GEOMATERIALS

将材料的破坏归结为剪切破坏,每种材料对应于特定的剪切滑动面,抗剪强度为滑动面上正应力的函数,基于不同材料的强度特性将一系列的剪切滑动面统一起来,建立了岩土材料的非线性统一强度模型.非线性统一强度模型的滑动面为β应力空间内的等倾面,在β应力空间内的强度面为圆锥面;在普通应力空间内的强度面为一系列连续光滑、外凸的锥面,在偏平面上强度曲线涵盖了从下限Matsuoka-Nakai曲线到上限Drucker-Prager圆之间的所有区域,子午面上强度线为直线.非线性统一强度模型只有3个材料参数,参数都具有明确的物理意义,通过与国内外学者已取得的岩土类材料真三轴强度试验结果的比较,表明模型可适用于多种类型的材料,并合理描述其非线性强度特性.     

岩土材料的非线性统一强度模型

将材料的破坏归结为剪切破坏,每种材料对应于特定的剪切滑动面,抗剪强度为滑动面上正应力的函数,基于不同材料的强度特性将一系列的剪切滑动面统一起来,建立了岩土材料的非线性统一强度模型.非线性统一强度模型的滑动面为β应力空间内的等倾面,在β应力空间内的强度面为圆锥面;在普通应力空间内的强度面为一系列连续光滑、外凸的锥面,在偏平面上强度曲线涵盖了从下限Matsuoka-Nakai曲线到上限Drucker-Prager圆之间的所有区域,子午面上强度线为直线.非线性统一强度模型只有3个材料参数,参数都具有明确的物理意义,通过与国内外学者已取得的岩土类材料真三轴强度试验结果的比较,表明模型可适用于多种类型的材料,并合理描述其非线性强度特性.      

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.     

框架预应力锚杆边坡支护结构的稳定性分析方法及其 应用

以边坡极限平衡理论及框架预应力锚杆边坡支护结构的圆弧滑动破坏模式为基础,在考虑框架、锚杆对土体边坡稳定性影响的情况下,基于圆弧滑动条分法的思想利用积分法建立了内部稳定性安全系数计算模型和最危险滑移面搜索模型。确定滑移面圆心坐标为几何控制参数,推导了安全系数计算模型中各变量与滑移面圆心坐标之间的函数关系式,从而获得了滑移面几何控制参数与内部稳定性安全系数之间的函数关系。利用网格法对框架预应力锚杆支护结构最危险滑移面的圆心坐标进行动态搜索和求解。最后采用Matlab语言编制了框架预应力锚杆边坡支护结构的稳定性分析程序,并结合工程实例进行了验算,讨论了此种方法的适用条件,结果表明该方法简明适用,较传统的经验分析方法更加合理。     

Critical state shear behavior of the soil-structure interface determined by discrete element modeling

The interface between soil and structure can be referred to as a soil-structure system, and its behavior plays an important role in many geotechnical engineering practices. In this study, results are presented from a series of monotonic direct shear tests performed on a sand-structure interface under constant normal stiffness using the discrete element method (DEM). Strain localization and dilatancy behavior of the interface is carefully examined at both macroscopic and microscopic scales. The effects of soil initial relative density and normal stress on the interface shear behavior are also investigated. The results show that a shear band progressively develops along the structural surface as shear displacement increases. At large shear displacement a unique relationship between stress ratio and void ratio is reached in the shear band for a certain normal stress, indicating that a critical state exists in the shear band. It is also found that the thickness and void ratio of the shear band at the critical state decreases with increasing normal stress. Comparison of the DEM simulation results with experimental results provides insight into the shear behavior of a sand-structure interface and offers a means for quantitative modeling of such interfaces based on the critical state soil mechanics.     

Analysis of shear band instabilities in sintered metals

For structural components made of sintered metals the porosity may be rather high, and this affects conditions for material failure. A simple model for analysis of shear band development in uniformly strained solids is used here to study the effect of rather high initial porosities on the onset of a material instability. The elastic plastic behaviour of the material is represented by a material model, which combines the Gurson model, relevant to rather low porosities, with the FKM model, developed for high porosity powder compacts. Predictions are shown for various levels of initial porosity and for various levels of initial material imperfections, considering both plane strain and axisymmetric conditions outside the shear band. Also a comparison of localization predictions by the Gurson model, the FKM model and the combined model is presented.     

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.     

Influence of void nucleation on ductile shear fracture at a free surface

An approximate continuum model of a ductile, porous material is used to study the influence of the nucleation and growth of micro-voids on the formation of shear bands and the occurrence of surface shear fracture in a solid subject to plane strain tension. Bifurcation into diffuse modes is analysed for a plane strain tensile specimen described by these constitutive relations, which account for a considerable plastic dilatancy due to void growth and for the possibility of non-normality of the plastic flow law. In particular, bifurcation into surface wave modes and the possible influence of such modes triggering shear bands is investigated. For solids with initial imperfactions such as a surface undulation, a local material inhomogeneity on an inclusion colony, the inception and growth of plastic flow localization is analysed numerically. Both the formation of void-sheets and the final growth of cracks in the shear bands is described numerically. Some special features of shear band development in the solid obeying non-normality are studied by a simple model problem.     

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.