A symmetric over-nonlocal microplane model M4 for fracture in concrete

This paper analyzes the effectiveness of a nonlocal integral-type formulation of a constitutive law such as microplane model M4 in which the yield limits soften as a function of the total strain for prediction of fracture propagation. For a correct regularization of the mathematical problems caused by the softening behavior, an “over-nonlocal” generalization of the type proposed by Vermeer and Brinkgreve [Vermeer, P.A., Brinkgreve, R.B.J., 1994. A new effective non-local strain measure for softening plasticity. In: Chambon, R., Desrues, J., Vardoulakis, I. (Eds.), Localization and Bifurcation Theory for Soil and Rocks, Balkema, Rotterdam, pp. 89–100.] is adopted. Moreover, the symmetric weight function, proposed by Borino et al. [Borino, G. Failla, B., Parrinello, F., 2003. A symmetric nonlocal damage theory. International Journal of Solids and Structure 40, 3621–3645.] for damage mechanics, is introduced for the calculation of the nonlocal averaging of the total strain upon which the yield limits depend. The capability of the proposed model for reproducing the stress and strain fields in the vicinity of a notch is also investigated. Finally, the symmetric over-nonlocal generalization of microplane model M4 has been applied for the simulation of a mixed-mode fracture test such as the four-point-shear test and the test of axial tension at constant shear force [Nooru-Mohamend, M.B., 1992. Mixed-mode fracture of concrete: an experimental approach. Doctoral Thesis Delft University of Thechnology, Delft, The Netherlands.]     

Nonlocal implicit gradient-enhanced elasto-plasticity for the modelling of softening behaviour

An improved gradient-enhanced approach for softening elasto-plasticity is proposed, which in essence is fully nonlocal, i.e. an equivalent integral nonlocal format exists. The method utilises a nonlocal field variable in its constitutive framework, but in contrast to the integral models computes this nonlocal field with a gradient formulation. This formulation is considered ‘implicit’ in the sense that it strictly incorporates the higher-order gradients of the local field variable indirectly, unlike the common (explicit) gradient approaches. Furthermore, this implicit gradient formulation constitutes an additional partial differential equation (PDE) of the Helmholtz type, which is solved in a coupled fashion with the standard equilibrium condition. Such an approach is particularly advantageous since it combines the long-range interactions of an integral (nonlocal) model with the computational efficiency of a gradient formulation. Although these implicit gradient approaches have been successfully applied within damage mechanics, e.g. for quasi-brittle materials, the first attempts were deficient for plasticity. On the basis of a thorough comparison of the gradient-enhancements for plasticity and damage this paper rephrases the problem, which leads to a formulation that overcomes most reported problems. The two-dimensional finite element implementation for geometrically linear plain strain problems is presented. One- and two-dimensional numerical examples demonstrate the ability of this method to numerically model irreversible deformations, accompanied by the intense localisation of deformation and softening up to complete failure.     

Over-nonlocal gradient enhanced plastic-damage model for concrete

Classical continuum models exhibit strong mesh dependency during softening. One method to regularize the problem is to introduce a length scale parameter via the nonlocal formulation. However, standard nonlocal enhancement (either by integral or gradient formulation) may serve only as a partial localization limiter for many material models. The “over-nonlocal” formulation, where the weight for the nonlocal value is greater than unity and the excesses compensated by assigning a negative weight to the local value, is able to fully regularize certain material models when standard nonlocal enhancement fails to do so. A plastic-damage model for concrete is formulated with this over-nonlocal enhancement via the gradient approach and the full regularizing capabilities demonstrated.     

Nonlocal microplane model with strain-softening yield limits

