基于能量等效原理的应变局部化分析:I.一维解析解

基于热力学第一定律和非局部塑性理论,提出了一种求解应变局部化问题的非局部方法.对材料的每一点定义了局部和非局部两种状态空间,局部状态空间的内变量通过非局部权函数映射到非局部空间,成为非局部内变量.在应变软化过程中,局部状态空间中的塑性变形服从正交流动法则,材料的软化律在非局部状态空间中被引入.通过两个状态空间的塑性应变能耗散率的等效,得到了应变软化过程中明确定义的局部化区域以及其中的塑性应变分布.应用本方法导出了一维应变局部化问题的解析解.解析解表明,应变局部化区域的尺寸只与材料内尺度有关;对于高斯型非局部权函数,局部化区域的尺寸大约是材料内尺度的6倍.一维算例表明,局部化区域的塑性应变分布以及载荷-位移曲线仅与材料参数和结构几何尺寸有关,变形局部化区域的尺寸随着材料内尺度的减小而减小,同时塑性应变也随着材料内尺度的减小变得更加集中.当内尺度趋近于零时,应用本文方法得到的解与采用传统的局部塑性理论得到的解相同.     

基于能量等效原理的应变局部化分析:Ⅱ.有限元解法

以非局部塑性理论为基础,应用状态空间理论,通过局部和非局部两个状态空间的塑性能量耗散率等效原理,提出了一种求解应变局部化问题的新方法,以得到与网格无关的数值解.针对二维问题的屈服函数和流动法则导出了求解非局部内变量的一般方程,并提出了在有限元环境中求解应变局部化问题的应力更新算法.为了验证所提出的方法,对1个一维拉杆和3个二维平面应变加载试件进行了有限元分析.数值结果表明,塑性应变的分布和载荷-位移曲线都随着网格的变小而稳定地收敛,应变局部化区域的尺寸只与材料内尺度有关,而对有限元网格的大小不敏感.对于一维问题,当有限元网格尺寸减小时,数值解收敛于解析解.对于二维剪切带局部化问题,数值解随着网格尺寸的减小而稳定地向唯一解收敛.当网格尺寸减小时,剪切带的宽度和方向基本上没有变化.而且得到的塑性应变分布和网格变形是平滑的.这说明,所提方法可以克服经典连续介质力学模型导致的网格相关性问题,从而获得具有物理意义的客观解.此模型只需要单元之间的位移插值函数具有C~0连续性,因而容易在现有的有限元程序中实现而无需对程序作大的修改.     

基于非局部塑性模型的应变局部化理论分析及数值模拟

通过求解一个第二类Fredholm方程,得到了基于非局部塑性软化模型的应变局部化问题理论解,结果表明,只有在当采用过非局部修正形式的非局部塑性软化模型才能得到应变局部化解,且得到的塑性应变分布和荷载响应依赖于所引入的特征长度及过非局部权参数。通过一维应变局部化有限元数值解,验证了非局部理论的引入能克服计算结果的网格敏感...     

岩石单轴压缩作用下变形局部化的梯度塑性解

采用梯度塑性理论研究单轴压缩作用下岩石变形局部化,得到了单轴压缩作用下岩石变形局部化带宽度的一维、二维解析解,为实验测定内部材料长度参数提供了理论依据.     

一个综合模糊裂纹和损伤的混凝土应变软化本构模型

本文研究就变软化材料的本构关系，提出了一个考虑损伤的粘塑性模型，损伤不仅影响材料的临界应力，而且影响材料的粘塑性，为模拟材料的应变软化行为，假设受损混凝土的破坏局部区域由模糊裂纹和损伤所统治，软化模量和局部区域尺度参量依赖于模糊裂纹扩展时释放的断裂能的参变量，用文中提出的模型计算了混凝土单轴压缩时不同应变速率下的瞬时应力应变响应以及等应力长期作用下的徐变，均得到很的结果。     

材料的率相关与梯度二、四阶耦合模型的内尺度律与稳定性研究

在率相关与梯度塑性二阶耦合本构模型的基础上,提出了二、四阶率相关与梯度塑性耦合模型。采用简谐波的分析方法对材料的应变局部化及材料的稳定性进行了研究,得到了二、四阶耦合模型在一维情况下的内尺度律的变化及材料稳定性的关系,得到了波长变化的上下界及材料稳定性的条件;并对其进行了对比性研究,得到材料稳定点移动的规律。     

断裂力学及边界元法在有限体应变局部化问题中的应用

微裂纹模型是研究砂土变形中应变局部化问题的微观力学模型。本文在前文的基础上建立了有限平面的徽裂纹模型的基本方程。数值解析表明,微裂纹模型可再现应变局部化及应变软化现象。同时,本文简要地讨论了边界约束、尺寸效应及侧压对应变局部化的影响。     

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

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

