Effective viscoplastic behavior of polycrystalline aggregates lacking four independent slip systems inferred from homogenization methods; application to olivine

Polycrystalline aggregates lacking four independent systems for the glide of dislocations can deform in a purely viscoplastic regime only if additional deformation mechanisms (such as grain boundary sliding and diffusion) are activated. We introduce an implementation of the self-consistent scheme in which this additional physical mechanism, considered as a stress relaxation mechanism, is represented by a nonlinear isotropic viscoplastic potential. Several nonlinear extensions of the self-consistent scheme, including the second-order method of Ponte-Castañeda, are used to provide an estimate of the effective viscoplastic behavior of such polycrystals. The implementation of the method includes an approximation of the isotropic potential to ensure convergence of the attractive fixed-point numerical algorithm. The method is then applied to olivine polycrystals, the main constituent of the Earth's upper mantle. Due to the extreme local anisotropy of the local constitutive behavior and the subsequent intraphase stress and strain-rate field heterogeneities, the second-order method is the only extension providing qualitative and quantitative accurate results. The effective viscosity is strongly dependent on the strength of the relaxation mechanism. For olivine, a linear viscous relaxation (e.g. diffusion) could be relevant; in that case, the polycrystal stress sensitivity is reduced compared to that of dislocation glide, and the most active slip system is not necessarily the one with the smallest reference stress due to stress concentrations. This study reveals the significant importance of the strength and stress sensitivity of the additional relaxation mechanism for the rheology and lattice preferred orientation in such highly anisotropic polycrystalline aggregates.     

Rate-independent crystalline and polycrystalline plasticity,application to FCC materials

This paper deals with the simulation of the mechanical response and texture evolution of cubic crystals and polycrystals for a rate-independent elastic–plastic constitutive law. No viscous effects are considered. An algorithm is introduced to treat the difficult case of multi-surface plasticity. This algorithm allows the computation of the mechanical response of a single crystal. The corresponding yield surface is made of the intersection of several hyper-planes in the stress space. The problem of the multiplicity of the slip systems is solved thanks to a pseudo-inversion method. Self and latent hardening are taken into account. In order to compute the response of a polycrystal, a Taylor homogenization scheme is used. The stress–strain response of single crystals and polycrystals is computed for various loading cases. The texture evolution predicted for compression, plane strain compression and simple shear are compared with the results given by a visco-plastic polycrystalline model.     

Boundary film shear elastic modulus effect in hydrodynamic contact. Part II: influence of operational parameters

The present paper is the subsequent research of the first part (Theor Comput Fluid Dyn, 2009). It investigates the boundary film shear elastic modulus effect in a hydrodynamic contact in different operating conditions. The hydrodynamic contact is one-dimensional, composed of two parallel plane surfaces, which are respectively rough rigid with rectangular micro projections in profile periodically distributed on the surface and ideally smooth rigid. The whole contact consists of cavitated area and hydrodynamic area. The hydrodynamic area consists of many micro Raleigh bearings which are discontinuously and periodically distributed in the contact. The hydrodynamic contact in a micro Raleigh bearing consists of boundary film area and fluid film area which, respectively, occur in the outlet and inlet zones. In boundary film area, the film slips at the upper contact surface due to the limited shear stress capacity of the film–contact interface, while the film does not slip at the lower contact surface due to the shear stress capacity of the film–contact interface large enough. In boundary film area, the viscosity, density, and shear elastic modulus of the film are varied across the film thickness due to the film–contact interactions, and their effective values are used in modeling which depends on the film thickness. In fluid film area, the film does not slip at either of the contact surfaces, and the shear elastic modulus of the film is neglected. It is found from the simulation results that the boundary film shear elastic modulus influences are normally negligible on the mass flow through the contact, the carried load of the contact and the overall film thickness of the contact, and the boundary film shear elastic modulus would normally influence the local film thickness in an elastic contact when the local film thickness is on the film molecule diameter scale. It is also found that the boundary film shear elastic modulus effect has the tendency of being increased with the reduction of the width of a micro contact. It is increased with the reduction of the boundary film–contact interfacial shear strength or with the increase of the critical boundary film thickness, while it is strongest at certain values of the contact surface roughness, the width ratio of fluid film area to boundary film area, and the lubricant film shear elastic modulus.

