Shakedown theorems for elastic–plastic solids in the framework of gradient plasticity

Static and kinematic shakedown theorems are given for a class of generalized standard materials endowed with a hardening saturation surface in the framework of strain gradient plasticity. The so-called residual-based gradient plasticity theory is employed. The hardening law admits a hardening potential, which is a C¹-continuous function of a set of kinematic internal variables and of their spatial gradients, and is required to satisfy a global sign restriction (but not to be necessarily convex). The totally produced, the accumulated and the freely moving dislocations per unit volume, distinguished as statistically stored and geometrically necessary ones, are in this way accounted for. Like for a generalized standard material, the shakedown safety factor is found to depend on the (generally size dependent) yield and saturation limits, but not on the particular differential-type hardening law of the material.     

Effect of large elastic strains on cavitation instability predictions for elastic–plastic solids

For an infinite solid containing a void, the cavitation instability limit is defined as the remote stress–and strain state, at which the void grows without bound, driven by the elastic energy stored in the surrounding material. Such cavitation limits have been analysed by a number of authors for metal plasticity as well as for nonlinear elastic solids. The analyses for elastic–plastic solids are here extended to consider the effect of a large initial yield strain, and it is shown how the critical stress value decays for increasing value of the yield strain. Analyses are carried out for remote hydrostatic tension as well as for more general axisymmetric remote stress field, with an initially spherical void. Different levels of strain hardening are considered.     

A unified expression of elastic–plastic notch stress–strain calculation in bodies subjected to multiaxial cyclic loading

In this paper, a simplified thermodynamics analysis of cyclic plastic deformation is performed in order to establish an energy transition relation for describing the elastic–plastic stress and strain behavior of the notch-tip material element in bodies subjected to multiaxial cyclic loads. Based on the thermodynamics analysis, it is deduced that in the case of elastic–plastic deformation, Neuber’s rule inevitably overestimates the actual stress and strain at the notch tip, while the equivalent strain energy density (ESED) method tends to underestimate the actual notch-tip stress and strain. According to the actual energy conversion occurring in the notch-tip material element during cyclic plastic deformation, a unified expression for estimating the elastic–plastic notch stress–strain responses in bodies subjected to multiaxial cyclic loads is developed, of which Neuber’s rule and the ESED method become two particular cases, i.e. upper and lower bound limits of the notch stress and strain estimations. This expression is verified experimentally under both proportional and non-proportional multiaxial cyclic loads and a good agreement between the calculated and the measured notch strains has been achieved. It is also shown that in the case of multiaxial cyclic loading, the unified expression distinctly improves the accuracy of the notch-tip stress–strain estimations in comparison with Neuber’s rule and the ESED method. The unified expression of the notch stress–strain calculation developed in this paper can thus provide a more logical approximate approach for estimating the elastic–plastic notch-tip stress and strain responses of components subjected to lengthy multiaxial cyclic loading histories for local strain approach-based fatigue-crack-initiation life prediction.     

A new analytical model for the elastic–plastic behaviour of spot welded joints subjected to orthogonal load

An analytical procedure for the evaluation of the elastic–plastic stiffness behaviour of spot welded joints is presented. The procedure is based on a new model of spot weld region: a circular plate having variable thickness with a central rigid nugget, which is resolved using an original analytical method.The closed-form solution allows to describe the displacement of a rigid nugget when an axial orthogonal load is applied on the plate while plasticity and large deflections are present. The goal is to reach a reliable spot weld region model which can be used as the basis to develop a spot weld element in FE analysis even when plasticity and large deflections are in effect.The procedure is as completely original as no other can be found in the technical literature, and it has been applied to some examples of plates usually employed for spot weld analysis. The analytical results obtained by using the new general relations precisely match those obtained modelling spot weld area by FEA.     

Determination of the effective elastic–plastic response of metal–ceramic composites

The mechanical response of metal–ceramic composites is analysed through a homogenization model accounting for the mechanical behaviour of the constituent materials. In order to achieve this purpose a nonlinear homogenization method based on the phase field approach has been suitably implemented into a numerical code. A prescribed homogenized strain state is applied to a unit volume element of a metal–ceramic composite with proportional loading in which all components of the strain tensor are proportional to one scalar parameter. The mechanical response of the material has been modeled by considering a von Mises plasticity model for the metal phase and a Drucker–Prager associative elastic–plastic material model for the ceramic phase. A two stages plasticity has been obtained in which inelastic strain develops in the metal phase followed by a fully plastic response. A comparison with a finite element model of the stress–strain response of an axisymmetric unit cell has been carried out with the purpose to validate the homogenization based modeling presented in the paper. Plastic parameters of a Drucker–Prager yield surface for the homogenized composite have been calculated at different materials compositions. Associative Drucker–Prager plasticity has been found to be accurate for high ceramic content.     

