Strain gradient plasticity analysis of elasto-plastic contact between rough surfaces

From a microscopic point of view, the real contact area between two rough surfaces is the sum of the areas of contact between facing asperities. Since the real contact area is a fraction of the nominal contact area, the real contact pressure is much higher than the nominal contact pressure, which results in plastic deformation of asperities. As plasticity is size dependent at size scales below tens of micrometers, with the general trend of smaller being harder, macroscopic plasticity is not suitable to describe plastic deformation of small asperities and thus fails to capture the real contact area and pressure accurately. Here we adopt conventional mechanism-based strain gradient plasticity (CMSGP) to analyze the contact between a rigid platen and an elasto-plastic solid with a rough surface. Flattening of a single sinusoidal asperity is analyzed first to highlight the difference between CMSGP and J₂ isotropic plasticity. For the rough surface contact, besides CMSGP, pure elastic and J₂ isotropic plasticity analysis is also carried out for comparison. In all cases, the contact area A rises linearly with the applied load, but with a different slope which implies that the mean contact pressures are different. CMSGP produces qualitative changes in the distributions of local contact pressures compared with pure elastic and J₂ isotropic plasticity analysis, furthermore, bounded by the two.     

An analysis of nanoindentation in elasto-plastic solids

This paper examines the accuracy of the extracted elastic properties using the nanoindentation technique on elasto-plastic materials. The application of the correction factor evaluated in the linearly elastic case [Poon, B., Rittel, D., Ravichandran, G., 2008. An analysis of nanoindentation in linearly elastic solids. Int. J. Solids Structures 45 (24), 6018–6033.] on elastic–plastic materials is critically examined. It is then established that the accurate determination of the projected area of contact is found to be crucial for the accurate determination of elastic material properties. The conventional methods for the accurate determination of contact area are generally limited to ratios of Young’s modulus over yield stress, E/σ_y < 30 for elastic-perfectly plastic materials, which is too stringent for most materials. Thus, a new electrical resistance method is proposed to measure directly the projected contact area. Using numerical simulations, it was found that with the accurate determination of A, the error associated with the extracted elastic material properties is reduced by more than 50% in some cases. Using the newly proposed procedure, the error is also found to be independent of E/σ_y and the tip radius, ρ, and it is only a function of Poisson’s ratio, ν. This suggests that the errors might be due to the residual stresses at the plastic imprint that were found to depend on ν as well.     

Phenomenological study of parabolic and spherical indentation of elastic-ideally plastic material

A phenomenological study of parabolic and spherical indentation of elastic ideally plastic materials was carried out by using precise results of finite elements calculations. The study shows that no “pseudo-Hertzian” regime occurs during spherical indentation. As soon as the yield stress of the indented material is exceeded, a deviation from the, purely elastic Hertzian contact behaviour is found. Two elastic–plastic regimes and two plastic regimes are observed for materials of very large Young modulus to Yield stress ratio, E/σ_y. The first elastic–plastic regime corresponds to a strong evolution of the indented plastic zone. The first plastic regime corresponds to the commonly called “fully plastic regime”, in which the average indentation pressure is constant and equal to about three times the yield stress of the indented material. In this regime, the contact depth to penetration depth ratio tends toward a constant value, i.e. h_c/h = 1.47. h_c/h is only constant for very low values of yield strain (σ_y/E lower than 5 × 10^?6) when aE¹/Rσ_y is higher than 10,000. The second plastic regime corresponds to a decrease in the average indentation pressure and to a steeper increase in the pile-up. For materials with very large E/σ_y ratio, the second plastic regime appears when the value of the non-dimensional contact radius a/R is lower than 0.01. In the case of spherical and parabolic indentation, results show that the first plastic regime exists only for elastic-ideally plastic materials having an E/σ_y ratio higher than approximately 2.000.     

Size effects in generalised continuum crystal plasticity for two-phase laminates