软化材料厚壁筒的解析解及其稳定性分析

将弹塑性材料的应力应变全过程曲线简化为三线性模型(弹性-线性软化-残余理想塑性),并假设材料服从Tresca屈服准则和关联流动法则,推导出受内压厚壁筒的解析解.在这个解析解的基础上,讨论了厚壁筒的平衡稳定性问题,内压达到临界载荷时,厚壁筒丧失稳定性,其临界载荷就是软化塑性材料厚壁筒的承载能力.     

薄板试样斜断口形成机制及损伤断裂过程分析

为分析不同材料和尺寸的薄板试样在室温下拉伸破坏后均形成与横截面夹角在20°~25°之间斜断口的原因,首先用统计方法对试样内随机分布微缺陷进行讨论,提出一种在宏观尺度上材料内微缺陷分布局部非均匀简化模型的假设.应用含孔材料损伤本构模型对含有不同方向微缺陷分布局部非均匀薄带区域的16MnNb薄板试样变形至破坏全过程进行数值模拟.结果表明,斜断口形成主要是由于试样内在与横截面夹角小于45°的带形区域内微缺陷分布局部非均匀造成,且与该带形区域在试样中位置无关;由于考虑微缺陷分布局部非均匀,得到试样的斜断口形成过程与试验现象完全一致;同时结合试验断口形貌,对变形过程中颈缩截面内损伤演化和破坏过程进行研究,进一步解释薄板试样的损伤破坏机制.     

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.     

Nonlocal and gradient rate plasticity

This paper deals with a formulation of nonlocal and gradient plasticity with internal variables. The constitutive model complies with local internal variables which govern kinematic hardening and isotropic softening and with a nonlocal corrective internal variable defined either as the sum between a new internal variable and its spatial weighted average or as the gradient of a measure of plastic strain. The rate constitutive problem is cast in the framework provided by the convex analysis and the potential theory for monotone multivalued operators which provide the suitable tools to perform a theoretical analysis of such nonlocal and gradient problems. The validity of the maximum dissipation theorem is assessed and constitutive variational formulations of the rate model are provided. The structural rate problem for an assigned load rate is then formulated. The related variational formulation in the complete set of state variable is contributed and the methodology to derive variational formulations, with different combinations of the state variables, is explicitly provided. In particular the generalization to the present nonlocal and gradient model of the principles of Prager–Hodge, Greenberg and Capurso–Maier is presented. Finally nonlocal variational formulations provided in the literature are derived as special cases of the proposed model.     

Spectral analysis of localization in nonlocal and over-nonlocal materials with softening plasticity or damage

The paper shows that spectral wave propagation analysis reveals in a simple and clear manner the effectiveness of various regularization techniques for softening materials, i.e., materials for which the yield limits soften as a function of the total strain. Both plasticity and damage models are considered. It is verified analytically in a simple way that the nonlocal integral-type model with degrading yield limit depending on the total strain works correctly if and only one adopts an unconventional nonlocal formulation introduced in 1994 by Vermeer and Brinkgreve (and in 1996 by Planas, and by Strömberg and Ristinmaa), which is here called, for the sake of brevity, ‘over-nonlocal’ because it uses a linear combination of local and nonlocal variables in which a negative weight imposed on the local variable is compensated by assigning to the nonlocal variable weight greater than 1 (this is equivalent to a nonlocal variable with a smooth positive weight function of total weight greater than 1, normalized by superposing a negative delta-function spike at the center). The spectral approach readily confirms that the nonlocal integral-type generalization of softening plasticity with an additive format gives correct localization properties only if an over-nonlocal formulation is adopted. By contrast, the nonlocal integral-type generalization of softening plasticity with a multiplicative format provides realistic localization behavior, just like the nonlocal integral-type damage model, and thus does not necessitate an over-nonlocal formulation. The localization behavior of explicit and implicit gradient-type models is also analyzed. A simple analysis shows that plasticity and damage models with gradient-type localization limiter, whether explicit or implicit, have very different localization behaviors.     

Novel variational formulations for nonlocal plasticity

A nonlocal structural model of softening plasticity is considered in the framework of the internal variable theories of inelastic behaviours of associative type. The finite-step nonlocal structural problem in a geometrically linear range is formulated according to a backward difference scheme for time integration of the flow rule. The related finite-step variational formulation in the complete set of local and nonlocal state variables is recovered. A family of mixed nonlocal variational formulations, with different combinations of state variables, is provided starting from the general variational formulation. The specialization of a mixed variational formulation to existing nonlocal models of softening plasticity, assuming both linear and nonlinear constitutive behaviour, is provided to show the effectiveness of the theory.     

