An interfacial energy incorporated couple stress crystal plasticity and the finite element simulation of grain subdivision

A couple stress crystal plasticity formulation that incorporates interfacial couple stress energy was proposed in terms of the virtual work-rate principle for finite element method. By applying the assumed constitutive models of couple stress at the grain boundary as well as the grain interior, finite element simulations were conducted for various crystal models, with different grain subdivision models to examine how plastic deformation work is affected by grain subdivision from the interfacial couple stress energy effect.Finite element simulation results showed that the amount of predicted plastic deformation work depends on grain subdivision, and that the amount of work can be minimized for a particular grain subdivision. We inferred from the simulation results that actual grain subdivision might correspond to the minimum amount of plastic deformation work and, if this correlation is validated, actual grain subdivision might be predicted based on the interfacial energy incorporated couple stress crystal plasticity.     

Finite element modeling of crystal plasticity with grains shaped as truncated octahedrons

Different modeling strategies are tested for the prediction of texture development and microscopic strain heterogeneity in cold-rolled ULC steel and in multiphase steel under uniaxial tension. The polycrystalline aggregate is represented by a finite element mesh that is loaded under periodic boundary conditions. Grains are shaped as cubes or as truncated octahedrons, defining three levels of mesh refinement. Simulations rely on a simplified implementation of crystal plasticity, in which elastic strains are considered infinitesimal. The constitutive law is integrated fully implicitly in a reference frame bounded to the crystal lattice. Accuracy of the time-integration procedure is assessed by referring to some predictions of the crystal plasticity algorithm developed by Kalidindi et al. [Crystallographic texture evolution in bulk deformation processing of FCC metals. J. Mech. Phys. Solid 40 (1992) 537–569]. Then, results of the micro–macro modeling are compared to experimental data. It is found that the simulations with truncated octahedral grains yield improved predictions compared to those with cuboidal grains.     

粘弹性大变形动力响应的有限元分析

本文基于ＴｏｔａｌＬａｇｒａｎｇｉａｎ增量叠加方法，采用Ｋｉｒｃｈｈｏｆｆ应力增量和Ｇｒｅｅｎ应变增量表示的动力虚功方程和Ｋｉｒｃｈｈｏｆｆ应力－Ｇｒｅｅｎ应变的单积分型本构关系，导出粘弹性大变形的动力变分方程。依此采用Ｎｅｗｍａｒｋ法和八节点轴对称等参数元与二十节点三维等参数元编制了轴对称及三维问题的动力响应计算程序，典型例题的计算结果表明分析符合结构的物理性质。     

闭孔泡沫金属变形模式的有限元分析

运用有限元软件ABAQUS/Explicit模拟了三维Voronoi闭孔泡沫金属在不同的冲击速度下的变形行为。随着冲击速度的提高,得到了3种变形模式:准静态均匀模式、过渡模式和冲击模式,并以相对密度和冲击速度为坐标建立了变形模式图。引入应力均匀性指标和变形局部化指标,确定了模式转化的临界速度,并与已有的冲击速度预测公式进行了比较。根据临界速度的数值和理论结果,提出了一种确定锁定应变的方案,结果介于压实应变和完全密实应变之间。     

Assumed stress quasi-conforming triangular element for couple stress theory

In this paper, a 3-node triangular element for couple stress theory is proposed based on the assumed stress quasi-conforming method. The formulation starts from polynomial approximation of stresses. Then the stress-function matrix is treated as the weighted function to weaken the strain-displacement equations. Finally, the string-net functions are introduced to calculate strain integration and the stress smooth technique is adopted to improve the stress accuracy. Numerical results show that the proposed new model can pass the C~(0-1) patch test with excellent precision, does not exhibit extra zero energy modes and can capture the scale effects of microstructure.     

Finite element simulation of localized plastic deformation

Summary A rational finite element algorithm for capturing localized plastic deformation due to softening and/or non-associated plastic flow is devised. The incremental relations are based on implicit integration. In each increment it is important to use a carefully designed modified Newton iteration procedure in conjunction with a start solution whose calculation is based on a diagnostic bifurcation analysis. The algorithm performs successfully for a finite element mesh that is biased in the sense that the element sides are prealigned with the anticipated localization zone which is demonstrated for a slope stability problem.

---

Finite-Element-Berechnung von örtlichen plastischen Verformungen

