Strain gradient plasticity,strengthening effects and plastic limit analysis

Within the framework of isotropic strain gradient plasticity, a rate-independent constitutive model exhibiting size dependent hardening is formulated and discussed with particular concern to its strengthening behavior. The latter is modelled as a (fictitious) isotropic hardening featured by a potential which is a positively degree-one homogeneous function of the effective plastic strain and its gradient. This potential leads to a strengthening law in which the strengthening stress, i.e. the increase of the plastically undeformed material initial yield stress, is related to the effective plastic strain through a second order PDE and related higher order boundary conditions. The plasticity flow laws, with the role there played by the strengthening stress, are addressed and shown to admit a maximum dissipation principle. For an idealized elastic perfectly plastic material with strengthening effects, the plastic collapse load problem of a micro/nano scale structure is addressed and its basic features under the light of classical plastic limit analysis are pointed out. It is found that the conceptual framework of classical limit analysis, including the notion of rigid-plastic behavior, remains valid. The lower bound and upper bound theorems of classical limit analysis are extended to strengthening materials. A static-type maximum principle and a kinematic-type minimum principle, consequences of the lower and upper bound theorems, respectively, are each independently shown to solve the collapse load problem. These principles coincide with their respective classical counterparts in the case of simple material. Comparisons with existing theories are provided. An application of this nonclassical plastic limit analysis to a simple shear model is also presented, in which the plastic collapse load is shown to increase with the decreasing sample size (Hall–Petch size effects).     

Shakedown analysis for a class of strengthening materials within the framework of gradient plasticity

The classical shakedown theory is extended to a class of perfectly plastic materials with strengthening effects (Hall–Petch effects). To this aim, a strain gradient plasticity model previously advanced by Polizzotto (2010) is used, whereby a featuring strengthening law provides the strengthening stress, i.e. the increase of the yield strength produced by plastic deformation, as a degree-zero homogeneous second-order differential form in the accumulated plastic strain with associated higher order boundary conditions. The extended static (Melan) and kinematic (Koiter) shakedown theorems are proved together with the related lower bound and upper bound theorems. The shakedown limit load problem is addressed and discussed in the present context, and its solution uniqueness shown out. A simple micro-scale structural system is considered as an illustrative example. The shakedown limit load is shown to increase with decreasing the structural size, which is a manifestation of the classical Hall–Petch effects in a context of cyclic loading.     

Size-dependent bending of thin metallic films

Size-dependent large curvature pure bending of thin metallic films has been analytically studied taking into account the associated strengthening mechanisms at different thickness scales. The classical plasticity theory is applicable to films thicker than 100 μm. Consequently, their bending capacity is governed by the competition between the material hardening and the thickness reduction. For films with a thickness ranging from fractions of a micron to a few microns, in addition to the above mechanisms, the strain gradient effect plays an important role and introduces an internal length scale. When the film thickness reduces to the nano-scale, the strain gradient effect is gradually replaced by the dominant surface stress/energy effect.     

Modelling the torsion of thin metal wires by distortion gradient plasticity

Under small strains and rotations, we apply a phenomenological higher-order theory of distortion gradient plasticity to the torsion problem, here assumed as a paradigmatic benchmark of small-scale plasticity. Peculiar of the studied theory, proposed about ten years ago by Morton E. Gurtin, is the constitutive inclusion of the plastic spin, affecting both the free energy and the dissipation. In particular, the part of the free energy, called the defect energy, which accounts for Geometrically Necessary Dislocations, is a function of Nye's dislocation density tensor, dependent on the plastic distortion, including the plastic spin. For the specific torsion problem, we implement this distortion gradient plasticity theory into a Finite Element (FE) code characterised by implicit (Backward Euler) time integration, numerically robust and accurate for both viscoplastic and rate-independent material responses. We show that, contrariwise to other higher-order theories of strain gradient plasticity (neglecting the plastic spin), the distortion gradient plasticity can predict some strengthening even if a quadratic defect energy is chosen. On the basis of the results of many FE analyses, concerned with (i) cyclic loading, (ii) switch in the higher-order boundary conditions during monotonic plastic loading, (iii) the use of non-quadratic defect energies, and (iv) the prediction of experimental data, we mainly show that (a) including the plastic spin contribution in a gradient plasticity theory is highly recommendable to model small-scale plasticity, (b) less-than-quadratic defect energies may help in describing the experimental results, but they may lead to anomalous cyclic behaviour, and (c) dissipative (unrecoverable) higher-order finite stresses are responsible for an unexpected mechanical response under non-proportional loading.     

