A Dedicated DIC Methodology for Characterizing Plastic Deformation in Single Crystals

This paper focuses on the development of an appropriate Digital Image Correlation (DIC) methodology based on Image Registration and dedicated for characterizing the plastic deformation in single crystals. A pure nickel single crystal specimen is plastically deformed in tension and investigated by DIC technique. Based on the measured kinematic fields, the proposed method enables to identify the slip activity on the crystal surface and to locate precisely the slip band interfaces at microscale which behave as kinematic discontinuities. The computed displacement data are projected on a well-defined physical basis containing slip details, then the strain fields can be derived directly from a set of analytical functions. The possible errors in displacement induced by this projection approach are evaluated. Finally, some results of the evaluated strain fields are presented. It demonstrates that the developed DIC methodology allows quantitative characterization of a heterogeneous deformation process and promotes further relationships to be established between slip activity and strain field evolution in single crystals.     

位移幅值对铜镁合金微动磨损行为的影响

在研制的多功能微动磨损试验机上,开展了不同位移幅值下铜镁合金微动磨损试验,以研究位移幅值对铜镁合金微动磨损行为的影响. 微动过程中记录摩擦系数曲线与F_t-D-N曲线,利用光学显微镜(OM)、扫描电镜(SEM)、能量色散X射线光谱仪(EDS)及三维形貌仪对损伤区域进行了微观分析. 结果显示:随着位移幅值的增加,铜镁合金微动运行状态由部分滑移进入完全滑移,未发现混合滑移状态;部分滑移区中呈现由弹性变形协调逐渐向塑性变形协调转变的趋势. 磨损体积随位移幅值的增加而增加,在完全滑移区中体积损失非常严重. 在弹性变形协调的部分滑移状态下,接触表面损伤轻微,而由塑性变形协调的部分滑移状态下,接触中心出现较大切应力,疲劳裂纹扩展至接触表面导致材料剥落,接触边缘有磨粒磨损和氧化磨损的痕迹. 在完全滑移状态下,接触表面损伤主要为疲劳剥层,磨粒磨损和氧化磨损.      

Non-convex rate dependent strain gradient crystal plasticity and deformation patterning

A rate dependent strain gradient crystal plasticity framework is presented where the displacement and the plastic slip fields are considered as primary variables. These coupled fields are determined on a global level by solving simultaneously the linear momentum balance and the slip evolution equation, which is derived in a thermodynamically consistent manner. The formulation is based on the 1D theory presented in Yalcinkaya et al. (2011), where the patterning of plastic slip is obtained in a system with non-convex energetic hardening through a phenomenological double-well plastic potential. In the current multi-dimensional multi-slip analysis the non-convexity enters the framework through a latent hardening potential presented in Ortiz and Repettto (1999) where the microstructure evolution is obtained explicitly via a lamination procedure. The current study aims the implicit evolution of deformation patterns due to the incorporated physically based non-convex potential.     

An experimental study on grain deformation and interactions in an Al-0.5%Mg multicrystal

Heterogeneous plastic deformation behavior of a coarse-grained Al-0.5%Mg multicrystal was investigated experimentally at the individual grain level. A flat uniaxial tensile specimen consisting of a single layer of millimeter-sized grains was deformed quasi-statically up to an axial strain of 15% at room temperature. The initial local crystallographic orientations of the grains and their evolutions after 5, 12, and 15% plastic strains were measured by electron backscattered diffraction pattern analysis in a scanning electron microscope. The local in-plane plastic strains and rigid body rotations of the grains were measured by correlation of digital optical video images of the specimen surface acquired during the tensile test. It is found that both intergranular and intragranular plastic deformation fields in the aluminum multicrystal specimen under uniaxial tension are highly heterogeneous. Single or double sets of slip-plane traces were predominantly observed on the electro-polished surfaces of the millimeter-sized grains after deformation. The active slip systems associated with these observed slip-plane traces were identified based on the grain orientation after deformation, the Schmid factor, and grain interactions in terms of the slip-plane trace morphology at grain boundaries. It is found that the aluminum multicrystal obeys neither the Sachs nor the Taylor polycrystal deformation models but deforms heterogeneously to favor easy slip transmission and accommodation among the grains.     

Anti-plane shear cracks in ideally plastic crystals

