摩擦因素对微压痕实验的影响

材料在发生微米亚微米尺度不均匀变形时表现出尺度效应,微压痕问题是典型代表。目前的微压痕研究中普遍采用无摩擦假设,而相对缺少对摩擦因素的深入研究。基于CMSG低阶应变梯度塑性理论,通过有限元分析研究了摩擦因素对微压痕实验的影响。结果表明对于目前普遍采用的较钝的Berkovich压头(简化为圆锥后顶角为140.6?),摩擦的影响并不重要,可以简化为无摩擦情况。     

Temperature dependence of material length scale for strain gradient plasticity and its effect on near‐tip opening displacement

This paper investigates the temperature dependence of the material length scale in the conventional mechanism‐based strain gradient (CMSG) plasticity theory. The work reported here also examines the plastic strain gradient effect on the opening displacement near a sharp crack tip. The study examines the mechanical properties of two typical structural steels (S355 and S690) in onshore and offshore structures at two different temperatures (20 and 300 °C) through both the uniaxial tension test and the indentation test. The CMSG‐based finite element analysis then confirms a constant material length scale for these two steels at the two tested temperatures, despite the apparent temperature dependence of the macroscopic material parameters measured from the tension test. Using the calibrated material length scale, the subsequent numerical study demonstrates that the magnitude of the near‐tip crack opening displacement computed by the CMSG theory remains significantly lower than that computed from the classical plasticity.     

Non‐uniform plastic deformation of micron scale objects

Significant increases in apparent flow strength are observed when non‐uniform plastic deformation of metals occurs at the scale ranging from roughly one to ten microns. Several basic plane strain problems are analysed numerically in this paper based on a new formulation of strain gradient plasticity. The problems are the tangential and normal loading of a finite rectangular block of material bonded to rigid platens and having traction‐free ends, and the normal loading of a half‐space by a flat, rigid punch. The solutions illustrate fundamental features of plasticity at the micron scale that are not captured by conventional plasticity theory. These include the role of material length parameters in establishing the size dependence of strength and the elevation of resistance to plastic flow resulting from constraint on plastic flow at boundaries. Details of the finite element method employed in the numerical analysis of the higher order gradient theory will be discussed and related to prior formulations having some of the same features. Copyright © 2003 John Wiley & Sons, Ltd.     

Interface crack problems with strain gradient effects

In this paper, the strain gradient theory proposed by Chen and Wang (2001a, 2002b) is used to analyze an interface crack tip field at micron scales. Numerical results show that at a distance much larger than the dislocation spacing the classical continuum plasticity is applicable; but the stress level with the strain gradient effect is significantly higher than that in classical plasticity immediately ahead of the crack tip. The singularity of stresses in the strain gradient theory is higher than that in HRR field and it slightly exceeds or equals to the square root singularity and has no relation with the material hardening exponents. Several kinds of interface crack fields are calculated and compared. The interface crack tip field between an elastic-plastic material and a rigid substrate is different from that between two elastic-plastic solids. This study provides explanations for the crack growth in materials by decohesion at the atomic scale.     

Mode III crack growth in linear hardening materials with strain gradient effects

The flow-theory version of couple stress strain gradient plasticity is adopted for investigating the asymptotic fields near a steadily propagating crack-tip, under Mode III loading conditions. By incorporating a material characteristic length, typically of the order of few microns for ductile metals, the adopted constitutive model accounts for the microstructure of the material and can capture the strong size effects arising at small scales. The effects of microstructure result in a substantial increase in the singularities of the skew-symmetric stress and couple stress fields, which occurs also for a small hardening coefficient. The symmetric stress field turns out to be non-singular according to the asymptotic solution for the stationary crack problem in linear elastic couple stress materials. The performed asymptotic analysis can provide useful predictions about the increase of the traction level ahead of the crack-tip due to the sole contribution of the rotation gradient, which has been found relevant and non-negligible at the micron scale.     

Constraint effects on the elastic–plastic fracture behaviour in strain gradient solids