---

     

A multiscale approach for modeling scale-dependent yield stress in polycrystalline metals

Modeling of scale-dependent characteristics of mechanical properties of metal polycrystals is studied using both discrete dislocation dynamics and continuum crystal plasticity. The initial movements of dislocation arc emitted from a Frank-Read type dislocation source and bounded by surrounding grain boundaries are examined by dislocation dynamics analyses system and we find the minimum resolved shear stress for the FR source to emit at least one closed loop. When the grain size is large enough compared to the size of FR source, the minimum resolved shear stress levels off to a certain value, but when the grain size is close to the size of the FR source, the minimum resolved shear stress shows a sharp increase. These results are modeled into the expression of the critical resolved shear stress of slip systems and continuum mechanics based crystal plasticity analyses of six-grained polycrystal models are made. Results of the crystal plasticity analyses show a distinct increase of macro- and microscopic yield stress for specimens with smaller mean grain diameter. Scale-dependent characteristics of the yield stress and its relation to some control parameters are discussed.     

Collective micro shear processes and plastic instabilities in crystalline and amorphous structures

The mechanisms that may lead to plastic instabilities are reviewed briefly, and their intimate connection with localization of plastic shear is emphasized. The combination of macro and micro experiments, ie, analysis of stress fluctuations and serrations, and of slip and shear band evolution during deformation, for investigations of these mechanisms is demonstrated with examples of LÜDERS and PORTEVIN-LeCHATELIER band propagation in single glide oriented crystals, and of shear band formation in predeformed mono- and polycrystals as well as amorphous materials. The problems connected with collective effects of the micro units of deformation (ie, dislocations in crystalline, local shear transformations in amorphous materials) are indicated, and the necessity to study the mechanism of transfer of slip (shear) to neighbouring regions is stressed.     

Homogenization estimates for the average behavior and field fluctuations in cubic and hexagonal viscoplastic polycrystals

The “second-order” homogenization procedure (J. Mech. Phys. Solids 50 (2002) 737) is used to compute estimates of the self-consistent type for the effective response of cubic and hexagonal viscoplastic polycrystals with isotropic textures. The method, which requires the computation of the averages of the stress field and the covariances of its fluctuations over the various grain orientations in an optimally selected “linear comparison polycrystal,” is also used to generate information on the heterogeneity of the stress and strain-rate fields within the polycrystals. In contrast with earlier estimates of the self-consistent type, such as those arising from the “incremental” and “tangent” schemes, the new estimates for the effective behavior are found to satisfy all known bounds, even in the strongly nonlinear, rate-insensitive limit. In addition, they are found to satisfy a recently proposed scaling law at large grain anisotropy. The fluctuations of the stresses and strain rates, which are nonzero for all grain orientations, are found to generally increase with decreasing strain-rate sensitivity (i.e., increasing nonlinearity) and with increasing grain anisotropy (which is typically higher for lower-symmetry systems).     

Ultrasonic attenuation in polycrystals using a self-consistent approach

The self-consistent method of averaging elastic moduli to define the effective medium of a polycrystal is used to investigate the dynamic problem of wave propagation. An alternative covariance tensor describing the elastic moduli fluctuations of the polycrystal containing self-consistent elastic properties is derived and found to be significantly smaller than the covariance tensor formed through traditional Voigt averaging. Attenuation curves are generated using the self-consistent elastic moduli and covariance tensors and these results are compared with previous Voigt-averaged estimates. The second-order polycrystalline dispersion relation for the self-consistent scheme is found for cases of low and high crystallite anisotropy. The attenuation coefficients and dispersion relations derived through the self-consistent scheme are considerably different than previous estimates. Experimentally measured longitudinal attenuation coefficients support the use of the self-consistent scheme for estimation of attenuation.     

AN ANALYSIS OF THE EFFECTS OF SYMMETRICAL TILT GRAIN BOUNDARIES ON THE PLASTIC DEFORMATION