Cracks in ductile single crystals are analyzed here for geometries and orientations such that two-dimensional states of anti-plane shear constitute possible deformation fields. The crystals are modelled as ideally plastic and yield according to critical resolved shear stresses on their slip systems. Restrictions on the asymptotic forms of stress and deformation fields at crack tips are established for anti-plane loading of stationary and quasistatically growing cracks, and solutions are presented for several specific orientations in f.c.c. and b.c.c. crystals. The asymptotic solutions are complemented by complete elastic-plastic solutions for stationary and growing cracks under small scale yielding, based on previous work by Rice (1967, 1984) and Freund (1979). Remarkably, the plastic zone at a stationary crack tip collapses into discrete planes of displacement and stress discontinuity emanating from the tip; plastic flow consists of concentrated shear on the displacement discontinuities. For the growing crack these same planes, if not coincident with the crack plane, constitute collapsed plastic zones in which velocity and plastic strain discontinuities occur, but across which the stresses and anti-plane displacement are fully continuous. The planes of discontinuity are in several cases coincident with crystal slip planes but it is shown that this need not be the case, e.g., for orientations in which anti-plane yielding occurs by multi-slip, or for special orientations in which the crack tip and the discontinuity planes are perpendicular to the activated slip plane.     

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. It is found that the near-tip deformation field, especially, the development of kink or slip shear bands, is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.     

基于边界位移的复合材料力学参数无损反演方法

非均质材料参数识别在工程学、医学以及生物力学等众多领域具有重要意义.目前求解材料参数识别这类反问题主要采用优化方法,通常需要已知结构的全场位移信息,使含有位移的目标函数最小化,从而获得材料参数分布.然而在实际工程中,结构内部的位移较难测量且测量精度低.因此,本文拟提出一类仅利用边界位移就能进行非均质材料参数分布反演的方...     

Computation of Full-field Strains Using Principal Component Analysis

The primary output from several full-field deformation measurement techniques, e.g., Digital Image Correlation (DIC), is the displacement vector at a dense grid of points covering the area of interest. Since such displacement data inherently contain noise, they are usually smoothed first and then differentiated to obtain strains. Another common approach is to use finite-element shape functions for the strains and compute them by treating the measured displacements as nodal displacements. In this paper, we propose a novel method for strain calculation from full-field data, based on the multivariate analysis technique of Principal Component Analysis (PCA) using which we first obtain the singular values and singular vectors for each component of the displacement field. By choosing only the dominant singular values and their corresponding singular vectors, we show that the dimensionality of the displacement data is sharply reduced and a significant portion of the noise is eliminated. Moreover, the shapes of the dominant singular vectors offer physical insight into dominant deformation patterns. We demonstrate the accuracy of the proposed technique by applying it to two cases each of homogeneous and inhomogeneous strain fields and show that in all cases the proposed method yields excellent results.     

Non-local crystal plasticity model with intrinsic SSD and GND effects

A strain gradient-dependent crystal plasticity approach is presented to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. In order to be capable of predicting scale dependence, the heterogeneous deformation-induced evolution and distribution of geometrically necessary dislocations (GNDs) are incorporated into the phenomenological continuum theory of crystal plasticity. Consequently, the resulting boundary value problem accommodates, in addition to the ordinary stress equilibrium condition, a condition which sets the additional nodal degrees of freedom, the edge and screw GND densities, proportional (in a weak sense) to the gradients of crystalline slip. Next to this direct coupling between microstructural dislocation evolutions and macroscopic gradients of plastic slip, another characteristic of the presented crystal plasticity model is the incorporation of the GND-effect, which leads to an essentially different constitutive behaviour than the statistically stored dislocation (SSD) densities. The GNDs, by their geometrical nature of locally similar signs, are expected to influence the plastic flow through a non-local back-stress measure, counteracting the resolved shear stress on the slip systems in the undeformed situation and providing a kinematic hardening contribution. Furthermore, the interactions between both SSD and GND densities are subject to the formation of slip system obstacle densities and accompanying hardening, accountable for slip resistance. As an example problem and without loss of generality, the model is applied to predict the formation of boundary layers and the accompanying size effect of a constrained strip under simple shear deformation, for symmetric double-slip conditions.     

Topology optimization of truss structures with systematic reliability constraints under multiple loading cases

In this paper, a mathematical model for topology optimization of truss structures with constraints of displacement and system reliability under multiple loading cases is constructed. In order to avoid the difficulty of computing the structure's system reliability, a solving approach is presented in which the failure probability of system is divided into the sum of all bars' failure probability by means of reliability distribution. In addition, by drawing into the reliability safety factor and the fundamental relationship in structural mechanics, all probability constraints of displacement and stress are equivalently displayed as conventional form and linear function of the design variables. The optimization problem with multiple constraints is treated by the compact constraint tactics and is solved by the improved simplex method. The examples show that the approach proposed in this paper is feasible and efficient. The project supported by the National Natural Science Foundation of China.     

A semi-implicit integration scheme for rate independent finite crystal plasticity