ABSTRACT Crack‐tip constraint effects (or T‐stress effects) on the elastic–plastic fracture behaviour in strain gradient materials are analysed in the present study. The T‐stress effects on the stress distributions along the plane ahead of the stationary and growing crack tip, respectively, are analysed by using the Fleck and Hutchinson strain gradient plasticity formation. For a steadily growing crack, the T‐stress effects on the steady‐state fracture toughness are analysed by adopting the embedded fracture process zone model. In addition, the analysis for the growing crack is applied to an interfacial cracking experiment for a metal/ceramic system, and the material length‐scale parameter appearing in the strain gradient plasticity theory is predicted. In the present analyses, a new finite element method specially designed for strain gradient problems by Wei and Hutchinson is adopted.     

A micro‐mechanical damage model based on gradient plasticity: algorithms and applications

As soon as material failure dominates a deformation process, the material increasingly displays strain softening and the finite element computation is significantly affected by the element size. Without remedying this effect in the constitutive model one cannot hope for a reliable prediction of the ductile material failure process. In the present paper, a micro‐mechanical damage model coupled to gradient‐dependent plasticity theory is presented and its finite element algorithm is discussed. By incorporating the Laplacian of plastic strain into the damage constitutive relationship, the known mesh‐dependence is overcome and computational results are uniquely correlated with the given material parameters. The implicit C¹ shape function is used and can be transformed to arbitrary quadrilateral elements. The introduced intrinsic material length parameter is able to predict size effects in material failure. Copyright © 2002 John Wiley & Sons, Ltd.     

Analytic and numerical studies on mode I and mode II fracture in elastic-plastic materials with strain gradient effects

This paper presents a study on fracture of materials at microscale (∼1 μm) by the strain gradient theory (Fleck and Hutchinson, 1993; Fleck et al., 1994). For remotely imposed classical K fields, the full-field solutions are obtained analytically or numerically for elastic and elastic-plastic materials with strain gradient effects. The analytical elastic full-field solution shows that stresses ahead of a crack tip are significantly higher than their counterparts in the classical K fields. The sizes of dominance zones for mode I and mode II near-tip asymptotic fields are 0.3l and 0.5l,while strain gradient effects are observed within land 2l to the crack tip, respectively, where l is the intrinsic material length in strain gradient theory and is on the order of microns in strain gradient plasticity (Fleck et al., 1994; Nix and Gao, 1998; Stolken and Evans, 1997). The Dugdale–Barenblatt type plasticity model is obtained to provide an estimation of plastic zone size for mode II fracture in materials with strain grain effects. The finite element method is used to investigate the small-scale-yielding solution for an elastic-power law hardening solid. It is found that the size of the dominance zone for the near-tip asymptotic field is the intrinsic material lengthl. For mode II fracture under the small-scale-yielding condition, transition from the remote classical K _IIfield to the near-tip asymptotic field in strain gradient plasticity goes through the HRR field only when K _IIis relatively large such that the plastic zone size is much larger than the intrinsic material length l. For mode I fracture under small-scale-yielding condition, however, transition from the remote classical K _I field to the near-tip asymptotic field in strain gradient plasticity does not go through the HRR field, but via a plastic zone. This revised version was published online in July 2006 with corrections to the Cover Date.     

Capturing strain localization in dense sands with random density

This paper presents a three‐invariant constitutive framework suitable for the numerical analyses of localization instabilities in granular materials exhibiting unstructured random density. A recently proposed elastoplastic model for sands based on critical state plasticity is enhanced with the third stress invariant to capture the difference in the compressive and extensional yield strengths commonly observed in geomaterials undergoing plastic deformation. The new three‐invariant constitutive model, similar to its two‐invariant predecessor, is capable of accounting for meso‐scale inhomogeneities as well as material and geometric nonlinearities. Details regarding the numerical implementation of the model into a fully nonlinear finite element framework are presented and a closed‐form expression for the consistent tangent operator, whose spectral form is used in the strain localization analyses, is derived. An algorithm based on the spectral form of the so‐called acoustic tensor is proposed to search for the necessary conditions for deformation bands to develop. The aforementioned framework is utilized in a series of boundary‐value problems on dense sand specimens whose density fields are modelled as exponentially distributed unstructured random fields to account for the effect of inhomogeneities at the meso‐scale and the intrinsic uncertainty associated with them. Copyright © 2006 John Wiley & Sons, Ltd.     