The paper deals with the problem of nonlocal generalization of constitutive models such as microplane model M4 for concrete, in which the yield limits, called stress–strain boundaries, are softening functions of the total strain. Such constitutive models call for a different nonlocal generalization than those for continuum damage mechanics, in which the total strain is reversible, or for plasticity, in which there is no memory of the initial state. In the proposed nonlocal formulation, the softening yield limit is a function of the spatially averaged nonlocal strains rather than the local strains, while the elastic strains are local. It is demonstrated analytically as well numerically that, with the proposed nonlocal model, the tensile stress across the strain localization band at very large strain does soften to zero and the cracking band retains a finite width even at very large tensile strain across the band only if one adopts an “over-nonlocal” generalization of the type proposed by Vermeer and Brinkgreve [In: Chambon, R., Desrues, J., Vardoulakis, I. (Eds.), Localisation and Bifurcation Theory for Soils and Rocks, Balkema, Rotterdam, 1994, p. 89] (and also used by Planas et al. [Basic issue of nonlocal models: uniaxial modeling, Tecnical Report 96-jp03, Departamento de Ciencia de Materiales, Universidad Politecnica de Madrid, Madrid, Spain, 1996], and by Strömberg and Ristinmaa [Comput. Meth. Appl. Mech. Eng. 136 (1996) 127]). Numerical finite element studies document the avoidance of spurious mesh sensitivity and mesh orientation bias, and demonstrate objectivity and size effect.     

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

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

Buckling and post-buckling of gradient and nonlocal plasticity columns experiencing softening

The buckling and the post-buckling behaviors of a perfect axially loaded column are analytically investigated through a global bilinear moment–curvature elastoplastic constitutive law. Three plasticity cases are studied, namely the linear hardening plasticity law, the perfect elastoplastic case and the softening case. The applications of such a study can be found in various structural engineering problems, including reinforced concrete, steel, timber or composite structures. It is analytically shown that for all kinds of elastoplastic behaviors, the plasticity phenomena lead to a global softening branch in the load–deflection diagram. The propagation of the plasticity zone during the post-buckling process is analytically characterized in case of linear hardening or softening plasticity laws. However, it is shown that the unphysical elastic unloading solution necessarily occurs in presence of local softening moment–curvature constitutive law. A nonlocal plasticity moment–curvature softening law is then used to control the localization branch in the post-buckling stage. This nonlocal plasticity law includes the explicit and the implicit gradient plasticity law. Higher-order plasticity boundary conditions are derived from an extended variational principle. Some parametric studies finally illustrate the main findings of this paper, including the plasticity modulus effect on the post-buckling behavior of these plasticity structural systems.     

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.     

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

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

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.     

Gradient-enhanced softening material models

Classical constitutive models exhibit strong mesh dependency during softening and the numerical responses tend towards perfectly brittle behavior upon mesh refinements. Such sensitivity can be avoided by adopting the gradient-enhanced formulation. The implicit approach incorporates the gradient contributions indirectly via an additional Helmholtz equation and requires only C⁰ continuity. The explicit approach computes the gradient terms directly from the local field variables. Assuming a weak satisfaction of the yield function, C¹ continuity or C⁰ continuity with additional degrees of freedoms in the penalty approach is required. This makes the explicit method less attractive computationally. However, the explicit approach is able to fully regularize some material models where the standard implicit method fails to perform. Drawing analogy to the over-nonlocal integral formulation, the over-implicit-gradient framework is proposed. In addition, an alternative framework for the explicit gradient method requiring only C⁰ continuity is proposed. The regularizing effects of the abovementioned two gradient frameworks show promising applications to strain-softening materials.     

A nonlocal model with strain-based damage

A thermodynamically consistent formulation of nonlocal damage in the framework of the internal variable theories of inelastic behaviours of associative type is presented. The damage behaviour is defined in the strain space and the effective stress turns out to be additively splitted in the actual stress and in the nonlocal counterpart of the relaxation stress related to damage phenomena. An important advantage of models with strain-based loading functions and explicit damage evolution laws is that the stress corresponding to a given strain can be evaluated directly without any need for solving a nonlinear system of equations. A mixed nonlocal variational formulation in the complete set of state variables is presented and is specialized to a mixed two-field variational formulation. Hence a finite element procedure for the analysis of the elastic model with nonlocal damage is established on the basis of the proposed two-field variational formulation. Two examples concerning a one-dimensional bar in simple tension and a two-dimensional notched plate are addressed. No mesh dependence or boundary effects are apparent.     