A numerical investigation on the simple polycrystals containing threesymmetrical tilt grain boundaries(GBs)is carried out within the framework of crystalplasticity which precisely considers the finite deformation and finite lattice rotation aswell as elastic anisotropy.The calculated results show that the slip geometry and theredistribution of stresses arising from the anisotropy and boundary constraint play animportant role in the plastic deformation in the simple polycrystals.The stress levelalong GB is sensitive to the load level and misorientation,and the stresses along GB aredistributed nonuniformly.The GB may exhibit a softening or strengthening feature,which depends on the misorientation angle.The localized deformation bands usuallydevelop accompanying the GB plastic deformation,the impingement of the localizedband on the GB may result in another localized deformation band.The yield stresseswith different misorientation angles are favorably compared with the experimentalresults.     

边界膜剪切弹性模量对于微接触性能的影响

建立接触模型,理论分析了微接触中边界膜剪切弹性模量对于接触性能的影响. 接触区由两平行平面形成,属一维接触. 上接触表面为粗糙表面,具有矩形微凸体. 下接触表面为光滑平面. 两接触表面均处理成刚性表面. 微接触区中充满流体. 它分成两个子区,在微接触的出口区由于极小的接触间隙充满边界膜,在微接触的入口区由于接触间隙较大充满流体膜. 边界膜和流体膜行为决定整个微接触性能. 当膜厚较大时,这里边界膜可看成纳米级薄膜. 由于上接触表面处有限的剪应力承受能力,边界膜可于上接触表面滑移. 设下接触表面处剪应力承受能力很大而边界膜在下接触表面不滑移. 由于边界膜-接触表面间相互作用,边界膜黏度、密度和剪切弹性模量均沿膜厚变化,在理论分析中使用它们的等效值,这些值与边界膜厚度有关. 流体膜在两个接触表面均不发生滑移,分析中不考虑流体膜剪切弹性模量. 流体膜采用传统分析法. 给出了理论分析和若干变工况参数下的计算结果.      

Estimation of overall properties of random polycrystals with the use of invariant decompositions of Hooke’s tensor

In the paper the theoretical analysis of bounds and self-consistent estimates of overall properties of linear random polycrystals composed of arbitrarily anisotropic grains is presented. In the study two invariant decompositions of Hooke’s tensors are used. The applied method enables derivation of novel expressions for estimates of the bulk and shear moduli, which depend on invariants of local stiffness tensor. With use of these expressions the materials are considered for which at the local level constraints are imposed on deformation or some stresses are unsustained.     

Boundary film shear elastic modulus effect in hydrodynamic contact. Part I: theoretical analysis and typical solution

Boundary film shear elastic modulus effect is analyzed in a hydrodynamic contact. The contact is one-dimensional composed of two parallel plane surfaces, which are, respectively, rough rigid with rectangular micro projections in profile periodically distributed on the surface and ideally smooth rigid. The whole contact is consisted of cavitated area and hydrodynamic area. The hydrodynamic area consists of many micro Raleigh bearings which are discontinuously and periodically distributed in the contact. Analysis is thus carried out for a micro Raleigh bearing in this contact. The hydrodynamic contact in this micro Raleigh bearing consists of boundary film area and fluid film area which, respectively, occur in the outlet and inlet zones. In boundary film area, the film slips at the upper contact surface due to the limited shear stress capacity of the film–contact interface, while the film does not slip at the lower contact surface due to the shear stress capacity large enough at the film–contact interface. In boundary film area, the viscosity, density and shear elastic modulus of the film are varied across the film thickness due to the film–contact interactions, and their effective values are used in modeling, which depend on the film thickness. The analytical approach proposed by Zhang (J Mol Liq 128:60–64, 2006) and Zhang et al. (Int J Fluid Mech Res 30:542–557, 2003) is used for boundary film area. In fluid film area, the film does not slip at either of the contact surfaces, and the shear elastic modulus of the film is neglected. Conventional hydrodynamic analysis is used for fluid film area. The present paper presents the theoretical analysis and a typical solution. It is found that for the simulated case the boundary film shear elastic modulus effects on the mass flow through the contact, the overall film thickness of the contact and the carried load of the contact are negligible but the boundary film shear elastic modulus effect on the local film thickness of the contact may be significant when the boundary film thickness is on the 1 nm scale and the contact surfaces are elastic. In Part II will be presented detailed results showing boundary film shear elastic modulus effects in different operating conditions.