A unified treatment of strain gradient plasticity

A theoretical framework is presented that has potential to cover a large range of strain gradient plasticity effects in isotropic materials. Both incremental plasticity and viscoplasticity models are presented. Many of the alternative models that have been presented in the literature are included as special cases. Based on the expression for plastic dissipation, it is in accordance with Gurtin (J. Mech. Phys. Solids 48 (2000) 989; Int. J. Plast. 19 (2003) 47) argued that the plastic flow direction is governed by a microstress q_ij and not the deviatoric Cauchy stress σ_ij′ that has been assumed by many others. The structure of the governing equations is of second order in the displacements and the plastic strains which makes it comparatively easy to implement in a finite element programme. In addition, a framework for the formulation of consistent boundary conditions is presented. It is shown that there is a close connection between surface energy of an interface and boundary conditions in terms of plastic strains and moment stresses. This should make it possible to study boundary layer effects at the interface between grains or phases. Consistent boundary conditions for an expanding elastic-plastic boundary are as well formulated. As examples, biaxial tension of a thin film on a thick substrate, torsion of a thin wire and a spherical void under remote hydrostatic tension are investigated.     

Bending and buckling of thin strain gradient elastic beams

Bending of strain gradient elastic thin beams is studied adopting Bernoulli-Euler principle. Simple linear strain gradient elastic theory with surface energy is employed. The governing beam equations with its boundary conditions are derived through a variational method. It turns out that new terms are introduced, indicating the importance of the cross-section area in bending of thin beams. Those terms are missing from the existing strain gradient beam theories. Those terms increase highly the stiffness of the thin beam. The buckling problem of the thin beams is also discussed.     

Nonlinear strain gradient elastic thin shallow shells

The governing equilibrium equations for strain gradient elastic thin shallow shells are derived, considering nonlinear strains and linear constitutive strain gradient elastic relations. Adopting Kirchhoff’s theory of thin shallow structures, the equilibrium equations, along with the boundary conditions, are formulated through a variational procedure. It turns out that new terms are introduced, indicating the importance of the cross-section area in bending of thin plates. Those terms are missing from the existing strain gradient shallow thin shell theories. Those terms highly increase the stiffness of the structures. When the curvature of the shallow shell becomes zero, the governing equilibrium for the plates is derived.     

Size-dependent yield strength of thin films

Biaxial strain and pure shear of a thin film are analysed using a strain gradient plasticity theory presented by Gudmundson [Gudmundson, P., 2004. A unified treatment of strain gradient plasticity. Journal of the Mechanics and Physics of Solids 52, 1379–1406]. Constitutive equations are formulated based on the assumption that the free energy only depends on the elastic strain and that the dissipation is influenced by the plastic strain gradients. The three material length scale parameters controlling the gradient effects in a general case are here represented by a single one. Boundary conditions for plastic strains are formulated in terms of a surface energy that represents dislocation buildup at an elastic/plastic interface. This implies constrained plastic flow at the interface and it enables the simulation of interfaces with different constitutive properties. The surface energy is also controlled by a single length scale parameter, which together with the material length scale defines a particular material.Numerical results reveal that a boundary layer is developed in the film for both biaxial and shear loading, giving rise to size effects. The size effects are strongly connected to the buildup of surface energy at the interface. If the interface length scale is small, the size effect vanishes. For a stiffer interface, corresponding to a non-vanishing surface energy at the interface, the yield strength is found to scale with the inverse of film thickness.Numerical predictions by the theory are compared to different experimental data and to dislocation dynamics simulations. Estimates of material length scale parameters are presented.     

Characterizing Torsional Properties of Microwires Using an Automated Torsion Balance

In order to characterize the torsional behavior of microwires, an automated torsion tester is established based on the principle of torsion balance. The main challenges in developing a torsion tester at small scales are addressed. An in-situ torsional vibration method for precisely calibrating the torque meter is developed. The torsion tester permits the measurement of torque to nN m, as a function of surface shear strain to a sensitivity of sub-microstrain. Using this technique, we performed (monotonic and/or cyclic) torsion tests on polycrystalline copper and gold wires. It is found that (i) a size effect appears in both the initial yielding and the plastic flow of torsional response; (ii) a reverse plasticity occurs upon unloading in cyclic torsion response; and (iii) the Hall-Petch effect and the strain gradient effect are synergistic. We also performed cyclic torsion tests on human hairs and spider silk which are natural protein fibers with a different morphological structure to metallic wires. It is shown that the single hair exhibits torsional recovery, and that the spider silk displays torsionally superelastic behavior whereby it is able to withstand great shear strain.     