Übersicht Die Konstruktion eines Algorithmus, der im Zusammenhang mit der Methode der Finiten Elemente örtliche plastische Verformungen als Folge von Materialentfestigung und/oder nichtassoziiertem plastischen Fließen erfaßt, wird beschrieben. Die inkrementellen Gleichungen ergeben sich aus der Verwendung eines impliziten Integrationsverfahrens. In jedem Inkrement wird ein sorgfältig entwickeltes, modifiziertes Newton-Iterationsverfahren verwendet, dessen erste Lösungsabschätzung auf einer diagnostischen Analyse des Verzweigungsproblems beruht. Am Beispiel einer Böschungsstabilitätsungtersuchung wird gezeigt, daß der Algorithmus erfolgreich ist, wenn das verwendete Finite-Elemente-Netz so ausgerichtet ist, daß die Elementseiten der erwarteten Versagenzone folgen.

---

---

Presented at the workshop on Numerical Methods for Localization and Bifurcation of Granular Bodies, held at the Technical University of Gdansk (Poland), September 25–30, 1989     

Finite element analysis of thermo-elastic behavior of piezoelectric structures under finite deformation

It is noted that the behavior of most piezoelectric materials is temperature dependent and such piezo-thermo-elastic coupling phenomenon has become even more pronounced in the case of finite deformation. On the other hand, for the purpose of precise shape and vibration control of piezoelectric smart structures, their deformation under external excitation must be ideally modeled. This demands a thorough study of the coupled piezo-thermo-elastic response under finite deformation. In this study, the governing equations of piezoelectric structures are formulated through the theory of virtual displacement principle and a finite element method is developed. It should be emphasized that in the finite element method the fully coupled piezo-thermo-elastic behavior and the geometric non-linearity are considered. The method developed is then applied to simulate the dynamic and steady response of a clamped plate to heat flux acting on one side of the plate to mimic the behavior of a battery plate of satelite irradiated under the sun. The results obtained are compared against classical solutions, whereby the thermal conductivity is assumed to be independent of deformation. It is found that the full-coupled theory predicts less transient response of the temperature compared to the classic analysis. In the steady state limit, the predicted temperature distribution within the plate for small heat flux is almost the same for both analyses. However, it is noted that increasing the heat flux will increase the deviation between the predictions of the temperature distribution by the full coupled theory and by the classic analysis. It is concluded from the present study that, in order to precisely predict the deformation of smart structures, the piezo-thermo-elastic coupling, geometric non-linearity and the deformation dependent thermal conductivity should be taken into account. Project supported by the National Natural Science Foundation of China (Nos. 10132010 and 50135030) and the Foundation of In-service Doctors of Xi'an Jiaotong University.     

A finite deformation theory of higher-order gradient crystal plasticity

For higher-order gradient crystal plasticity, a finite deformation formulation is presented. The theory does not deviate much from the conventional crystal plasticity theory. Only a back stress effect and additional differential equations for evolution of the geometrically necessary dislocation (GND) densities supplement the conventional theory within a non-work-conjugate framework in which there is no need to introduce higher-order microscopic stresses that would be work-conjugate to slip rate gradients. We discuss its connection to a work-conjugate type of finite deformation gradient crystal plasticity that is based on an assumption of the existence of higher-order stresses. Furthermore, a boundary-value problem for simple shear of a constrained thin strip is studied numerically, and some characteristic features of finite deformation are demonstrated through a comparison to a solution for the small deformation theory. As in a previous formulation for small deformation, the present formulation applies to the context of multiple and three-dimensional slip deformations.     

有限元表面应力计算

用有限元[1]通用程序进行结构计算时,最常用的是位移法,因而计算得到的位移有较高的精度。由位移计算应力时,有限元法应用的是应力-应变关系和应变-位移关系,其中应变-位移是微商关系。在数值计算中,微商只能转化为差商等用插值近似处理。这样,虽然位移精度高,但应力的计算精度就被大打折扣。本文应用弹性力学辛体系理论[2],解析求解了位移和应力的影响函数。利用有限元程序计算得到的位移,由功互等定理,不需要微分插值,就可以得到指定点的应力,应力精度大大提高。工程实际中有许多问题的最大应力往往发生在构件表面。针对表面应力问题,本文给出了半平面表面应力的影响函数,进行了数值算例计算。计算结果表明,用本文提出的影响函数法求解一点的应力,其精度明显提高,并且计算结果有很好的稳定性。用本文的影响函数法编制成子程序,可作为有限元软件应力计算的一个模块,可以更好地发挥有限元程序的功效。     

Finite deformation analysis of mechanism-based strain gradient plasticity: torsion and crack tip field

A finite deformation theory of mechanism-based strain gradient (MSG) plasticity is developed in this paper based on the Taylor dislocation model. The theory ensures the proper decomposition of deformation in order to exclude the volumetric deformation from the strain gradient tensor since the latter represents the density of geometrically necessary dislocations. The solution for a thin cylinder under large torsion is obtained. The numerical method is used to investigate the finite deformation crack tip field in MSG plasticity. It is established that the stress level around a crack tip in MSG plasticity is significantly higher than its counterpart (i.e. HRR field) in classical plasticity.     

