The boundary layer solutions of the interface problem considering the strain gradient

According to the classical elastic theory, there is always a discontinuity of rotation angle on the interface between different materials. This illogic result can be overcome by the strain gradient plasticity theory. In the light of this theory, there is a group of boundary layer solutions near the interface, which have made important adjustment of the classical results. This project is supported by National Natural Science Foundation of China (19891180).     

理想正交各向异性弹塑性材料平面应变条件下静止裂纹尖端场

以Hill唯象理论为基础,建立正交各向异性弹塑性材料的本构关系,给出理想正交各向异性弹塑性材料在平面应变条件下混合型静止裂纹尖端的弹塑性场.与J.Pan的解不同,采用自相似假定,可以用解析方法求得不存在应力间断的应力场.对满塑性区条件和应变的奇异性加以讨论,这些为建立断裂准则提供了理论的依据.     

A mathematical basis for strain-gradient plasticity theory. Part II: Tensorial plastic multiplier

A phenomenological, flow theory version of gradient plasticity for isotropic and anisotropic solids is constructed along the lines of Gudmundson [Gudmundson, P., 2004. A unified treatment of strain-gradient plasticity. J. Mech. Phys. Solids 52, 1379-1406]. Both energetic and dissipative stresses are considered in order to develop a kinematic hardening theory, which in the absence of gradient terms reduces to conventional J₂ flow theory with kinematic hardening. The dissipative stress measures, work-conjugate to plastic strain and its gradient, satisfy a yield condition with associated plastic flow. The theory includes interfacial terms: elastic energy is stored and plastic work is dissipated at internal interfaces, and a yield surface is postulated for the work-conjugate stress quantities at the interface. Uniqueness and extremum principles are constructed for the solution of boundary value problems, for both the rate-dependent and the rate-independent cases. In the absence of strain gradient and interface effects, the minimum principles reduce to the classical extremum principles for a kinematically hardening elasto-plastic solid. A rigid-hardening version of the theory is also stated and the resulting theory gives rise to an extension to the classical limit load theorems. This has particular appeal as previous trial fields for limit load analysis can be used to generate immediately size-dependent bounds on limit loads.     

颗粒增强复合材料的界面开裂与尺度效应

采用基于Huang等提出的塑性应变梯度传统理论发展的有限元方法,模拟了颗粒增强金属基复合材料的界面开裂与颗粒尺度效应.分别针对考虑颗粒与基体间界面开裂和不开裂两种情况进行分析,并将考虑界面开裂的模拟结果与实验结果进行比较,证明了模型的有效性,同时也获得应变梯度理论中所包含的材料特征尺度参量的取值.     

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.     

Finite deformation analysis of mechanism-based strain gradient plasticity: torsion and crack tip field

A finite deformation theory of mechanism-based strain gradient (MSG) plasticity is developed in this paper based on the Taylor dislocation model. The theory ensures the proper decomposition of deformation in order to exclude the volumetric deformation from the strain gradient tensor since the latter represents the density of geometrically necessary dislocations. The solution for a thin cylinder under large torsion is obtained. The numerical method is used to investigate the finite deformation crack tip field in MSG plasticity. It is established that the stress level around a crack tip in MSG plasticity is significantly higher than its counterpart (i.e. HRR field) in classical plasticity.     

数据驱动梯度结构材料弹塑性本构

梯度结构材料因其优异的力学性能被广泛应用于工程结构中.论文整合塑性理论和人工神经网络技术,发展了一种构建梯度结构材料弹塑性本构模型的新方法.该方法基于梯度结构材料不同位置的微结构,构建不同代表性体积单元,进而生成应力应变数据,应用生成的数据训练人工神经网络,建立基于神经网络的材料本构模型.应用该方法,论文开展了针对实际...     

Post-necking behaviour modelled by a gradient dependent plasticity theory