Determination of material intrinsic length and strain gradient hardening in microbending process

Microbending experiments of pure aluminum show that the springback angles increase with the decrease of foil thickness, which indicates obvious size effects and attributes to plastic strain gradient hardening. Then a constitutive model, taking into accounts both plastic strain and plastic strain gradient hardening, is proposed to analyze the microbending process of thin foil. The model is based on the relationship between shear yield stress and dislocation density, in which the material intrinsic length is related to material properties and average grain numbers along the characteristic scale direction of part. It is adopted in analytical model to calculate the non-dimensional bending moment and predict the springback angle after microbending. It is confirmed that the predictions by the proposed hardening model agree well with the experimental data, while those predicted by the classical plasticity model cannot capture such size effects. The contribution of plastic strain gradient increases with the decrease of foil thickness and is independent on the bending angle.     

On bending of strain gradient elastic micro-plates

Bending of strain gradient elastic thin plates is studied, adopting Kirchhoff’s theory of plates. Simple linear strain gradient elastic theory with surface energy is employed. The governing plate equation with its boundary conditions are derived through a variational method. It turns out that new terms are introduced, indicating the importance of the cross-section area in bending of thin plates. Those terms are missing from the existing strain gradient plate theories; however, they strongly increase the stiffness of the thin plate.     

Formulation of strain gradient plasticity with interface energy in a consistent thermodynamic framework

In this work, the strain gradient formulation is used within the context of the thermodynamic principle, internal state variables, and thermodynamic and dissipation potentials. These in turn provide balance of momentum, boundary conditions, yield condition and flow rule, and free energy and dissipative energies. This new formulation contributes to the following important related issues: (i) the effects of interface energy that are incorporated into the formulation to address various boundary conditions, strengthening and formation of the boundary layers, (ii) nonlocal temperature effects that are crucial, for instance, for simulating the behavior of high speed machining for metallic materials using the strain gradient plasticity models, (iii) a new form of the nonlocal flow rule, (iv) physical bases of the length scale parameter and its identification using nano-indentation experiments and (v) a wide range of applications of the theory. Applications to thin films on thick substrates for various loading conditions and torsion of thin wires are investigated here along with the appropriate length scale effect on the behavior of these structures. Numerical issues of the theory are discussed and results are obtained using Matlab and Mathematica for the nonlinear ordinary differential equations (NODE) which constitute the boundary value problem.This study reveals that the micro-stress term has an important effect on the development of the boundary layers and hardening of the material at both hard and soft interface boundary conditions in thin films. The interface boundary conditions are described by the interfacial length scale and interfacial strength parameters. These parameters are important to control the size effect and hardening of the material. For more complex geometries the generalized form of the boundary value problem using the nonlocal finite element formulation is required to address the problems involved.     

Experimental techniques for testing unstiffened plates in compression and bending

We present details of a dual-actuator rig developed for testing rectangular plates supported on three sides, with the remaining (longitudinal) edge free, under combined uniaxial compression and in-plane bending. Particular attention is given to ensuring a constant strain gradient at the loaded ends, as opposed to a constant load eccentricity, in order to determine the post-buckling behavior and ultimate load and moment capacities of unstiffened thin-walled elements. Strain gradients varying from pure compression to pure bending are facilitated. Attention is also given to ensuring simply supported boundary conditions, and the methods used for anchoring the tensile stresses that develop at the loaded edges as a result of large plate deflections. Details of the methods for controlling the applied displacements are given, for which a system of four laser displacement devices was employed in order to achieve the required strain gradient. The operation of the rig is verified against established theoretical solutions.     

A nonlocal strain gradient plasticity theory for finite deformations