Variational multiscale methods to embed the macromechanical continuum formulation with fine‐scale strain gradient theories

A variational basis is presented to link fine‐scale theories of material behaviour with the classical, macromechanical continuum theory. The approach is based on the weak form of the linear momentum balance equations, and a separation of the weighting function and displacement fields into coarse and fine‐scale components. Coarse and fine‐scale weak forms are defined. The latter is used to introduce a strain gradient theory that operates at finer scales of deformation. Attention is focused upon applications requiring the enhanced physical accuracy of the fine‐scale strain gradient theory, without the computational cost of discretization that spans the range from coarse to fine scales. A variationally consistent method is developed to embed the fine‐scale strain gradient theory in the macromechanical formulation. The embedding is achieved by eliminating the fine‐scale displacement field from the problem. Two examples demonstrate the numerical efficiency of the method, while retaining physical and mathematical properties of the fine‐scale strain gradient theory. Copyright © 2003 John Wiley & Sons, Ltd.     

Regularization of material instabilities by meshfree approximations with intrinsic length scales

Meshfree approximation, such as Moving Least Square (MLS) and Reproducing Kernel (RK) approximations, possess intrinsic non‐local properties. These non‐local properties of meshfree approximations are exploited to incorporate an intrinsic length scale which regularizes problems with material instabilities. The discrete equilibrium equation is obtained by employing an assumed strain method in the Galerkin approximation. This proposed method is essentially uniformly non‐local, but in contrast to non‐local finite elements, no kinematic modes are observed. Gradient‐type regularization can also be modelled by this method without the additional boundary conditions and other complications of the conventional gradient methods. Numerical examples show that the displacement‐based MLS/RK formulation (1‐level regularization) is sufficient to remedy mesh‐sensitivity in damage‐induced strain localization. For strain localization associated with plasticity, a two‐level MLS/RK regularization in displacement and strain shown to be effective. Copyright © 2000 John Wiley & Sons, Ltd.     

Micromechanical modeling of intrinsic and specimen size effects in microforming

Size effect is a crucial phenomenon in the microforming processes of metallic alloys involving only limited amount of grains. At this scale intrinsic size effect arises due to the size of the grains and the specimen/statistical size effect occurs due to the number of grains where the properties of individual grains become decisive on the mechanical behavior of the material. This paper deals with the micromechanical modeling of the size dependent plastic response of polycrystalline metallic materials at micron scale through a strain gradient crystal plasticity framework. The model is implemented into a Finite Element software as a coupled implicit user element subroutine where the plastic slip and displacement fields are taken as global variables. Uniaxial tensile tests are conducted for microstructures having different number of grains with random orientations in plane strain setting. The influence of the grain size and number on both local and macroscopic behavior of the material is investigated. The attention is focussed on the effect of the grain boundary conditions, deformation rate and the grain size on the mechanical behavior of micron sized specimens. The model is intrinsically capable of capturing both experimentally observed phenomena thanks to the incorporated internal length scale and the crystallographic orientation definition of each grain.     

A finite element displacement formulation for gradient elastoplasticity

We present a second gradient elastoplastic model for strain‐softening materials based entirely on a finite element displacement formulation. The stress increment is related to both the strain increment and its Laplacian. The displacement field is the only field needed to be discretized using a C¹ continuity element. The required higher‐order boundary conditions arise naturally from the displacement field. The model is developed to regularize the ill‐posedness caused by strain‐softening material behaviour. The gradient terms in the constitutive equations introduce an extra material parameter with dimensions of length allowing robust modelling of the post‐peak material behaviour leading to localization of deformation. Mesh insensitivity is demonstrated by modelling localization of deformation in biaxial tests. It is shown that both the thickness and inclination of the shear‐band zone are insensitive to the mesh directionality and refinement and agree with the expected theoretical and experimental values. Copyright © 2001 John Wiley & Sons, Ltd.     