The delay of the onset of localization and the post-necking behaviour for stretched thin sheets are determined by three-dimensional effects. Thus, a 2-D finite element analysis based on a local plasticity theory will give a physically unrealistic mesh dependent solution. This, in spite of the fact that the stress state, is essentially two-dimensional. By incorporating a length scale with relation to the thickness of the sheet, it is demonstrated how a 2-D finite element analysis based on a gradient dependent plasticity theory can give a good approximation of the post-necking behaviour. This is illustrated by numerical comparison of results from a full 3-D finite element analysis, with results from a 2-D finite element model based on a finite strain version of a gradient dependent J₂-flow theory. Some numerical problems in the modeling will be discussed briefly.     

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.     

Solution of some periodic elastoplastic problems

The article considers an elastoplastic problem for a plane weakened by an infinite number of round openings. It is assumed that the level of the stresses and the distance between the openings are such that the round openings are completely enveloped by the corresponding plastic zone; under these circumstances, the adjacent plastic regions do not coalesce. The article also considers the inverse elastoplastic problem under conditions of plane strain for an unbounded plane, weakened by a periodic series of openings. A number of communications have been devoted to periodic problems in the theory of elasticity and plasticity with an unknown boundary [1–8]. In distinction from [1–8], in which the method of perturbations was used, another method is used to solve periodic elastoplastic problems, making it possible to obtain a solution with any arbitrary relative dimensions of the region.     

Transient thermoelastoplastic stress analysis for a hollow sphere using the incremental theory of plasticity

An analysis based on the incremental strain theory is formulated for solving the problem of an elastoplastic hollow sphere subjected to a transient temperature distribution. Thermal and material properties are assumed to be temperature dependent and the behaviour of the medium to be characterized by the Ramberg-Osgood stress-strain relation. A method of successive elastic solutions is used to obtain a numerical solution. An illustrative example shows that the effective stress is not a monotonie function of the radius, but is much dependent on the history, gradient, and distribution of the temperature in the hollow sphere. In addition, unloading in the plastically deformed region is confirmed from the detailed discussion on the distribution of strains. As a result, the analysis based on the total strain theory is not permissible for solving this kind of elastoplastic problems subjected to transient thermal loading. In the following analysis the problem is treated in a quasi-static sense and the inertia terms in the thermoelastoplastic equations are neglected.     

Fracture of bicrystal metal/ceramic interfaces: A study via the mechanism-based strain gradient crystal plasticity theory

Two continuum mechanical models of crystal plasticity theory namely, conventional crystal plasticity theory and mechanism-based crystal plasticity theory, are used to perform a comparative study of stresses that are reached at and ahead of the crack tip of a bicrystal niobium/alumina specimen. Finite element analyses are done for a stationary crack tip and growing cracks using a cohesive modelling approach. Using mechanism-based strain gradient crystal plasticity theory the stresses reached ahead of the crack tip are found to be two times larger than the stresses obtained from conventional crystal plasticity theory. Results also show that strain gradient effects strongly depend on the intrinsic material length to the size of plastic zone ratio (l/R₀). It is found that the larger the (l/R₀) ratio, the higher the stresses reached using mechanism-based strain gradient crystal plasticity theory. An insight into the role of cohesive strength and work of adhesion in macroscopic fracture is also presented which can be used by experimentalists to design better bimaterials by varying cohesive strength and work of adhesion.     

A mathematical basis for strain-gradient plasticity theory—Part I: Scalar plastic multiplier

Strain-gradient plasticity theories are reviewed in which some measure of the plastic strain rate is treated as an independent kinematic variable. Dislocation arguments are invoked in order to provide a physical basis for the hardening at interfaces. A phenomenological, flow theory version of gradient plasticity is constructed in which stress measures, work-conjugate to plastic strain and its gradient, satisfy a yield condition. Plastic work is also done at internal interfaces and a yield surface is postulated for the work-conjugate stress quantities at the interface. Thereby, the theory has the potential to account for grain size effects in polycrystals. Both the bulk and interfacial stresses are taken to be dissipative in nature and due attention is paid to ensure that positive plastic work is done. It is shown that the mathematical structure of the elasto-plastic strain-gradient theory has similarities to conventional rigid-plasticity theory. Uniqueness and extremum principles are constructed for the solution of boundary value problems.     

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

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