---

     

Bounds and self-consistent estimates for elastic constants of random polycrystals with hexagonal, trigonal, and tetragonal symmetries

Peselnick, Meister, and Watt have developed rigorous methods for bounding elastic constants of random polycrystals based on the Hashin-Shtrikman variational principles. In particular, a fairly complex set of equations that amounts to an algorithm has been presented previously for finding the bounds on effective elastic moduli for polycrystals having hexagonal, trigonal, and tetragonal symmetries. A more analytical approach developed here, although based on the same ideas, results in a new set of compact formulas for all the cases considered. Once these formulas have been established, it is then straightforward to perform what could be considered an analytic continuation of the formulas (into the region of parameter space between the bounds) that can subsequently be used to provide self-consistent estimates for the elastic constants in all cases. This approach is very similar in spirit but differs in its details from earlier work of Willis, showing how Hashin-Shtrikman bounds and certain classes of self-consistent estimates may be related. These self-consistent estimates always lie within the bounds for physical choices of the crystal elastic constants and for all the choices of crystal symmetry considered. For cubic symmetry, the present method reproduces the self-consistent estimates obtained earlier by various authors, but the formulas for both bounds and estimates are generated in a more symmetric form. Numerical values of the estimates obtained this way are also very comparable to those found by the Gubernatis and Krumhansl coherent potential approximation (or CPA), but do not require computations of scattering coefficients.     

Bounds on the Elastic Moduli of Completely Random Two-Dimensional Polycrystals

Explicit bounds on the elastic moduli of completely random planar polycrystals, the shape and crystalline orientations of the constituent grains of which are uncorrelated, are derived and calculated for a number of crystals of general two-dimensional anisotropy. The bounds on the elastic two-dimensional bulk modulus happen to coincide with the simple third order (in anisotropy contrast) bounds for the subclass of idealistic circular cell polycrystals. The bounds on the shear modulus are close to the much simpler bounds for circular cell polycrystals, which approximate aggregates of equiaxed grains.     

Study of the antiplane deformation of linear 2-D polycrystals with different microstructures

The effective behavior and the distribution of local mechanical fields of linearly viscous 2-D polycrystals under antiplane shear is investigated. Several microstructures are considered, and a full-field approach based on the Fast Fourier Transform technique is applied. First, the accuracy of this technique is evaluated on a strictly isotropic 2-phase microstructure. Voronoi tessellation is then used to generate artificial microstructures, and a real (fully recrystallized) polycrystalline microstructure is obtained by electron back-scattering diffraction. Ensemble averages over several configurations using eight crystalline orientations (phases) are performed. Although a slight anisotropy is obtained for the effective behavior of each individual configuration, statistically, the results are in very good agreement with the available analytical isotropic solution. At phase level, a marked asymmetry is obtained for the distribution of local stresses. The intraphase first- and second-order moments of the stress field, calculated for both microstructures are compared with corresponding self-consistent predictions.     

Evaluation of self-consistent polycrystal plasticity models for magnesium alloy AZ31B sheet

Various self-consistent polycrystal plasticity models for hexagonal close packed (HCP) polycrystals are evaluated by studying the deformation behavior of magnesium alloy AZ31B sheet under different uniaxial strain paths. In all employed polycrystal plasticity models both slip and twinning contribute to plastic deformation. The material parameters for the various models are fitted to experimental uniaxial tension and compression along the rolling direction (RD) and then used to predict uniaxial tension and compression along the traverse direction (TD) and uniaxial compression in the normal direction (ND). An assessment of the predictive capability of the polycrystal plasticity models is made based on comparisons of the predicted and experimental stress responses and R values. It is found that, among the models examined, the self-consistent models with grain interaction stiffness halfway between those of the limiting Secant (stiff) and Tangent (compliant) approximations give the best results. Among the available options, the Affine self-consistent scheme results in the best overall performance. Furthermore, it is demonstrated that the R values under uniaxial tension and compression within the sheet plane show a strong dependence on imposed strain. This suggests that developing anisotropic yield functions using measured R values must account for the strain dependence.     

