An endochronic model of yield surface accounting for deformation induced anisotropy

In this article, an endochronic model of yield surface is proposed. Based on this model, the yield surface is simulated such that the forward and rear parts of the yield surface are described by different ellipses which are characterized by corresponding aspect ratio functions, respectively. Verification of the endochronic theory used the experimental results of yield surfaces obtained by Wu and Yeh for 304 stainless steel (Wu, H.C., Yeh, W.C., 1991. On the experimental determination of yield surfaces and some results of annealed 304 stainless steel. Int. J. Plasticity 7, 803–826). The experiments were performed cyclically under uniaxial, torsional, and combined axial–torsional loading conditions. The result has shown that the agreement between the prediction and experiments is quite satisfactory. In addition to the distortion of the yield surface plastically behaving a sharp front accompanied by a blunt rear, the anisotropic kinematic hardening effect has been addressed in this investigation. Although the experimental results of yield surfaces subjected to non-proportional loading conditions can be found in the literature, lack of information about the plastic strain history makes it impossible to verify the theory under such complicated loading conditions. The domain of applicability and validity of the theory, which is defined in terms of plastic strain increments, need be further investigated with the aim to set up related experiments.     

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.     

Directional distortional hardening in metal plasticity within thermodynamics

This paper presents a complete theory for metal plasticity that includes isotropic, kinematic, and directional distortional hardening, within the framework of thermodynamics. Directional distortion is defined here as the formation of a region of high curvature on the yield surface, approximately in the direction of loading, and a region of flattening approximately in the opposite direction, as observed in experiments on various types of metals. The distinguishing features of this theory are the introduction of a fourth order tensor-valued internal variable, whose evolution in conjunction with a directional scalar multiplier describes the evolving directional distortion, and the fact that the hardening laws for all internal variables are derived on the basis of sufficient conditions to satisfy the thermodynamic requirement of positive dissipation. The applicability of the theory is illustrated by fitting experimental data on distorted yield surfaces in the course of plastic deformation.     

Calcul à la rupture en présence d'un écoulement : formulation cinématique avec un champ de pression approchéYield design formulation for porous media subjected to flow,using approximate pressure field

We attempt here to use the kinematic method of yield design in the case of a porous medium subjected to flow (with or without free surface), without looking for the exact solution of the pressure field. The method proposed here is based on the use of approximate pressure fields. In this paper, we show how, under different conditions concerning the yield criterion and the velocity field, the use of such approximate fields allows one to obtain a necessary condition for stability without having to find the real pressure field. To cite this article: A. Corfdir, C. R. Mecanique 334 (2006).     

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.     

Evolution of subsequent yield surfaces and elastic constants with finite plastic deformation. Part II: A very high work hardening aluminum alloy (annealed 1100 Al)

Results are presented on the evolution of subsequent yield surfaces with finite deformation in a very high work hardening annealed 1100 aluminum alloy. In Part I [Khan, A.S., Kazmi, R., Stoughton, T., Pandey, A., 2009a. Evolution of subsequent yield surfaces and elastic constants with finite plastic deformation. Part 1: a very low work hardening aluminum alloy (Al-6061–T6511) 25, 1611–1625.] of this paper, similar results are presented for a very low work hardening aluminum alloy. Those results were very different from the present ones, and all the results were for proportional loading paths. The subsequent yield surfaces are determined in tension, free end torsion and combined tension–torsion proportional and non-proportional loading paths, using 10 με deviation from linearity definition of yield. Yield surfaces are also determined after linear, bi-linear, and non-linear unloading paths after finite deformation under tension, free end torsion, and combined tension–torsion loading. The initial yield surface is closer to the von-Mises surface and the subsequent yield surfaces show distortion, expansion, positive cross-effect, and “nose” in the loading direction. Additionally, the subsequent yield surfaces after non-proportional loading paths show shrinkage and compounded distortion. The yield surfaces after unloading depict strong anisotropy, positive cross-effect and exhibits different proportion of distortion in each loading conditions. The Young’s and shear modulus decrease with plastic deformation and this decrease is much less than those reported in the published literature.     

A weakened form of Ilyushin's postulate and the structure of self-consistent Eulerian finite elastoplasticity