AN INVESTIGATION ON ASYMPTOTIC NEAR-TIP FIELDS OF DYNAMIC CRACK IN AN ELASTIC-PERFECTLY PLASTIC COMPRESSIBLE MATERIAL

The stress and deformation fields near the tip of a mode-I dynamic crack steadilypropagating in an elastic-perfectly plastic compressible material are considered under plane strain condi-tions. Within the framework of infinitesimal displacement gradient theory, the material is character-ized by the Von Mises yield criterion and the associated J_2 flow theory of plasticity. Through rigorousmathematical analysis, this paper eliminates the possibilities of elastic unloading and continuousasymptotic fields with singular deformation, and then constructs a fully continuous and boundedasymptotic stress and strain field. It is found that in this solution there exists a parameter (?)_0 whichcannot be determined by asymptotic analysis but may characterize the effect of the far field. Lastly thevariations of continuous stresses, velocities and strains around the crack tip are given numerically fordifferent values of (?)_0.     

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.     

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.     

基于统一强度理论的岩石弹塑性本构模型及其数值实现

基于弹塑性力学理论,以统一强度准则为屈服准则,建立了考虑硬化/软化行为和应变率效应的岩石弹塑性本构模型;采用Fortran语言通过LS-DYNA的用户自定义材料接口(Umat)对该弹塑性本构模型进行编程,并把该程序生成求解器以达到对该模型进行应用的目的;通过岩石的单轴压缩实验和SHPB实验对所建的弹塑性本构模型进行验证,结果表明,该弹塑性本构模型能够反映岩石在准静态和动态下的力学行为。     

Two-way coupled analysis of lithium diffusion and diffusion induced finite elastoplastic bending of bilayer electrodes in lithium-ion batteries

A fully coupling model for the diffusion induced finite elastoplastic bending of bilayer electrodes in lithium-ion batteries is proposed. The effect of the mechanical stress on the lithium diffusion is accounted for by the mechanical part of the chemical potential derived from the Gibbs free energy along with the logarithmic stress and strain. Eight dimensionless parameters, governing the stress-assisted diffusion and the diffusion induced elastoplastic bending, are identified. It is found that the finite plasticity starting from the interface of the bilayer increases the chemical potential gradient and thereby facilitates the lithium diffusion. The full plastic flow makes the abnormal lithium concentration distribution possible, i.e., the concentration at the lithium inlet can be lower than the concentration at the interface (downstream). The increase in the thickness of the active layer during charging is much larger than the eigen-stretch due to lithiation, and this excess thickening is found to be caused by the lithiation induced plastic yield.     

A multiscale model for the ductile fracture of crystalline materials

In this paper, a multiscale model that combines both macroscopic and microscopic analyses is presented for describing the ductile fracture process of crystalline materials. In the macroscopic fracture analysis, the recently developed strain gradient plasticity theory is used to describe the fracture toughness, the shielding effects of plastic deformation on the crack growth, and the crack tip field through the use of an elastic core model. The crack tip field resulting from the macroscopic analysis using the strain gradient plasticity theory displayes the 1/2 singularity of stress within the strain gradient dominated region. In the microscopic fracture analysis, the discrete dislocation theory is used to describe the shielding effects of discrete dislocations on the crack growth. The result of the macroscopic analysis near the crack tip, i.e. a new K-field, is taken as the boundary condition for the microscopic fracture analysis. The equilibrium locations of the discrete dislocations around the crack and the shielding effects of the discrete dislocations on the crack growth at the microscale are calculated. The macroscopic fracture analysis and the microscopic fracture analysis are connected based on the elastic core model. Through a comparison of the shielding effects from plastic deformation and the discrete dislocations, the elastic core size is determined.