Micromorphic continuum. Part II: Finite deformation plasticity coupled with damage

It is demonstrated how a micromorphic plasticity theory may be formulated on the basis of multiplicative decompositions of the macro- and microdeformation gradient tensor, respectively. The theory exhibits non-linear isotropic and non-linear kinematic hardening. The yield function is expressed in terms of Mandel stress and double stress tensors, appropriately defined for micromorphic continua. Flow rules are derived from the postulate of Il’iushin and represent generalized normality conditions. Evolution equations for isotropic and kinematic hardening are introduced as sufficient conditions for the validity of the second law of thermodynamics in every admissible process. Finally, it is sketched how isotropic damage effects may be incorporated in the theory. This is done for the concept of effective stress combined with the hypothesis of strain equivalence.     

On the partition of rate of deformation in crystal plasticity

Two different partitions of the rate of deformation tensor into its elastic and plastic parts are derived for elastic–plastic crystals in which crystallographic slip is the only cause of plastic deformation. One partition is associated with the Jaumann, and the other with convected rate of the Kirchhoff stress. Different expressions for the plastic part of the rate of deformation are obtained, and corresponding constitutive inequalities discussed. Relationship with the plastic part of the rate of the Lagrangian strain is also given.     

单晶体和双晶体微观层次变形行为的有限元分析

从微观层次上研究金属材料的变形行为,将位错引入到本构关系中,用硬化函数描述材料的硬化规律,考虑了变形的率相关性,采用三维模型用大变形有限单元法对单晶体在单向拉伸载荷和循环载荷作用下的变形行为、双晶体在单向拉伸作用下滑移系的开动进行了模拟计算,得到了与实验一致的计算结果。     

Using texture components in crystal plasticity finite element simulations

We discuss methods to map crystallographic textures in crystal plasticity finite element simulations. Fourier-type series expansion methods which use spherical harmonic library functions as well as the direct pole figure inversion methods are not well suited to reproduce texture information in a sufficiently localized spherical form onto finite element grids. Mathematically compact Gauss-shaped spherical texture components represent a better approach for including textures in finite element models since they represent an excellent compromise between discreteness (spherical localization), compactness (simple functions), mathematical precision (very good approximation also of complex orientation distribution functions already with small sets of texture components), scalability (the number of used texture components can be systematically varied according to the desired exactness of the texture fit), conceptual simplicity (simple mathematical handling), and physical significance (texture components can be directly linked to characteristic metallurgical mechanisms). The use of texture component functions has also advantages over the use of large sets of discrete single orientations with equal scatter and height since they are more compact, practical, and provide better physical insight into microstructural mechanisms and composition sensitive effects. The article presents a new approach for the mathematical reproduction of such crystallographic texture components in crystal plasticity finite element simulations. It explains in some detail why they are particularly suited for this purpose and how they can be used to map and recover textures in/from plasticity simulations.     

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.     

On large-strain finite element solutions of higher-order gradient crystal plasticity

A finite-strain higher-order gradient crystal plasticity model accounting for the backstress effect originating from the existence of geometrically necessary dislocations (GNDs) is applied to plane strain finite element analysis. Different element types are tested to seek out an element formulation that is reliable and useful for solving problems involving severe plastic deformation. In the present finite element formulation, the GND density rates are chosen to be additional nodal degrees of freedom. Different orders of shape functions are employed for the interpolation of displacement rates and GND density rates. Their effects on solutions are examined in detail by considering three boundary value problems: a simple shear of a constrained layer (a film), a compression problem with loading surfaces impenetrable to dislocations, and a tension problem involving shear band formation. In all the cases, the formulation in which eight-node elements with reduced integration and four-node elements with full integration are used respectively for displacement rates and the GND density rates gives reasonable solutions. In addition to the discussion on the choice of finite elements, detailed behavior in gradient-dependent solids, such as the accumulation of GND density and the distribution of backstress on each slip system, is investigated by utilizing the reliable computational results obtained.     

Finite element method for the stress analysis of isotropic cylindrical helical spring

This paper presents an efficient two nodes finite element with six degrees of freedom per node, capable to model the total behaviour of a helical spring. The formulation, which includes the shear deformation effects, is based on the assumed forces hybrid approach. The resultant forces approximation verifies exactly the resultant equilibrium equations. The developed model proves its accuracy compared with other elements. This element permits to get the distribution of different stresses along the spring and through the wire surface by only one element.