From a general standpoint in terms of internal variables, we formulate a general theory of self-consistent Eulerian finite elastoplasticity based on the additive decomposition of the Eulerian strain rate, i.e., D=D^e+D^p, as well as two consistency criteria. In this theory, the elastic behaviour is characterized by an exactly integrable elastic rate equation for D^e with a general form of complementary elastic potential. It is assumed that the yield function depends in a general manner on the Kirchhoff stress and the internal variables. Moreover, the plastic rate equation for D^p and the evolution equation for each internal variable are allowed to assume general forms relying on the just-mentioned variables and the stress rate. It is indicated that two consistency criteria, i.e., the self-consistency for the elastic rate equation and Prager's yielding stationarity, lead to the unique choice of objective rates, i.e., the logarithmic rate.The structure of the above theory is further studied and examined by virtue of a weakened form of Ilyushin's postulate. In a spinning frame defining the logarithmic rate, we introduce the notion of standard elastoplastic strain cycle, which starts at a point not on but inside a yield surface and incorporates only one infinitesimal plastic subpath. We show that this type of strain cycle is always possible. Then, by ruling out strain cycles starting at points on yield surfaces we propose a weakened form of Ilyushin's postulate, which says that the changing rate of the stress work done along every standard strain cycle should be non-negative, whenever the incorporated plastic subpath tends to vanish. By virtue of simple, rigorous procedures, we demonstrate that this weakened form of Ilyushin's postulate is adequate to ensure direct results concerning the normality rule and the convexity of the yield surface in the context of the foregoing Eulerian finite elastoplasticity theory. Specifically, with an exactly integrable elastic rate equation defining D^e, we prove that, in the space of the Kirchhoff stresses, the difference (D−D^e) is just the gradient of the yield function multiplied by a plastic multiplier, and thus bears the very kinematical and physical feature of plastic strain rate. Furthermore, we prove that, in the space of the Kirchhoff stresses, the elastic domain bounded by each yield surface should be convex. The main results are derived in a self-contained manner within the context of an Eulerian theory of finite elastoplasticity, without involving issues concerning how to define intermediate stress-free states and plastic strains, etc.     

Yield criteria for porous media in plane strain: second-order estimates versus numerical resultsCritère de plasticité des matériaux poreux en déformation plane : Estimations du second ordre et résultats numériques

This Note presents a comparison of some recently developed “second-order” homogenization estimates for two-dimensional, ideally plastic porous media subjected to plane strain conditions with corresponding yield analysis results using a new linearization technique and systematically optimized finite elements meshes. Good qualitative agreement is found between the second-order theory and the yield analysis results for the shape of the yield surfaces, which exhibit a corner on the hydrostatic axis, as well as for the dependence of the effective flow stress in shear on the porosity, which is found to be non-analytic in the dilute limit. Both of these features are inconsistent with the predictions of the standard Gurson model. To cite this article: J. Pastor, P. Ponte Castañeda, C. R. Mecanique 330 (2002) 741–747.     

Chargements axisymétriques d'un bicouche transversalement isotropeAxisymmetric loading on a transversely isotropic two-layered medium

An exact three-dimensional analysis is developped for an axisymmetric loading acting on the surface of a semi-infinite medium composed by two transversely isotropic materials. The loading is assumed to be parallel to the elastic symmetry axis of the upper layer. The solutions of a concentrated force and a uniform loading distributed on a circle are obtained by exact integral expressions. The numerical results are performed to show the anisotropic effect with isovalue curves of stress. To cite this article: C. Ruimy, M. Dahan, C. R. Mecanique 330 (2002) 469–473.     

Structure des développements de distorsion rapide à petits temps en turbulence homogène

For small times following the distortion of an isotropic state, rapid-distortion theory provides the tensorial form of one-point correlation expansions for homogeneous rotational turbulent mean flows. It is considered that the consistency with such expansions must be satisfied by any closure model. The Note describes the general structure of these expansions as well as some of their properties. It is shown how cumulated effects (strain and rotation) are involved. To cite this article: J. Piquet, C. R. Mecanique 333 (2005).     

Analyses of steady quasistatic crack growth in plane strain tension in elastic-plastic materials with non-isotropic hardening