The solutions of a boundary value problem are explored for various classes of generalised crystal plasticity models including Cosserat, strain gradient and micromorphic crystal plasticity. The considered microstructure consists of a two-phase laminate containing a purely elastic and an elasto-plastic phase undergoing single or double slip. The local distributions of plastic slip, lattice rotation and stresses are derived when the microstructure is subjected to simple shear. The arising size effects are characterised by the overall extra back stress component resulting from the action of higher order stresses, a characteristic length l_c describing the size-dependent domain of material response, and by the corresponding scaling law lⁿ as a function of microstructural length scale, l. Explicit relations for these quantities are derived and compared for the different models. The conditions at the interface between the elastic and elasto-plastic phases are shown to play a major role in the solution. A range of material parameters is shown to exist for which the Cosserat and micromorphic approaches exhibit the same behaviour. The models display in general significantly different asymptotic regimes for small microstructural length scales. Scaling power laws with the exponent continuously ranging from 0 to −2 are obtained depending on the values of the material parameters. The unusual exponent value −2 is obtained for the strain gradient plasticity model, denoted “curl H^p” in this work. These results provide guidelines for the identification of higher order material parameters of crystal plasticity models from experimental data, such as precipitate size effects in precipitate strengthened alloys.     

Experimental investigation of initial and subsequent yield surfaces for laminated metal matrix composites

The aim of this work is to construct yield surfaces to describe initial yielding and characterize hardening behavior of a highly anisotropic material. A methodology for constructing yield surfaces for isotropic materials using axial–torsion loading is extended to highly anisotropic materials. The technique uses a sensitive definition of yielding based on permanent strain rather than offset strain, and enables multiple yield points and multiple yield surfaces to be conducted on a single specimen. A target value of 20 × 10⁻⁶ is used for Al₂O₃ fiber reinforced aluminum laminates having a fiber volume fraction of 0.55. Sixteen radial probes are used to define the yield locus in the axial–shear stress plane. Initial yield surfaces for [0₄], [90₄], and [0/90]₂ fibrous aluminum laminates are well described by ellipses in the axial–shear stress plane having aspect ratios of 10, 2.5, and 3.3, respectively. For reference, the aspect ratio of the Mises ellipse for an isotropic material is 1.73. Initial yield surfaces do not have a tension–compression asymmetry. Four overload profiles (plus, ex, hourglass, and zee) are applied to characterize hardening of a [0/90]₂ laminate by constructing 30 subsequent yield surfaces. Parameters to describe the center and axes of an ellipse are regressed to the yield points. The results clearly indicate that kinematic hardening dominates so that material state evolution can be described by tracking the center of the yield locus. For a nonproportional overload of (σ, τ) = (500, 70) MPa, the center of the yield locus translated to (σ, τ) = (430, 37) MPa and the ellipse major axis was only 110 MPa.     

Yield criterion and elasto-plastic damage constitutive model for frozen sandy soil

A series of triaxial compression tests was carried out on a frozen sandy soil under confining pressures of 0–18 MPa at −6 °C. The experimental results indicate that, the strength of frozen sandy soil increases versus the increase in the confining pressures when σ₃ ? 3 MPa, but decreases when σ₃ > 3 MPa. This phenomenon is called the strengthening and weakening effects of confining pressures. A yield function, considering both effects, is proposed using the experimental method according to Drucker’s postulate, and the mathematical expression of the hardening parameter, which can describe the softening and hardening phenomenon, is provided. An elasto-plastic constitutive model for frozen sandy soil is developed. Based on the continuum damage theory, the cross anisotropic damage variables are deduced and their change regularities are investigated. Then the elasto-plastic damage constitutive model is proposed by introducing damage variables into elasto-plastic constitutive model. The validity of the model is verified by comparing its modeling results with experimental results obtained from triaxial tests. It is found that, this model can predict the deformation regularity of frozen soil exactly. It can simulate the stress–strain process under high confining pressures when pressure melting phenomena appear especially well.     

Matched asymptotic expansions in nonlinear fracture mechanics—III. In-plane loading of an elastic perfectly-plastic symmetric specimen