An efficient new numerical integration scheme is presented for rate independent crystal plasticity theory. A key feature of this approach is the ability to identify active slip systems prior to determining their shearing rate. Options are described for various cases of slip system hardening, including self hardening and latent hardening. Alternatives for the constitutive update are explored, including hyperelasticity based on the multiplicative decomposition of the deformation gradient as well as application of the consistency condition in a much more efficient hypoelastic formulation. Several conclusions are drawn concerning the influences of elastic and plastic properties on the activation of slip systems and their subsequent shearing rates. Key among these is the fact that, once activated, shearing rates are independent of the levels of shear flow resistance on the slip systems, provided that the plastic hardening moduli are much less in magnitude than the elastic moduli, as is usually the case. Determination of active slip systems and their shearing rates depend on the degree of elastic anisotropy of the crystal, but not on the magnitude of elastic stiffness.     

Determination of the midsurface of a deformed shell from prescribed surface strains and bendings via the polar decomposition

We show how to determine the midsurface of a deformed thin shell from known geometry of the undeformed midsurface as well as the surface strains and bendings. The latter two fields are assumed to have been found independently and beforehand by solving the so-called intrinsic field equations of the non-linear theory of thin shells. By the polar decomposition theorem the midsurface deformation gradient is represented as composition of the surface stretch and 3D finite rotation fields. Right and left polar decomposition theorems are discussed. For each decomposition the problem is solved in three steps: (a) the stretch field is found by pure algebra, (b) the rotation field is obtained by solving a system of first-order PDEs, and (c) position of the deformed midsurface follows then by quadratures. The integrability conditions for the rotation field are proved to be equivalent to the compatibility conditions of the non-linear theory of thin shells. Along any path on the undeformed shell midsurface the system of PDEs for the rotation field reduces to the system of linear tensor ODEs identical to the one that describes spherical motion of a rigid body about a fixed point. This allows one to use analytical and numerical methods developed in analytical mechanics that in special cases may lead to closed-form solutions.     

解决约束违背问题的一种自适应调整方法

带约束的优化问题的目的是要找到满足等式或者不等式约束的最优点。在某些情况下,优化求解得到的"最优点"可能会使得某个或某几个约束条件超出目标约束限,或者在所有约束条件中的最大值远远小于目标约束限。针对这一类问题,本文提出一种在寻优过程的每一次迭代中自适应调整约束限的方法,通过动态调整迭代过程中迭代模型约束限的值,将约束条件中最大值的约束条件变为等式约束,使得迭代解始终在可行域范围内,且收敛后的最优解不违背任何约束条件。本文将该方法成功应用于位移约束下结构重量最小化拓扑优化模型,原来不满足约束条件的情况在使用该方法后都能使约束得到满足,解决了约束条件被违背的问题。     

Dislocation pile-ups in bicrystals within continuum dislocation theory

Within continuum dislocation theory the plastic deformation of bicrystals under a mixed deformation of plane constrained uniaxial extension and shear is investigated with regard to the nucleation of dislocations and the dislocation pile-up near the phase boundaries of a model bicrystal with one active slip system within each single crystal. For plane uniaxial extension, we present a closed-form analytical solution for the evolution of the plastic distortion and of the dislocation network in the case of symmetric slip planes (i.e. for twins), which exhibits an energetic as well as a dissipative threshold for the dislocation nucleation. The general solution for non-symmetric slip systems is obtained numerically. For a combined deformation of extension and shear, we analyze the possibility of linearly superposing results obtained for both loading cases independently. All solutions presented in this paper also display the Bauschinger effect of translational work hardening and a size effect typical to problems of crystal plasticity.     

Comparisons of crystal hardening laws in multiple slip