Numerical analysis and elastic–plastic deformation behavior of anisotropically damaged solids

The present paper is concerned with the numerical modelling of the large elastic–plastic deformation behavior and localization prediction of ductile metals which are sensitive to hydrostatic stress and anisotropically damaged. The model is based on a generalized macroscopic theory within the framework of nonlinear continuum damage mechanics. The formulation relies on a multiplicative decomposition of the metric transformation tensor into elastic and damaged-plastic parts. Furthermore, undamaged configurations are introduced which are related to the damaged configurations via associated metric transformations which allow for the interpretation as damage tensors. Strain rates are shown to be additively decomposed into elastic, plastic and damage strain rate tensors. Moreover, based on the standard dissipative material approach the constitutive framework is completed by different stress tensors, a yield criterion and a separate damage condition as well as corresponding potential functions. The evolution laws for plastic and damage strain rates are discussed in some detail. Estimates of the stress and strain histories are obtained via an explicit integration procedure which employs an inelastic (damage-plastic) predictor followed by an elastic corrector step. Numerical simulations of the elastic–plastic deformation behavior of damaged solids demonstrate the efficiency of the formulation. A variety of large strain elastic–plastic-damage problems including severe localization is presented, and the influence of different model parameters on the deformation and localization prediction of ductile metals is discussed.     

A new approach for computing the steady state fluid–structure interaction response of periodic problems

A special type of fluid–structure interaction (FSI) problems are problems with periodic boundary conditions like in turbomachinery. The steady state FSI response of these problems is usually calculated with similar techniques as used for transient FSI analyses. This means that, when the fluid and structure problem are not simultaneously solved with a monolithic approach, the problem is partitioned into a fluid and structural part and that each time step coupling iterations are performed to account for strong interactions between the two sub-domains. This paper shows that a time-partitioned FSI computation can be very inefficient to compute the steady state FSI response of periodic problems. A new approach is introduced in which coupling iterations are performed on periodic level instead of per time step. The convergence behaviour can be significantly improved by implementing existing partitioned solution methods as used for time step coupling (TSC) algorithms in the time periodic coupling (TPC) framework. The new algorithm has been evaluated by comparing the convergence behaviour to TSC algorithms. It is shown that the number of fluid–structure evaluations can be considerably reduced when a TPC algorithm is applied instead of a TSC. One of the most appealing advantages of the TPC approach is that the structural problem can be solved in the frequency domain resulting in a very efficient algorithm for computing steady state FSI responses.     

A new formulation of the effective elastic–plastic response of two-phase particulate composite materials

In this paper the double-inclusion model, originally developed to determine effective linear elastic properties of composite materials, is reformulated and extended to predict the effective nonlinear elastic–plastic response of two-phase particulate composites reinforced with spherical particles. The resulting problem of elastic–plastic deformation of a double-inclusion embedded in an infinite reference medium subjected to an incrementally applied far-field strain is solved by the finite element method. The proposed double-inclusion model is evaluated by comparison of the model predictions to the available exact results obtained by the direct approach using representative volume elements containing many particles. It is found that the double-inclusion formulation is capable of providing accurate prediction of the effective elastic–plastic response of two-phase particulate composites at moderate particle volume fractions.     

Dual finite element methods in homogenization for elastic–plastic fibrous composite material

The problem of homogenization for a periodic, elastic–perfectly plastic, fiber reinforced, composite material is considered. The overall mechanical behavior of the material is described using the anisotropic model of elastic–plastic body with kinematic hardening. The appropriate initial–boundary value problem, set for one repeatable cell of the composite, is solved in order to find the effective constitutive relations. The cell problem is solved using the finite element method formulated in two dual forms: in displacements and in stresses. Stress functions are used in the latter formulation.     

Effect of deviatoric nonassociativity on the failure prediction for elastic–plastic materials

This paper presents a theoretical study of the effect of nonassociativity of the plastic flow rule on the critical plastic modulus for discontinuous bifurcation in an elastic–plastic material. Nonassociativity in both the spherical and the deviatoric spaces are considered, with an emphasis on the effect of nonassociativity in the deviatoric space. A particular form of nonassociativity in the deviatoric space is introduced, where the projections of the plastic flow direction and the normal to the yield surface are assumed to have the same length but the projection of plastic flow direction is allowed to lag that of the normal by an angle. It is shown that even for the simple yield surface of von Mises, nonassociativity in the deviatoric space can lead to a bifurcation for a load parameter significantly lower than the value predicted with an associated flow rule.     

A critical state sand model with elastic–plastic coupling