The matched expansion method, introduced by the authors in two earlier papers (1976) devoted to mode III loading, is applied to the practically important case of mode I loading of a symmetric specimen. The method allows the linear elastic far-field to be considered separately from the elasto-plastic near-tip field, except for coupling through a set of parameters that are determined explicitly in the matching. The effects of the plasticity are thus found, once and for all, from the solution of a set of standard elasto-plastic problems for a semi-infinite crack in an infinite body, whose properties may be tabulated. The solution for any particular specimen geometry and loading then follows from a small set of linear elastic solutions for the specimen, which define, through coefficients γ_ij appearing in their near-tip expansions, all the parameters in the “inner” and “outer” solutions. The effects of plasticity appear in these parameters only through a set of constants C_ti that define the far-field expansions of the “inner” (near-tip) solutions: they are material constants, depending upon the constitutive relation for the material, but not upon specimen geometry and loading. The J-integral, being obtainable from the far-field, is expressed as an explicit asymptotic series in the loading parameter ε, whose coefficients are given as functions of the “elastic” parameters γ_ij and the material constants C_i. It is demonstrated that a plastic-zone correction term, r_y, can be chosen to yield a two-term asymptotic expansion for J; the value of r_y depends upon the yielding model only through the constant C₁.The Dugdale (1960) model of yielding is treated, as a simple example for which all calculations can be performed analytically, and for which exact solutions are available for comparison.Finally, the near-tip solutions are constructed for a material obeying the Mises yield criterion and associated flow-rule, using a specially developed finite element program. The first eight of the constants C_i are tabulated, which suffice to define the J-integral up to terms of order ε⁶ (where ε is a loading parameter) and some representative near-tip features are displayed graphically. The computed value of C₁ shows that the conventionally adopted value for the plastic-zone correction r_y is too large by a factor of roughly 2.8, if it is to yield a genuine asymptotic estimate for J. As an example, the “elastic” parameters γ_ij are found, from a boundary collocation program, for a centre-cracked square plate subjected to tensile loading; and a plot of J versus load, and the plastic-zone shape at a particular load level, are displayed.     

Anisotropic hardening – numerical application of a cubic yield theory and consideration of variable r-values for sheet metal

The first part of this paper contains a polynomial yield condition of third order connected with evolution equations for material tensors of higher orders. They are formulated by formal generalisation of an approach by Danilov. The second part presents a possibility of taking into account the rotation of the yield surface as a result of a variable planar anisotropy (r-value) in sheet metal. This is done by an extension of the evolution equations, based on a quadratic yield function. The corresponding deformation law and the set of evolution equations are numerically integrated for selected loading paths in the subspaces σ₁,σ₂ and σ,τ. Some of the experimentally observed effects, such as the increasing curvature of the yield locus curve in the loading direction or the specific rotation of the yield surface, are correctly reproduced.     

A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming

In this paper an anisotropic material model based on non-associated flow rule and mixed isotropic–kinematic hardening was developed and implemented into a user-defined material (UMAT) subroutine for the commercial finite element code ABAQUS. Both yield function and plastic potential were defined in the form of Hill’s [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. R. Soc. Lond. A 193, 281–297] quadratic anisotropic function, where the coefficients for the yield function were determined from the yield stresses in different material orientations, and those of the plastic potential were determined from the r-values in different directions. Isotropic hardening follows a nonlinear behavior, generally in the power law form for most grades of steel and the exponential law form for aluminum alloys. Also, a kinematic hardening law was implemented to account for cyclic loading effects. The evolution of the backstress tensor was modeled based on the nonlinear kinematic hardening theory (Armstrong–Frederick formulation). Computational plasticity equations were then formulated by using a return-mapping algorithm to integrate the stress over each time increment. Either explicit or implicit time integration schemes can be used for this model. Finally, the implemented material model was utilized to simulate two sheet metal forming processes: the cup drawing of AA2090-T3, and the springback of the channel drawing of two sheet materials (DP600 and AA6022-T43). Experimental cyclic shear tests were carried out in order to determine the cyclic stress–strain behavior and the Bauschinger ratio. The in-plane anisotropy (r-value and yield stress directionalities) of these sheet materials was also compared with the results of numerical simulations using the non-associated model. These results showed that this non-associated, mixed hardening model significantly improves the prediction of earing in the cup drawing process and the prediction of springback in the sidewall of drawn channel sections, even when a simple quadratic constitutive model is used.     