Strain gradient plasticity for finite deformations is addressed within the framework of nonlocal continuum thermodynamics, featured by the concepts of (nonlocality) energy residual and globally simple material. The plastic strain gradient is assumed to be physically meaningful in the domain of particle isoclinic configurations (with the director vector triad constant both in space and time), whereas the objective notion of corotational gradient makes it possible to compute the plastic strain gradient in any domain of particle intermediate configurations. A phenomenological elastic–plastic constitutive model is presented, with mixed kinematic/isotropic hardening laws in the form of PDEs and related higher order boundary conditions (including those associated with the moving elastic/plastic boundary). Two fourth-order projection tensor operators, functions of the elastic and plastic strain states, are shown to relate the skew-symmetric parts of the Mandel stress and back stress to the related symmetric parts. Consistent with the thermodynamic restrictions therein derived, the flow laws for rate-independent associative plasticity are formulated in a six-dimensional tensor space in terms of symmetric parts of Mandel stresses and related work-conjugate generalized plastic strain rates. A simple shear problem application is presented for illustrative purposes.     

An evaluation of higher-order plasticity theories for predicting size effects and localisation

Conventional plasticity theories are unable to capture the observed increase in strength of metallic structures with diminishing size. They also give rise to ill-posed boundary value problems at the onset of material softening. In order to overcome both deficiencies, a range of higher-order plasticity theories have been formulated in the literature. The purpose of this paper is to compare existing higher-order theories for the prediction of a size effect and the handling of localisation effects. To this end, size effect predictions for foils in bending are compared with existing experimental data. Furthermore, a study of one-dimensional harmonic incremental solutions from a uniform reference state allows one to assess the nature of material localisation as predicted by these competing higher-order theories. These analyses show that only one of the theories considered—the Fleck–Hutchinson strain gradient plasticity theory based upon the Toupin–Mindlin strain gradient framework [Fleck, N.A., Hutchinson, J.W., 1997. Strain gradient plasticity. Adv. Appl. Mech. 33, 295–361]—allows one to describe both phenomena. The other theories show either nonphysical size effects or a pathologically localised post-peak response.     

Wrinkling of tubes in bending from finite strain three-dimensional continuum theory

The problem of a tube under pure bending is first solved as a generalised plane strain problem. This then provides the prebifurcation solution, which is uniform along the length of the tube. The onset of wrinkling is then predicted by introducing buckling modes involving a sinusoidal variation of the displacements along the length of the tube. Both the prebuckling analysis and the bifurcation check require only a two-dimensional finite element discretisation of the cross-section with special elements. The formulation does not rely on any of the approximations of a shell theory, or small strains. The same elements can be used for pure bending and local buckling a prismatic beam of arbitrary cross-section. Here the flow theory of plasticity with isotropic hardening is used for the prebuckling solution, but the bifurcation check is based on the incremental moduli of a finite strain deformation theory of plasticity.For tubes under pure bending, the results for limit point collapse (due to ovalisation) and bifurcation buckling (wrinkling) are compared to existing analysis and test results, to see whether removing the approximations of a shell theory and small strains (used in the existing analyses) leads to a better prediction of the experimental results. The small strain analysis results depend on whether the true or nominal stress–strain curve is used. By comparing small and finite strain analysis results it is found that the small strain approximation is good if one uses (a) the nominal stress–strain curve in compression to predict bifurcation buckling (wrinkling), and (b) the true stress–strain curve to calculate the limit point collapse curvature.In regard to the shell theory approximations, it is found that the three-dimensional continuum theory predicts slightly shorter critical wrinkling wavelengths, especially for lower diameter-to-wall-thickness (D/t) ratios. However this difference is not sufficient to account for the significantly lower wavelengths observed in the tests.     

Strain gradient plasticity analysis of transformation induced plasticity in multiphase steels

“To what extent do plastic strain gradients affect the strengthening resulting from the transformation of small metastable inclusions into hard inclusions within a plastically deforming matrix?” is the central question addressed here. Though general in the approach, the focus is on the behavior of TRIP-assisted multiphase steels. A two-dimensional embedded cell model of a simplified microstructure composed of a single metastable austenitic inclusion surrounded by a soft ferritic matrix is considered. The cell is inserted in a large homogenized medium. The transformation of a fraction of the austenite into a hard martensite plate is simulated, accounting for a transformation strain, and leading to complex elastic and plastic accommodation. The size of a transforming plate in real multiphase steels is typically between 0.1 and 2 μm, a range of size in which plastic strain gradient effects are expected to play a major role. The single parameter version of the Fleck–Hutchinson strain gradient plasticity theory is used to describe the plasticity in the austenite, ferrite and martensite phases. The higher order boundary conditions imposed on the plastic flow have a large impact on the predicted strengthening. Using realistic values of the intrinsic length parameter setting the scale at which the gradients effects have an influence leads to a noticeable increase of the strengthening on top of the increase due to the transformation of a volume fraction of the retained austenite. The geometrical parameters such as the volume fraction of retained austenite and of the transforming zone also bring significant strengthening. Strain gradient effects also significantly affect the stress state inside the martensite plate during and after transformation with a potential impact on the damage resistance of these steels.     