Decaying Turbulence in the Generalised Burgers Equation

This work is concerned with the derivation of optimal scaling laws, in the sense of matching lower and upper bounds on the energy, for a solid undergoing ductile fracture. The specific problem considered concerns a material sample in the form of an infinite slab of finite thickness subjected to prescribed opening displacements on its two surfaces. The solid is assumed to obey deformation-theory of plasticity and, in order to further simplify the analysis, we assume isotropic rigid-plastic deformations with zero plastic spin. When hardening exponents are given values consistent with observation, the energy is found to exhibit sublinear growth. We regularize the energy through the addition of nonlocal energy terms of the strain-gradient plasticity type. This nonlocal regularization has the effect of introducing an intrinsic length scale into the energy. Under these assumptions, ductile fracture emerges as the net result of two competing effects: whereas the sublinear growth of the local energy promotes localization of deformation to failure planes, the nonlocal regularization stabilizes this process, thus resulting in an orderly progression towards failure and a well-defined specific fracture energy. The optimal scaling laws derived here show that ductile fracture results from localization of deformations to void sheets, and that it requires a well-defined energy per unit fracture area. In particular, fractal modes of fracture are ruled out under the assumptions of the analysis. The optimal scaling laws additionally show that ductile fracture is cohesive in nature, that is, it obeys a well-defined relation between tractions and opening displacements. Finally, the scaling laws supply a link between micromechanical properties and macroscopic fracture properties. In particular, they reveal the relative roles that surface energy and microplasticity play as contributors to the specific fracture energy of the material.     

A model for high temperature creep of single crystal superalloys based on nonlocal damage and viscoplastic material behavior

A model for high temperature creep of single crystal superalloys is developed, which includes constitutive laws for nonlocal damage and viscoplasticity. It is based on a variational formulation, employing potentials for free energy, and dissipation originating from plasticity and damage. Evolution equations for plastic strain and damage variables are derived from the well-established minimum principle for the dissipation potential. The model is capable of describing the different stages of creep in a unified way. Plastic deformation in superalloys incorporates the evolution of dislocation densities of the different phases present. It results in a time dependence of the creep rate in primary and secondary creep. Tertiary creep is taken into account by introducing local and nonlocal damage. Herein, the nonlocal one is included in order to model strain localization as well as to remove mesh dependence of finite element calculations. Numerical results and comparisons with experimental data of the single crystal superalloy LEK94 are shown.     

Suppressed plastic deformation at blunt crack-tips due to strain gradient effects

Large deformation gradients occur near a crack-tip and strain gradient dependent crack-tip deformation and stress fields are expected. Nevertheless, for material length scales much smaller than the scale of the deformation gradients, a conventional elastic–plastic solution is obtained. On the other hand, for significant large material length scales, a conventional elastic solution is obtained. This transition in behaviour is investigated based on a finite strain version of the Fleck–Hutchinson strain gradient plasticity model from 2001. The predictions show that for a wide range of material parameters, the transition from the conventional elastic–plastic to the elastic solution occurs for length scales ranging from 0.001 times the size of the plastic zone to a length scale of the same order of magnitude as the plastic zone.     

Modelling of localization and propagation of phase transformation in superelastic SMA by a gradient nonlocal approach

In this work, a nonlocal phenomenological behavior model is proposed in order to describe the localization and propagation of stress-induced martensite transformation in shape memory alloy (SMA) wires and thin films. It is a nonlocal extension of an existing local model that was derived from a micromechanical-inspired Gibbs free energy expression. The proposed model uses, besides the local field of the internal variable, namely the martensite volume fraction, a nonlocal counterpart. This latter acts as an additional degree of freedom, which is determined by solving an additional partial differential equation (PDE), derived so as to be equivalent to the integral definition of a nonlocal quantity. This PDE involves an internal length parameter, dictating the global scale at which the nonlocal interactions of the underlying micromechanisms are manifested during phase transformation. Moreover, to account for the unstable softening behavior, the transformation yield force parameter is considered as a gradually decreasing function of the martensite fraction. Possible material and geometric imperfections that are responsible for localization initiation are also considered in this analysis. The obtained constitutive equations are implemented in the Abaqus® finite element code in one and two dimensions. This requires the development of specific finite elements having the nonlocal volume fraction variable as an additional degree of freedom. This implementation is achieved through the UEL user’s subroutine. The effect of martensitic localization on the superelastic global behavior of SMA wire and holed thin plate, subjected to tension loading, is analyzed. Numerical results show that the developed tool correctly captures the commonly observed unstable superelastic behavior characterized by nucleation and propagation of martensitic phase. In particular, they show the influence of the internal length parameter, appearing in the nonlocal model, on the size of the localization area and the stress nucleation peak.