Multiaxial yield behaviour of Al replicated foam

Multiaxial experiments are performed on replicated aluminium foam using a custom-built apparatus. The foam structure is isotropic, and features open monomodal pores in average diameter. Plane stress (σ₁, σ₂, σ₃=0) and axisymmetric (σ₁, σ₂=σ₃) yield envelopes are measured using cubical specimens, supplemented by tests on hollow cylindrical and uniaxial samples. In addition to the three stress components at 0.2% offset strain, the computer-controlled testing apparatus also measures the three instantaneous displacement vectors.Results show that the shape of the yield surface is independent of the relative density of the foam in the explored range (13-28%). Strain increment vectors lie, within error, roughly normal to the line traced through data points in stress space. Replicated foams feature asymmetric yield behaviour between tension and compression. The data additionally show an influence on the yield surface of the third stress tensor invariant (i.e., of the Lode angle). Simple general expressions for the yield surface are fitted to the data, leading to conclude that their behaviour is slightly better captured by parabolic rather than elliptic expressions dependent on all three stress invariants.     

On relations for general two-dimensional ideal plasticity problems under the full plasticity condition

Relations for two-dimensional ideal plasticity problems under the full plasticity condition are determined with material anisotropy, inhomogeneity, and compressibility properties taken into account. These properties are determined by the direction cosines of the principal stress, the coordinates of points in space, and the mean stress.For the yield strength we take a function of the form k = k(σ, n ₁, n ₂, n ₃, x, y, z). The desired relations are determined for the general plane ideal plasticity problem. The relations thus obtained are generalized to the cases of axisymmetric and spherical plasticity problems.     

Endochronic analysis for rate-dependentelasto-plastic deformation

This work presents a novel formulation of the scaling function of the intrinsic timescale. The new formulated scaling function and the rate-sensitivity function are used in theendochronic theory to simulate the rate-dependent elasto-plastic behaviors of materials. Severalcases of rate-dependent elasto-plastic material responses are discussed separately. Experimentaldata found in literature of 304 stainless steels and Ti₇Al₂Cb₁Tatitanium alloy for rate-dependent elasto-plastic responses are used for comparison. It is shownthat the theory adequately simulates the experimental results.     

Power-law creep and residual stresses in a carbopol gel

We report on the interplay between creep and residual stresses in a carbopol microgel. When a constant shear stress σ is applied below the yield stress σ _y, the strain is shown to increase as a power law of time, γ(t) = γ ₀ + (t/τ) ^α , with an exponent α = 0.39 ± 0.04 that is strongly reminiscent of Andrade creep in hard solids. For applied shear stresses lower than some typical value σ _c ? 0.2σ _y, the microgel experiences a more complex, anomalous creep behavior, characterized by an initial decrease of the strain, that we attribute to the existence of residual stresses of the order of σ _c that persist after a rest time under a zero shear rate following preshear. The influence of gel concentration on creep and residual stresses are investigated as well as possible aging effects. We discuss our results in light of previous works on colloidal glasses and other soft glassy systems.     

On the characteristic equations of elasto-plastic flow

Three dimensional characteristic surfaces (slip surfaces) of elasto-plastic Navier's equations and the criteria for their existence are discussed, and the solutions are also applied to two dimensional cases. By making use of isotropic yield function, the following results are proved. If, and only if the plastic/elastic moduli ratio is zero and Det$\text{(φ}{}_{\mathrm{ij}}\text{)=0,(φ}{}_{\mathrm{ij}}\text{=}\text{?φ}\text{?σ}{}_{\mathrm{ij}}$, φ: yield function,σ_ij: stress tensor), characteristic surfaces exist. There are two and only two characteristic surface elements at each point, and they are identical with the surfaces of maximum shearing stress.     

State vector approach to analysis of multilayered magneto-electro-elastic plates