A symmetric nonlocal damage theory

The paper presents a thermodynamically consistent formulation for nonlocal damage models. Nonlocal models have been recognized as a theoretically clean and computationally efficient approach to overcome the shortcomings arising in continuum media with softening. The main features of the presented formulation are: (i) relations derived by the free energy potential fully complying with nonlocal thermodynamic principles; (ii) nonlocal integral operator which is self-adjoint at every point of the solid, including zones near to the solid’s boundary; (iii) capacity of regularizing the softening ill-posed continuum problem, restoring a meaningful nonlocal boundary value problem. In the present approach the nonlocal integral operator is applied consistently to the damage variable and to its thermodynamic conjugate force, i.e. nonlocality is restricted to internal variables only. The present model, when associative nonlocal damage flow rules are assumed, allows the derivation of the continuum tangent moduli tensor and the consistent tangent stiffness matrix which are symmetric. The formulation has been compared with other available nonlocal damage theories.Finally, the theory has been implemented in a finite element program and the numerical results obtained for 1-D and 2-D problems show its capability to reproduce in every circumstance a physical meaningful solution and fully mesh independent results.     

Localization properties of strain-softening gradient plasticity models. Part II: Theories with gradients of internal variables

The second part of this paper compares and evaluates enhancements of the conventional plasticity theory by gradients of internal variables. Attention is focused on their performance as localization limiters. Both explicit and implicit gradient formulations are considered. It is shown that certain models suffer by serious mathematical deficiencies that would complicate their numerical implementation. Some other models are appropriate only at early stages of the softening process but later exhibit locking accompanied by a spurious expansion of the localized plastic zone. The comparative study indicates that a convenient and robust tool for regularized modeling of the entire localization process is provided by the implicit gradient approach combined with a suitable form of the hardening/softening law.     

Bending of an elastoplastic Hencky bar-chain: from discrete to nonlocal continuous beam models

The static behavior of an elastoplastic one-dimensional lattice system in bending, also called a microstructured elastoplastic beam or elastoplastic Hencky bar-chain (HBC) system, is investigated. The lattice beam is loaded by concentrated or distributed transverse monotonic forces up to the complete collapse. The phenomenon of softening localization is also included. The lattice system is composed of piecewise linear hardening–softening elastoplastic hinges connected via rigid elements. This physical system can be viewed as the generalization of the elastic HBC model to the nonlinear elastoplasticity range. This lattice problem is demonstrated to be equivalent to the finite difference formulation of a continuous elastoplastic beam in bending. Solutions to the lattice problem may be obtained from the resolution of piecewise linear difference equations. A continuous nonlocal elastoplastic theory is then built from the lattice difference equations using a continualization process. The new nonlocal elastoplastic theory associated with both a distributed nonlocal elastoplastic law coupled to a cohesive elastoplastic model depends on length scales calibrated from the spacing of the lattice model. Differential equations of the nonlocal engineering model are solved for the structural configurations investigated in the lattice problem. It is shown that the new micromechanics-based nonlocal elastoplastic beam model efficiently captures the scale effects of the elastoplastic lattice model, used as the reference. The hardening–softening localization process of the nonlocal continuous model strongly depends on the lattice spacing which controls the size of the nonlocal length scales.     

Modeling fracture in large scale shell structures

We consider a problem of modeling fracture and failure preceded by large scale yielding of ductile shells from the point of view of large-scale structural analysis. We place a special emphasis on the computational efficiency of the constitutive formulation. In this context, we seek the formulation embedded in the shell mechanics framework, which is both theoretically sound and easily implementable into a large-scale explicit dynamic finite element code without precluding vectorization or parallelization. This is achieved through the elasto-plastic damage constitutive model for finite-element analysis of plates and shells. The proposed damage model is purely phenomenological with a scalar damage parameter, which has no physical interpretation, except that it represents on a global scale the micromechanical changes the material undergoes during the process of necking and fracture. The localization leading to softening and fracture is represented by the damage calibration function with exponential damage growth after the onset of necking. The proposed phenomenological damage model uses a general plasticity and shell mechanics frameworks which makes it general and easily implementable into existing finite element codes. The proposed formulation has been implemented into the explicit dynamic finite element software code EPSA (Atkatsh et al., 1980, Atkatsh et al., 1983).     