This paper brings together and concisely reviews results from recent analytical investigations on single crystals (variously done alone or with students) in which predictions from different theoretical hardening laws are contrasted and compared with experimental studies. Finitely deforming f.c.c. crystals in both constrained and unconstrained multiple-slip configurations are considered. Four crystal hardening laws are given prominence. Two of these belong to a class of theories in which the physical hardening moduli relating rates-of-change of critical strengths (in the 24 crystallographically equivalent slip systems) to slip-rates are taken as symmetric. These are G.I. Taylor's classic isotropic hardening rule (proposed in 1923), which is almost universally adopted in the metallurgical literature for various approximate analyses of single and poly-crystal deformation, and a 2-parameter modification of Taylor's rule that has an empirical basis in the qualitative features of experimentally determined latent hardening in single slip. The other two hardening laws featured here belong to a class of theories that were introduced in 1977 by this author. This class requires the above moduli to be nonsymmetric and explicity dependent upon the current stress state in such a manner that the following consequences are assured. (1) The deformation-dependent hardening of latent slip systems necessarily develops anisotropically if there is relative rotation of gross material and underlying crystal lattice. (2) The theories admit self-adjoint boundary value problems for crystalline aggregates, hence a variational formulation. (The fact that symmetric physical hardening moduli do not permit variational formulations of polycrystalline problems was shown at the 1972 Warsaw Symposium.) The two members of this class considered here are the original (and simplest p possible) theory of rotation-dependent anisotropy, which was proposed by this author in 1977 and commonly has been referred to as the “simple theory,” and a modification of this theory introduced in 1982 by Peirce, Asaro and Needleman that lies between Taylor's rule and the simple theory in its predictions for finitely deforming f.c.c. crystals. (In a series of five papers during 1977–1979, the simple theory was shown to universally account for the experimental phenomenon of “overshooting” in single slip in both f.c.c. and b.c.c. crystals.) Theoretical results from the various hardening rules are contrasted and compared with finite strain experiments in the metallurgical literature. Both tensile-loaded crystals in 4, 6 and 8-fold symmetry orientations and compressively loaded crystals under conditions of channel die constraint are treated. A postulate of minimum plastic work introduced in 1981 plays a prominent role in the theoretical analyses, in many cases providing a unique solution to the slip system inequalities and deformation constraints (where applicable). The rather remarkable ability of the simple theory to reconcile diverse qualitative features of both constrained and unconstrained finited deformation of f.c.c. crystals is demonstrated. Finally, conditions for total loading (all systems active) in 6-fold symmetry are investigated, and certain concepts regarding the selection of active systems under prescribed straining are critically assessed.     

Cylindrical void in a rigid-ideally plastic single crystal III: Hexagonal close-packed crystal

The analytical solution is derived for the plane strain stress field around a cylindrical void in a hexagonal close-packed single crystal with three in-plane slip systems oriented at the angle π/3 with respect to one another. The critical resolved shear stress on each slip system is assumed to be equal. The crystal is loaded by both internal pressure and a far-field equibiaxial compressive stress. The deformation field takes the form of angular sectors, called slip sectors, within which only one slip system is active; the boundaries between different sectors are radial lines. The stress fields are derived by enforcing equilibrium and a rigid, ideally plastic constitutive relationship, in the spirit of anisotropic slip line theory. The results show that each slip sector is divided into smaller regions denoted as stress sectors and the stress state valid within each stress sector is derived. It is shown that stresses are unique and are continuous within stress sectors and across stress sector boundaries, but the gradient of stresses is not continuous across the boundaries between stress sectors. The solution shows self-similarity in that the stresses over the entire domain can be determined from the stresses within a small region adjacent to the void by invoking certain scaling and symmetry properties. In addition, the stress state exhibits periodicity along logarithmic spirals which emanate from the void. The results predict that the mean value of in-plane pressure required to activate plastic deformation around a void in a single crystal can be higher than that necessary for a void in an isotropic material and is sensitive to the orientation of the slip systems relative to the void.     

Experimental analysis of crack tip fields in rubber materials under large deformation

A three-nested-deformation model is proposed to describe crack-tip fields in rubber-like materials with large deformation.The model is inspired by the distribution of the measured in-plane and out-of-plane deformation.The inplane displacement of crack-tip fields under both Mode I and mixed-mode(Mode I-II) fracture conditions is measured by using the digital Moire’ method.The deformation characteristics and experimental sector division mode are investigated by comparing the measured displacement fields under different fracture modes.The out-of-plane displacement field near the crack tip is measured using the three-dimensional digital speckle correlation method.     

“Finite-Element” Displacement Fields Analysis from Digital Images: Application to Portevin–Le Châtelier Bands

A new methodology is proposed to estimate displacement fields from pairs of images (reference and strained) that evaluates continuous displacement fields. This approach is specialized to a finite-element decomposition, therefore providing a natural interface with a numerical modeling of the mechanical behavior used for identification purposes. The method is illustrated with the analysis of Portevin–Le Chatelier bands in an aluminum alloy sample subjected to a tensile test. A significant progress with respect to classical digital image correlation techniques is observed in terms of spatial resolution and uncertainty.     

Experimental and numerical investigations of the deformation of soft materials under tangential loading

The surface ‘tensile test’, in which tangential loads are applied through surface mounted adhesive tapes, is a viable method for the assessment of mechanical properties of soft materials, particularly biological soft materials in vivo. In the present work the deformation pattern and force–displacement relationship in the surface tensile test were experimentally investigated using surface displacement analysis (SDA) and numerically simulated using finite element modelling. The experimental and FE results showed close agreement using silicone rubber as a model material. The force–displacement relationship was found to be dependent on the tape separations. SDA measurements and FE simulation showed that the displacement and strain fields were not uniform and the distribution pattern varies with tape separation. A combined experimental–numerical approach to inversely extract material properties using multiple tests with different length scales is proposed and assessed using a model material.