Multi‐scale methods

In this paper four multiple scale methods are proposed. The meshless hierarchical partition of unity is used as a multiple scale basis. The multiple scale analysis with the introduction of a dilation parameter to perform multiresolution analysis is discussed. The multiple field based on a 1‐D gradient plasticity theory with material length scale is also proposed to remove the mesh dependency difficulty in softening/localization problems. A non‐local (smoothing) particle integration procedure with its multiple scale analysis are then developed. These techniques are described in the context of the reproducing kernel particle method. Results are presented for elastic‐plastic one‐dimensional problems and 2‐D large deformation strain localization problems to illustrate the effectiveness of these methods. Copyright © 2000 John Wiley & Sons, Ltd.     

Determination of nanoindentation size effects and variable material intrinsic length scale for body-centered cubic metals

It is well-known by now that the micro and nanoindentation hardness of metallic materials displays a strong size effect. The objective of this work is to formulate a micromechanical-based model for Temperature and Rate Indentation Size Effects (TRISE) for body centered cubic (BCC) metals encountered in nanoindentation experiments. In this regard, two physically based models are proposed here in order to capture the TRISE in single and polycrystalline materials by considering different expressions of the geometrical necessary dislocations (GNDs) density.The gradient plasticity theory formulates a constitutive framework on the continuum level that bridges the gap between the micromechanical plasticity and the classical continuum plasticity by incorporating the material length scale. A micromechanical-based model of variable material intrinsic length scale is also developed in the present work. The proposed length scale allows for variations in temperature and strain rate and its dependence on the grain size and accumulated plastic strain.The results of indentation experiments performed on niobium, tungsten, and single- and polycrystalline commercially pure iron (very similar to iron alloys) are used here to implement the aforementioned framework in order to predict simultaneously the TRISE and variable length scale at different temperatures, strain rates and various distances from the grain boundary. Numerical analysis is performed using the ABAQUS/VUMAT software with a physically based viscoplastic constitutive model.     

Experimental study on tensile property of AZ31B magnesium alloy at different high strain rates and temperatures

As the lightest metal material, magnesium alloy is widely used in the automobile and aviation industries. Due to the crashing of the automobile is a process of complicated and highly nonlinear deformation. The material deformation behavior has changed significantly compared with quasi-static, so the deformation characteristic of magnesium alloy material under the high strain rate has great significance in the automobile industry. In this paper, the tensile deformation behavior of AZ31B magnesium alloy is studied over a large range of the strain rates, from 700 s⁻¹ to 3 × 10³ s⁻¹ and at different temperatures from 20 to 250 °C through a Split-Hopkinson Tensile Bar (SHTB) with heating equipment. Compared with the quasi-static tension, the tensile strength and fracture elongation under high strain rates is larger at room temperature, but when at the high strain rates, fracture elongation reduces with the increasing of the strain rate at room temperature, the adiabatic temperature rising can enhance the material plasticity. The morphology of fracture surfaces over wide range of strain rates and temperatures are observed by Scanning Electron Microscopy (SEM). The fracture appearance analysis indicates that the fracture pattern of AZ31B in the quasi-static tensile tests at room temperature is mainly quasi-cleavage pattern. However, the fracture morphology of AZ31B under high strain rates and high temperatures is mainly composed of the dimple pattern, which indicates ductile fracture pattern. The fracture mode is a transition from quasi-cleavage fracture to ductile fracture with the increasing of temperature, the reason for this phenomenon might be the softening effect under the high strain rates.     

Assumed‐deformation gradient finite elements with nodal integration for nearly incompressible large deformation analysis