面内变刚度薄板弯曲问题的挠度?弯矩耦合神经网络方法

发展了一种求解面内变刚度功能梯度薄板弯曲问题的神经网络方法. 面内变刚度薄板弯曲问题的偏微分控制方程为一复杂的4阶偏微分方程, 传统的基于强形式的神经网络解法在求解该偏微分方程时可能会遇到难以收敛、边界条件难以处理的情况. 本文基于Kirchhoff薄板弯曲理论, 提出了一种直角坐标系下任意面内变刚度薄板弯曲问题的神经网络解法. 神经网络模型包含挠度网络与弯矩网络, 分别用于预测薄板的挠度与弯矩, 从而将求解4阶偏微分方程转换为求解一系列二阶偏微分方程组, 通过对挠度、弯矩试函数的构造可使得神经网络计算结果严格满足边界条件. 在误差的反向传播中, 根据本文提出的误差函数公式计算训练误差并结合Adam优化算法更新模型的内部参数. 求解了不同边界条件、形状的面内变刚度薄板弯曲问题, 并将所得计算结果与理论解、有限元解进行对比. 研究表明, 本文模型对于求解面内变刚度薄板弯曲问题具备适应性, 虽然模型中的弯矩网络收敛较挠度网络要慢, 但本文方法在试函数的构造上更为简单、适应性更强.      

A thermodynamic based higher-order gradient theory for size dependent plasticity

A physically motivated and thermodynamically consistent formulation of small strain higher-order gradient plasticity theory is presented. Based on dislocation mechanics interpretations, gradients of variables associated with kinematic and isotropic hardenings are introduced. This framework is a two non-local parameter framework that takes into consideration large variations in the plastic strain tensor and large variations in the plasticity history variable; the equivalent (effective) plastic strain. The presence of plastic strain gradients is motivated by the evolution of dislocation density tensor that results from non-vanishing net Burgers vector and, hence, incorporating additional kinematic hardening (anisotropy) effects through lattice incompatibility. The presence of gradients in the effective (scalar) plastic strain is motivated by the accumulation of geometrically necessary dislocations and, hence, incorporating additional isotropic hardening effects (i.e. strengthening). It is demonstrated that the non-local yield condition, flow rule, and non-zero microscopic boundary conditions can be derived directly from the principle of virtual power. It is also shown that the local Clausius–Duhem inequality does not hold for gradient-dependent material and, therefore, a non-local form should be adopted. The non-local Clausius–Duhem inequality has an additional term that results from microstructural long-range energy interchanges between the material points within the body. A detailed discussion on the physics and the application of proper microscopic boundary conditions, either on free surfaces, clamped surfaces, or intermediate constrained surfaces, is presented. It is shown that there is a close connection between interface/surface energy of an interface or free surface and the microscopic boundary conditions in terms of microtraction stresses. Some generalities and utility of this theory are discussed and comparisons with other gradient theories are given. Applications of the proposed theory for size effects in thin films are presented.     

A conventional theory of mechanism-based strain gradient plasticity

There exist two frameworks of strain gradient plasticity theories to model size effects observed at the micron and sub-micron scales in experiments. The first framework involves the higher-order stress and therefore requires extra boundary conditions, such as the theory of mechanism-based strain gradient (MSG) plasticity [J Mech Phys Solids 47 (1999) 1239; J Mech Phys Solids 48 (2000) 99; J Mater Res 15 (2000) 1786] established from the Taylor dislocation model. The other framework does not involve the higher-order stress, and the strain gradient effect come into play via the incremental plastic moduli. A conventional theory of mechanism-based strain gradient plasticity is established in this paper. It is also based on the Taylor dislocation model, but it does not involve the higher-order stress and therefore falls into the second strain gradient plasticity framework that preserves the structure of conventional plasticity theories. The plastic strain gradient appears only in the constitutive model, and the equilibrium equations and boundary conditions are the same as the conventional continuum theories. It is shown that the difference between this theory and the higher-order MSG plasticity theory based on the same dislocation model is only significant within a thin boundary layer of the solid.