The state vector equations for three dimensional, orthotropic and linearly magneto-electro-elastic media are derived from the governing equations by eliminating σ_x, σ_y·τ_xy, B_x, B_y, D_x and D_y. An efficient method is presented for analysis of multilayered magneto-electro-elastic plates. The methodology is based on the mixed formulation, in which basic unknowns are formed by collecting not only displacements, electrical potential and magnetic potential but also some of stresses, electrical displacements, and magnetic induction. As special case, simply supported and multilayered rectangular plate is analyzed under the surface loading. Numerical results are presented graphically. The procedure of numerical calculation shows that the formulation presented here is simple and direct.     

Yielding and intrinsic plasticity of Ti–Zr–Ni–Cu–Be bulk metallic glass

Bulk metallic glass with composition Ti₄₀Zr₂₅Ni₈Cu₉Be₁₈ exhibits considerably high compressive yield stress, significant plasticity (with a concomitant vein-like fracture morphology) and relatively low density. Yielding and intrinsic plasticity of this alloy are discussed in terms of its thermal and elastic properties. An influence of normal stresses acting on the shear plane is evidenced by: (i) the fracture angle (<45°) and (ii) finite-element simulations of nanoindentation curves, which require the use of a specific yield criterion, sensitive to local normal stresses acting on the shear plane, to properly match the experimental data. The ratio between hardness and compressive yield strength (constraint factor) is analyzed in terms of several models and is best adjusted using a modified expanding cavity model incorporating a pressure-sensitivity index defined by the Drucker–Prager yield criterion. Furthermore, comparative results from compression tests and nanoindentation reveal that deformation also causes strain softening, a phenomenon which is accompanied with the occurrence of serrated plastic flow and results in a so-called indentation size effect (ISE). A new approach to model the ISE of this metallic glass using the free volume concept is presented.     

The coefficient of normal restitution for the Hertzian contact of two rough spheres colliding in a viscous fluid

We use previous theoretical results for the added mass, history and lubrication forces between two spheres colliding in a fluid with viscosity ν to investigate the effect of viscous dissipation on the coefficient of restitution during contact. We assume that the mechanical interaction is governed by Hertzian mechanical contact of small duration τ and that the minimum approach distance between particles is approximately equal to the height σ of surface micro-asperities. A non-dimensionalization of the equation of motion indicates that the contact dynamics is governed by two parameters – the ratio ϵ of the surface roughness σ and the sphere radius a, and a contact Stokes number defined as St_c = σ²/ντ. An asymptotic solution of the equation of motion in the limit of small ϵ/St_c is used to obtain an explicit expression for the coefficient of restitution during contact and the latter is compared with estimates based on numerical solutions of the non-linear equation of motion.     

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.     

Large strain deformation of polycrystalline metals at low homologous temperatures

Large strain compression data (true strains to about ?3.0) are presented for polycrystalline α U and α Fe at room temperature. The results, together with other published data at low homologous temperatures (≈0.2 T_m), where T_m is the absolute melting temperature, suggest that a steady-state flow stress σ_s is approached after extensive strain-hardening, α U exhibits a very high strain-hardening rate, with σ_s ≈ 2900 MPa (420 ksi) indicating that cold-working is a very potent method of strengthening this metal. All the data evaluated can be fit by the stress-strain relation σ = σ_s? exp (?(Nε)^p)(σ_s? σ_y), where σ_y is the yield stess, p is a constant equal to a for the metals analyzed, N is a constant associated with the strain-hardening characteristics of a material, σ is true stress, and ε is true strain.     

Optimal layout of trusses with finite numbers of joints2

The paper concerns the design of a truss with a given finite number of joints to carry a given system of loads in such a manner that the axial stress in any bar remains within the given allowable range [?σ₀, cσ₀], σ₀ > 0 being a reference tensile stress and c > 0 an arbitrary factor. The efficiency of a truss of this kind is defined as the quotient of the minimum amount of material required for this type of design by the amount of material needed for the considered truss. A method is described for determining a layout whose efficiency has at least a given value.