An assumed‐strain finite element technique for non‐linear finite deformation is presented. The weighted‐residual method enforces weakly the balance equation with the natural boundary condition and also the kinematic equation that links the elementwise and the assumed‐deformation gradient. Assumed gradient operators are derived via nodal integration from the kinematic‐weighted residual. A variety of finite element shapes fits the derived framework: four‐node tetrahedra, eight‐, 27‐, and 64‐node hexahedra are presented here. Since the assumed‐deformation gradients are expressed entirely in terms of the nodal displacements, the degrees of freedom are only the primitive variables (displacements at the nodes). The formulation allows for general anisotropic materials and no volumetric/deviatoric split is required. The consistent tangent operator is inexpensive and symmetric. Furthermore, the material update and the tangent moduli computation are carried out exactly as for classical displacement‐based models; the only deviation is the consistent use of the assumed‐deformation gradient in place of the displacement‐derived deformation gradient. Examples illustrate the performance with respect to the ability of the present technique to resist volumetric locking. A constraint count can partially explain the insensitivity of the resulting finite element models to locking in the incompressible limit. Copyright © 2008 John Wiley & Sons, Ltd.     

Simulation of effects of particle size and volume fraction on Al alloy strength,elongation, and toughness by using strain gradient plasticity concept

By using the constitutive equation based on the mechanism-based strain gradient plasticity with finite element software, the yield strength, uniform elongation, and toughness of aluminum alloy 6063 with different grain sizes, different particle diameters and volume fractions were studied numerically. The toughness is defined as the product of yield strength and uniform elongation. The calculation results indicate that the grain refinement and particle refinement cannot substantially improve the uniform elongation but can increase the yield strength of Al alloy when the grain size is on the order of the micron and submicron scale. When the grain size less than 2 μm, Al alloys usually exhibit high strength and low uniform elongation, and when the grain size greater than 5 μm, the materials exhibit low strength and high elongation; in either case the toughness is low. However, in the grain size of several micrometers, the toughness of Al alloy is the highest.     

An arbitrary Lagrangian–Eulerian finite element method for finite strain plasticity

This paper presents a new arbitrary Lagrangian–Eulerian (ALE) finite element formulation for finite strain plasticity in non‐linear solid mechanics. We consider the models of finite strain plasticity defined by the multiplicative decomposition of the deformation gradient in an elastic and a plastic part ( F = F ^e F ^p), with the stresses given by a hyperelastic relation. In contrast with more classical ALE approaches based on plastic models of the hypoelastic type, the ALE formulation presented herein considers the direct interpolation of the motion of the material with respect to the reference mesh together with the motion of the spatial mesh with respect to this same reference mesh. This aspect is shown to be crucial for a simple treatment of the advection of the plastic internal variables and dynamic variables. In fact, this advection is carried out exactly through a particle tracking in the reference mesh, a calculation that can be accomplished very efficiently with the use of the connectivity graph of the fixed reference mesh. A staggered scheme defined by three steps (the smoothing, the advection and the Lagrangian steps) leads to an efficient method for the solution of the resulting equations. We present several representative numerical simulations that illustrate the performance of the newly proposed methods. Both quasi‐static and dynamic conditions are considered in these model examples. Copyright © 2003 John Wiley & Sons, Ltd.     

Two‐point constitutive equations and integration algorithms for isotropic‐hardening rate‐independent elastoplastic materials in large deformation

This paper presents alternative forms of hyperelastic–plastic constitutive equations and their integration algorithms for isotropic‐hardening materials at large strain, which are established in two‐point tensor field, namely between the first Piola–Kirchhoff stress tensor and deformation gradient. The eigenvalue problems for symmetric and non‐symmetric tensors are applied to kinematics of multiplicative plasticity, which imply the transformation relationships of eigenvectors in current, intermediate and initial configurations. Based on the principle of plastic maximum dissipation, the two‐point hyperelastic stress–strain relationships and the evolution equations are achieved, in which it is considered that the plastic spin vanishes for isotropic plasticity. On the computational side, the exponential algorithm is used to integrate the plastic evolution equation. The return‐mapping procedure in principal axes, with respect to logarithmic elastic strain, possesses the same structure as infinitesimal deformation theory. Then, the theory of derivatives of non‐symmetric tensor functions is applied to derive the two‐point closed‐form consistent tangent modulus, which is useful for Newton's iterative solution of boundary value problem. Finally, the numerical simulation illustrates the application of the proposed formulations. Copyright © 2008 John Wiley & Sons, Ltd.