Steady state crack propagation problems of elastic-plastic materials in Mode I, plane strain under small scale yielding conditions were investigated with the aid of the finite element method. The elastic-perfectly plastic solution shows that elastic unloading wedges subtended by the crack tip in the plastic wake region do exist and that the stress state around the crack tip is similar to the modified Prandtl fan solution. To demonstrate the effects of a vertex on the yield surface, the small strain version of a phenomenological J₂, corner theory of plasticity (Christoffersen, J. and Hutchinson, J. W. J. Mech. Phys. Solids,27, 465 C 1979) with a power law stress strain relation was used to govern the strain hardening of the material. The results are compared with the conventional J₂ incremental plasticity solution. To take account of Bauschinger like effects caused by the stress history near the crack tip, a simple kinematic hardening rule with a bilinear stress strain relation was also studied. The results are again compared with the smooth yield surface isotropic hardening solution for the same stress strain curve. There appears to be more potential for steady state crack growth in the conventional J₂ incremental plasticity material than in the other two plasticity laws considered here if a crack opening displacement fracture criterion is used. However, a fracture criterion dependent on both stress and strain could lead to a contrary prediction.     

Multiaxial yield behavior of 1100 aluminum following various magnitudes of prestrain

The multiaxial yield and flow behavior of metals has been of interest for many years. Recently, the experimental work of Phillips & Lee [1979], Shiratori et al. [1979] and Ohashi [1982] has been quite notable in this field. These authors have concentrated their efforts in measuring yield loci after small to moderate prestrains (≤0.06). In this paper we discuss small strain yield loci we have measured after prestrains between 0.03 and 0.05 in torsion. These experiments on 1100 aluminum are in general agreement with the literature. They show a translation, distortion and expansion of the yield loci. A rounded nose forms in the direction of prestrain with the yield locus flattering opposite the prestrain. We observed that the distortions change to match the strain direction after very small reversals in prestrain.The subsequent yield locus has also been measured after a large torsional prestrain of γ=0.5. Using a 5 × 10^?6 offset criterion for yielding, the shape, distortion and translation of the yield locus was very similar to that found after the smaller prestrains. In addition a large-strain yield locus, using a back extrapolation technique, was determined for the same sample. This yield locus exhibited close to von Mises isotropic expansion. The observed deviations, while slight are extremely important. They match those predicted by a polycrystal slip model. Thus, the small-strain yield locus, after a large prestrain, appears to be determined largely from dislocation considerations only, where as the large-strain yield locus is determined by the developing texture. Finally, aluminum sheet was deformed by rolling to larger prestrains ?^{von Mises} = 0.5, 1.0, 1.5, 2.0 and 2.5 and subsequently tested in plane strain compression. Two types of compression experiments were done, one such that there was no deformation mode change from rolling, the other rotating the direction of zero strain by 90° producing a stress path change. The large strain yield and flow behavior of these experiments was again predicted using the relaxed constraint polycrystal model of Kocks & Canova [1981]. For these very large prestrains the experiments and texture theory differ. Micrstructural observations have shown the presence of micro-shear bands which resulted from the rolling prestrain. We speculate that these features are responsible for the deviation from crystal plasticity theory.We believe that this work points to several operative mechanisms of deformation. Small-strain yielding (5 × 10^?6) appears to be controlled purely by dislocation mechanisms and interactions even after relatively large prestrains. Large-strain yielding, on the other hand, is controlled by texture after moderate prestrains (at least to γ = 0.5). After large prestrains, obtained by rolling, the experiments deviate from texture based predictions. This is possibly the result of microstructural deformation mechanisms, for example micro-shear bands, playing a role in the deformation process.     

On loading,yield and quasi-yield hypersurfaces in plasticity theory

In the context of the author's previously published “simple” theory of plasticity[1] in which no loading or yield surfaces are assumed to exist, it is shown that (a) loading surfaces must exist for a plastic material as a result of Caratheodory's theorem on Pfaffian forms, and that (b) a yield hypersurface in state space may be defined as the boundary of the region in which no loading surfaces exist (the elastic region) if this region has a positive volume, otherwise this region degenerates into the quasi-yield hypersurface. The significance of loading and yield (or quasi-yield) hypersurfaces is further explored for one-component loadings, with particular attention to the Bauschinger effect and kinematic hardening.     

L'approche variationnelle de la fatigue : des premiers résultats