Higher-order stress and grain size effects due to self-energy of geometrically necessary dislocations

The higher-order stress work-conjugate to slip gradient in single crystals at small strains is derived based on the self-energy of geometrically necessary dislocations (GNDs). It is shown that this higher-order stress changes stepwise as a function of in-plane slip gradient and therefore significantly influences the onset of initial yielding in polycrystals. The higher-order stress based on the self-energy of GNDs is then incorporated into the strain gradient plasticity theory of Gurtin [2002. A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 50, 5-32] and applied to single-slip-oriented 2D and 3D model crystal grains of size D. It is thus found that the self-energy of GNDs gives a D^-1-dependent term for the averaged resolved shear stress in such a model grain under yielding. Using published experimental data for several polycrystalline metals, it is demonstrated that the D^-1-dependent term successfully explains the grain size dependence of initial yield stress and the dislocation cell size dependence of flow stress in the submicron to several-micron range of grain and cell sizes.     

A micromechanical theory of grain-size dependence in metal plasticity

The effect of grain-size on the elastoplastic behavior of metals is investigated from the micromechanics standpoint. First, based on the observations that dislocation pile-ups, formation of cell structures, and other inelastic activities influenced by the presence of grain boundary actually take place transcrystallinely, a grain-size dependent constitutive equation is proposed for the slip deformation of slip systems. By means of a modified Hill's self-consistent relation the local stress of a grain is calculated, and used in conjunction with this constitutive equation to evaluate the plastic strain of each constituent grain. The grain-size effect on the plastic flow of polycrystals then can be determined by an averaging process. To check the validity of the proposed theory it was finally applied to predict the stress-strain curves and flow stresses of a copper at various grain-sizes. The obtained results were found to be in good agreement with experimental data.     

A study of directionality of latent hardening behavior in aluminium single crystals

By defining the yield stress in a latent hardening test as the first deviation from the elastic straight line, the yielding and hardening behavior on a latent system in the positive and negative slip direction was studied in aluminium single crystals. It is shown that the yield stresses on both the positive and negative latent systems are about equal to or a little lower than the maximum resolved shear stress in the primary test, but much higher than that of the active system. The influence of relative orientation and prestrain on latent hardening and initial work-hardening in the secondary test was also investigated, and it was found that there is a considerable effect on initial work-hardening, but none on latent hardening. With reasonable approximation, a hardening rule for single crystal could be proposed from the experimental results, that is, except for the yield stress on the system negative to the active system that is very low, hardening on the other systems is nearly the same as self-hardening.     

铝单晶体潜在硬化的指向性行为研究

应用高纯铝单晶体，采用偏离弹性线法定义其在潜在滑移系统在屈服应力，研究了潜在滑移系在正负两个滑移方向上的屈服及应变硬化行为。结果表明，潜在滑移系负行为的差异要远远小于动滑移系的Ｂａｕｓｃｈｉｎｇｅｒ效应，其正负方向的屈服应力一般均等于或略小于预应变时的最大分切应力，大大高于动滑移系的负向屈服应力，潜在滑移系和原始滑移系的相对取向及预应变的大小对单晶体潜在移系在潜在硬化的影响不大，但对其起始过渡区应     

A constitutive model for orthotropic elasto-plasticity at large strains

Summary The paper presents a thermodynamically consistent constitutive model for elasto-plastic analysis of orthotropic materials at large strain. The elastic and plastic anisotropies are assumed to be persistent in the material but the anisotropy axes can undergo a rigid rotation due to large plastic deformations. The orthotropic yield function is formulated in terms of the generally nonsymmetric Mandel stress tensor such that its skew-symmetric part is additionally taken into account. Special attention is focused on the convexity of the yield surface resulting in the nine-dimensional stress space. Of particular interest are new convexity conditions which do not appear in the classical theory of anisotropic plasticity. They impose additional constraints on the material constants governing the plastic spin. The role of the plastic spin is further studied in simple shear accompanied by large elastic and large plastic deformations. If the plastic spin is neglected, the shear stress response is characterized by oscillations with an amplitude strictly dependent on the degree of the plastic anisotropy.accepted for publication 2 March 2004