Soil elastic moduli are highly pressure-dependent. Experimental findings have indicated that the elastic shear modulus of sands depends on p^χ, where p is mean principal effective stress and χ is a non-dimensional parameter. χ practically remains unchanged for shear strains less than 10⁻⁵ where the mechanical behavior is purely elastic. However, experiments have revealed that the emergence of plasticity for shear strains larger than 10⁻⁵ provokes a gradual increase in χ. Technically, this observation is an elastic–plastic coupling effect in which plasticity causes to change the elastic characteristics. Here, this issue is considered in hyper-elasticity framework in conjunction with a critical state compatible bounding surface plasticity platform for granular soils. To this aim, constitutive equations linking χ to a proper kinematic hardening parameter are presented. Then, using the proposed approach, a hyper-elastic theory is modified to consider the mentioned elastic–plastic coupling effect in the whole domain of the elastoplastic behavior. Adopting the improved hyper-elasticity necessitates the modification of a number of basic plasticity platform elements. In this regard, dilatancy and plastic hardening modulus of the bounding surface platform are modified. Successful performance of the modified constitutive model is presented against experimental data of loading/unloading triaxial tests.     

An elastic–plastic constitutive law for the description of uniaxial and multiaxial ratchetting

An elastic–plastic constitutive model is proposed to describe 1-D and 2-D ratchetting. The model is also able to give correct results for 2-D ratchetting when only uniaxial identification is used, while no special threshold or parameter is used for the case of non-proportional loading. The original feature of this model consist in the introduction of a ratchetting stress (material characteristic) along with the maximal stress supported in the history of loading and the plastic strain at the last unloading. In this paper uniaxial and 3-D formulations have been described based on a numerical implementation in the software Code_Aster. Uniaxial and also multiaxial identifications have been used. Simulations have been realized for proportional and non-proportional homogeneous cases, as well as for structures under anisothermal thermomechanical loading. The results of a benchmark on a structure, comparing experiment, simulations by this model and some other phenomenological models, and a polycrystalline model are presented. An analysis of error margin due to the choice of Mises criterion is exposed.     

Orthotropic elastic–plastic material model for paper materials

Paper and paperboard generally exhibit anisotropic and non-linear mechanical material behaviour. In this work, the development of an orthotropic elastic–plastic constitutive model, suitable for modelling of the material behaviour of paper is presented. The anisotropic material behaviour is introduced into the model by orthotropic elasticity and an isotropic plasticity equivalent transformation tensor. A parabolic stress–strain relation is adopted to describe the hardening of the material. The experimental and numerical procedures for evaluation of the required material parameters for the model are described. Uniaxial tensile testing in three different inplane material directions provides the calibration of the material parameters under plane stress conditions. The numerical implementation of the material model is presented and the model is shown to perform well in agreement with experimentally observed mechanical behaviour of paper.     

Micromorphic approach for gradient-extended thermo-elastic–plastic solids in the logarithmic strain space

The coupled thermo-mechanical strain gradient plasticity theory that accounts for microstructure-based size effects is outlined within this work. It extends the recent work of Miehe et al. (Comput Methods Appl Mech Eng 268:704–734, 2014) to account for thermal effects at finite strains. From the computational viewpoint, the finite element design of the coupled problem is not straightforward and requires additional strategies due to the difficulties near the elastic–plastic boundaries. To simplify the finite element formulation, we extend it toward the micromorphic approach to gradient thermo-plasticity model in the logarithmic strain space. The key point is the introduction of dual local–global field variables via a penalty method, where only the global fields are restricted by boundary conditions. Hence, the problem of restricting the gradient variable to the plastic domain is relaxed, which makes the formulation very attractive for finite element implementation as discussed in Forest (J Eng Mech 135:117–131, 2009) and Miehe et al. (Philos Trans R Soc A Math Phys Eng Sci 374:20150170, 2016).     

Spherical impression of thin elastic films on elastic–plastic substrates

The mechanical behavior of thin elastic films deposited onto structural alloys plays a critical role in determining film durability. This paper presents analysis of an impression experiment designed to evaluate some of the relevant properties of these films. The modeling provides quantitative strain information which can be used to estimate the fracture toughness of the film, the static friction coefficient of the surface and the constitutive behavior of the substrate. Results are presented for radial and circumferential strain distributions in the film relevant to the interpretation of cracking patterns. Additionally, load-displacement curves are provided that may be used to evaluate the plastic properties of the substrate. To facilitate estimates of the film cracking strain through correlation with experiments, the radial strain distributions are presented as functions of impression depth, yield strain and hardening exponent.     

Shakedown kinematic theorem for elastic–perfectly plastic bodies

The classical shakedown kinematic theorem due to Koiter for elastic–perfectly plastic bodies is re-examined and divided into separated shakedown and nonshakedown theorems. While the shakedown theorem is based on the set of Koiter's plastic strain rate cycles, the non-shakedown one involves a broader set of admissible plastic strain rate cycles, the end-cycle accumulated strains of which are deviatoric parts of compatible strain fields. For certain broad classes of practical problems the two statements are unified to yield the unique theorem in Koiter's sense.