By using a principle of least energy and a Dugdale surface energy with an irreversibility condition, we build a debonding model of thin films valid both for monotone and cyclic loading. We show that, if the internal length introduced in Dugdale model is small in comparison to the film length, then the growth of the debonding follows Griffith's law under monotone loading and a Paris-type law under cycling loading. To cite this article: A. Jaubert, J.-J. Marigo, C. R. Mecanique 333 (2005).     

A new approach to the pseudo-orthogonal properties of eigenfunction expansion form of the crack-tip complex potential function in anisotropic and piezoelectric fracture mechanics

Over the past twenty years, the well-known weight function theory based on the Bueckner work conjugate integral has been widely used to calculate crack tip fracture dominant parameter such as the stress intensity factor, the energy release rate (or the J-integral) and the T-stress in various kinds of cracked materials (e.g. isotropic materials, anisotropic materials and piezoelectric materials). Meanwhile, the pseudo-orthogonal property of the eigenfunction expansion form of the crack tip stress complex potential function has been proved to play a very important role in the theory. In this paper, we provide a new approach to establish the pseudo-orthogonal properties for crack problems in anisotropic and/or piezoelectric materials. In the latter case associated with mechanical-electric coupling, the electrical boundary conditions under both impermeable and permeable crack models are considered. The approach developed is much simpler than the classical complex variable separation technique proposed by previous researchers and hence the cumbersome and lengthy manipulations are avoided. Moreover, it is shown that, unlike previous works, the orthogonal properties of the material characteristic matrices A and B induced by the Stroh theory are no longer necessary in establishing the pseudo-orthogonal properties of eigenfunction expansion form in cracked piezoelectric materials. The approach can be easily extended to treat many other different crack problems concerning the Bueckner integral involving several complex arguments.     

A superposition principle in optimal plastic design for alternative loads

The paper is concerned with the optimal plastic design of sandwich beams, frames and trusses for alternative loading conditions. Upper and lower bounds for the optimal weight of a beam are derived, for single as well as for alternative loading conditions. These bounding theorems are used to establish a superposition principle. If no explicit bounds on the cross-sectional areas are prescribed, the optimal design for alternative loading conditions P₁ and P₂ can be obtained by superposition of the optimal designs for the single loading conditions and . If the cross-sections are to have at least given non-zero values, the principle furnishes upper and lower bounds to the optimal weight.The principle is illustrated by a simple example.     

Représentation à l'échelle microscopique d'un mécanisme de rupture dans un milieu a joints

. A representation of a class of failure mechanisms defined on a homogenized jointed medium by a velocity jump along a planar surface, is analyzed at the scale of the joints. The approach consists in representing a velocity jump at the macroscopic level by the kinematics of a finite band in the heterogeneous medium characterized by a piecewise rigid body motion of elementary blocs. The equivalence (in the sense of the yield design theory) between the macroscopic and macroscopic representations of such failure mechanisms is ensured by the equality of the corresponding maximum resisting rates of work. To cite this article: S. Maghous et al., C. R. Mecanique 333 (2005).     

On loading criteria in plasticitySur des critères de chargement en plasticité

The loading criteria of the Lagrangian strain-space formulation of rate-independent plasticity are compared with those of Nguyen and Bui and those of Kuhn–Tucker type. When the latter two sets of conditions are expressed in a fully strain-space form, their relationship to the loading criteria of the strain-space formulation becomes transparent. To cite this article: J. Casey, C. R. Mecanique 330 (2002) 285–290.     

Modified kinematic hardening rule for multiaxial ratcheting prediction

A modified kinematic hardening rule is proposed in which one biaxial loading dependent parameter δ′ connecting the radial evanescence term [(α:n)ndp] in the Burlet–Cailletaud model with the dynamic recovery term of Ohno–Wang kinematic hardening rule is introduced into the framework of the Ohno–Wang model. Compared with multiaxial ratcheting experimental data obtained on 1Cr18Ni9Ti stainless steel in the paper and CS1026 steel conducted by Hassan et al. [Int. J. Plasticity 8 (1992) 117], simulation results by modified model are quite well in all loading paths. The simulations of initial nonlinear part in ratcheting curves can be improved greatly while the evolutional parameter δ′ related to plastic strain accumulation is added into the modified model.