On coupled gradient-dependent plasticity and damage theories with a view to localization analysis

Combinations of gradient plasticity with scalar damage and of gradient damage with isotropic plasticity are proposed and implemented within a consistently linearized format. Both constitutive models incorporate a Laplacian of a strain measure and an internal length parameter associated with it, which makes them suitable for localization analysis.The theories are used for finite element simulations of localization in a one-dimensional model problem. The physical relevance of coupling hardening/softening plasticity with damage governed by different damage evolution functions is discussed. The sensitivity of the results with respect to the discretization and to some model parameters is analyzed. The model which combines gradient-damage with hardening plasticity is used to predict fracture mechanisms in a Compact Tension test.     

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.     

Nonlocal gradient-dependent modeling of plasticity with anisotropic hardening

This work addresses the formulation of the thermodynamics of nonlocal plasticity using the gradient theory. The formulation is based on the nonlocality energy residual introduced by Eringen and Edelen (1972). Gradients are introduced for those variables associated with isotropic and kinematic hardening. The formulation applies to small strain gradient plasticity and makes use of the evanescent memory model for kinematic hardening. This is accomplished using the kinematic flux evolution as developed by Zbib and Aifantis (1988). Therefore, the present theory is a four nonlocal parameter-based theory that accounts for the influence of large variations in the plastic strain, accumulated plastic strain, accumulated plastic strain gradients, and the micromechanical evolution of the kinematic flux. Using the principle of virtual power and the laws of thermodynamics, thermodynamically-consistent equations are derived for the nonlocal plasticity yield criterion and associated flow rule. The presence of higher-order gradients in the plastic strain is shown to enhance a corresponding history variable which arises from the accumulation of the plastic strain gradients. Furthermore, anisotropy is introduced by plastic strain gradients in the form of kinematic hardening. Plastic strain gradients can be attributed to the net Burgers vector, while gradients in the accumulation of plastic strain are responsible for the introduction of isotropic hardening. The equilibrium between internal Cauchy stress and the microstresses conjugate to the higher-order gradients frames the yield criterion, which is obtained from the principle of virtual power. Microscopic boundary conditions, associated with plastic flow, are introduced to supplement the macroscopic boundary conditions of classical plasticity. The nonlocal formulation developed here preserves the classical assumption of local plasticity, wherein plastic flow direction is governed by the deviatoric Cauchy stress. The theory is applied to the problems of thin films on both soft and hard substrates. Numerical solutions are presented for bi-axial tension and simple shear loading of thin films on substrates.     

Nonlocal continuum theory of anisotropically damaged metals

The paper deals with a consistent and systematic general framework for the development of anisotropic continuum damage in ductile metals based on thermodynamic laws and nonlocal theories. The proposed model relies on finite strain kinematics based on the consideration of damaged as well as fictitious undamaged configurations related via metric transformation tensors which allow for the interpretation of damage tensors. The formulation is accomplished by rate-independent plasticity using a nonlocal yield condition of Drucker–Prager type, anisotropic damage based on a nonlocal damage growth criterion as well as non-associated flow and damage rules. The nonlocal theory of inelastic continua is established to be able to take into account long-range microstructural interaction. The approach incorporates macroscopic interstate variables and their higher-order gradients which properly describe the change in the internal structure and investigate the size effect of statistical inhomogeneity of the heterogeneous material. The idea of bridging length-scales is made by using higher-order gradients in the evolution equations of the equivalent inelastic strain measures which leads to a system of elliptic partial differential equations which is solved using the finite difference method at each iteration of the loading step and the displacement-based finite element procedure is governed by the standard principle of virtual work. Numerical simulations of the elastic–plastic deformation behavior of damaged solids demonstrate the efficiency of the formulation. Tension tests undergoing large strains are used to investigate the damage growth in high strength steel. The influence of various model parameters on the prediction of the deformation and localization of ductile metals is discussed.     

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.