Associative coupled thermoplasticity at finite strain with temperature-dependent material parameters

The paper deals with associative coupled thermoplasticity at finite strains. J2 plasticity model is based on hyperelastic formulation with multiplicative decomposition of deformation gradient into elastic and plastic parts. As a novel aspect, temperature dependence of all material parameters is introduced. For such model, variational procedure is carried out and discretization by the mixed finite element method is performed. Consistent linearization is carried out and algorithmic elasto-plastic tangent moduli are given. Verification of algorithm is provided by two examples.     

Finite element implementation of a generalised non-local rate-dependent crystallographic formulation for finite strains

This work describes the finite element implementation of a generalised strain gradient and rate-dependent crystallographic formulation for finite strains and general anisothermal conditions based on a multiplicative decomposition of the deformation gradient. The implementation involved the development of both a novel finite element formulation to determine the spatial slip rate gradients at each material point, and an implicit numerical integration scheme at the constitutive level to update the stresses and solution dependent variables. The time-integration procedure uses a Newton–Raphson scheme with a single level of iteration to solve the incremental non-linear equations associated with the non-local constitutive formulation. Closed-form solutions for the relevant fourth-order Jacobian tensors are given. The proposed numerical scheme is formulated in a general form and hence should be applicable to most existing crystallographic models. The crystallographic formulation is then used to investigate the effect of the morphology and volume fraction of the reinforcing phase of a two-phase single crystal on its macroscopic behaviour.     

A fully coupled diffusional-mechanical formulation: numerical implementation,analytical validation,and effects of plasticity on equilibrium

A macroscopic coupled stress-diffusion theory which accounts for the effects of nonlinear material behaviour, based on the framework proposed by Cahn and Larché, is presented and implemented numerically into the finite element method. The numerical implementation is validated against analytical solutions for different boundary valued problems. Particular attention is payed to the open system elastic constants, i.e. those derived at constant diffusion potential, since they enable, under circumstances, the equilibrium composition field for any generic chemical-mechanical coupled problem to be obtained through the solution of an equivalent elastic problem. Finally, the effects of plasticity on the overall equilibrium state of the coupled problem solution are discussed.     

A formulation of anisotropic continuum elastoplasticity at finite strains. Part I: Modelling

A constitutive model for anisotropic elastoplasticity at finite strains is developed together with its numerical implementation. An anisotropic elastic constitutive law is described in an invariant setting by use of structural tensors and the elastic strain measure C_e. The elastic strain tensor as well as the structural tensors are assumed to be invariant in relation to superimposed rigid body rotations. An anisotropic Hill-type yield criterion, described by a non-symmetric Eshelby-like stress tensor and further structural tensors, is developed, where use is made of representation theorems for functions with non-symmetric arguments. The model also considers non-linear isotropic hardening. Explicit results for the specific case of orthotropic anisotropy are given. The associative flow rule is employed and the features of the inelastic flow rule are discussed in full. It is shown that the classical definition of the plastic material spin is meaningless in conjunction with the present formulation. Instead, the study motivates an alternative definition, which is based on the demand that such a quantity must be dissipation-free, as the plastic material spin is in the case of isotropy. Equivalent spatial formulations are presented too. The full numerical treatment is considered in Part II.     

On a formulation for anisotropic elastoplasticity at finite strains invariant with respect to the intermediate configuration

The paper is concerned with a formulation of anisotropic finite strain inelasticity based on the multiplicative decomposition of the deformation gradient F=F_eF_p. A major feature of the theory is its invariance with respect to rotations superimposed on the inelastic part of the deformation gradient. The paper motivates and shows how such an invariance can be achieved. At the heart of the formulation is the mixed-variant transformation of the structural tensor, defined as the tensor product of the privileged directions of the material as given in a reference configuration, under the action of F_p. Issues related to the plastic material spin are discussed in detail. It is shown that, in contrast to the isotropic case, any flow function formulated purely in terms of stress quantities, necessarily exhibits a non-vanishing plastic material spin. The possible construction of spin-free rates is discussed as well, where it is shown that the flow rule must then depend not only on the stress but on the strain as well.     

Shakedown reduced kinematic formulation,separated collapse modes,and numerical implementation

Our shakedown reduced kinematic formulation is developed to solve some typical plane stress problems, using finite element method. Whenever the comparisons are available, our results agree with the available ones in the literature. The advantage of our approach is its simplicity, computational effectiveness, and the separation of collapse modes for possible different treatments. Second-order cone programming developed for kinematic plastic limit analysis is effectively implemented to study the incremental plasticity collapse mode. The approach is ready to be used to solve general shakedown problems, including those for elastic–plastic kinematic hardening materials and under dynamic loading.     

基于内聚力理论的二维二次界面单元在ABAQUS中的UEL程序实现

依据有限元理论,结合内聚力模型法则,推导出二维二次粘结界面单元在大位移情况下的数值格式,得到用形函数表示的单元位移模式、载荷向量和刚度矩阵,并进行了离散化。基于ABAQUS软件的自定义扩展模块,编制了相应的用户单元子程序UEL,通过数值算例验证了该程序的准确性和有效性。这一成果能为在ABAQUS软件中开展相关数值研究,以及开发其他类型的内聚力界面有限单元提供思路和参考。     

Framework for non-coherent interface models at finite displacement jumps and finite strains

This paper deals with a novel constitutive framework suitable for non-coherent interfaces, such as cracks, undergoing large deformations in a geometrically exact setting. For this type of interface, the displacement field shows a jump across the interface. Within the engineering community, so-called cohesive zone models are frequently applied in order to describe non-coherent interfaces. However, for existing models to comply with the restrictions imposed by (a) thermodynamical consistency (e.g., the second law of thermodynamics), (b) balance equations (in particular, balance of angular momentum) and (c) material frame indifference, these models are essentially fiber models, i.e. models where the traction vector is collinear with the displacement jump. This constraints the ability to model shear and, in addition, anisotropic effects are excluded. A novel, extended constitutive framework which is consistent with the above mentioned fundamental physical principles is elaborated in this paper. In addition to the classical tractions associated with a cohesive zone model, the main idea is to consider additional tractions related to membrane-like forces and out-of-plane shear forces acting within the interface. For zero displacement jump, i.e. coherent interfaces, this framework degenerates to existing formulations presented in the literature. For hyperelasticity, the Helmholtz energy of the proposed novel framework depends on the displacement jump as well as on the tangent vectors of the interface with respect to the current configuration – or equivalently – the Helmholtz energy depends on the displacement jump and the surface deformation gradient. It turns out that by defining the Helmholtz energy in terms of the invariants of these variables, all above-mentioned fundamental physical principles are automatically fulfilled. Extensions of the novel framework necessary for material degradation (damage) and plasticity are also covered.     

Thermomechanical modeling of metals at finite strains: First and mixed order finite elements

The aim of this paper is to describe an updated EAS (Enhanced Assumed Strain) finite element formalism developed to model the thermomechanical behavior of metals submitted to large strains. We will also expose the use of mixed order elements (first order mechanical elements strongly coupled with quadratic thermal elements) which, as we will show, is of particular interest for modeling fast processes inducing important temperature gradients. The features of this formalism, used jointly with an Updated Lagrangian approach and an hypoelastic anisothermal constitutive formulation, will be described. Three applications involving finite strains and important thermomechanical couplings will be studied. The results obtained will be compared with the results given by the now classical SRI (Selective Reduced Integration) formalism.     

A three dimensional field formulation,and isogeometric solutions to point and line defects using Toupin's theory of gradient elasticity at finite strains

We present a field formulation for defects that draws from the classical representation of the cores as force dipoles. We write these dipoles as singular distributions. Exploiting the key insight that the variational setting is the only appropriate one for the theory of distributions, we arrive at universally applicable weak forms for defects in nonlinear elasticity. Remarkably, the standard, Galerkin finite element method yields numerical solutions for the elastic fields of defects that, when parameterized suitably, match very well with classical, linearized elasticity solutions. The true potential of our approach, however, lies in its easy extension to generate solutions to elastic fields of defects in the regime of nonlinear elasticity, and even more notably for Toupin's theory of gradient elasticity at finite strains (Toupin Arch. Ration. Mech. Anal., 11 (1962) 385). In computing these solutions we adopt recent numerical work on an isogeometric analytic framework that enabled the first three-dimensional solutions to general boundary value problems of Toupin's theory (Rudraraju et al. Comput. Methods Appl. Mech. Eng., 278 (2014) 705). We first present exhaustive solutions to point defects, edge and screw dislocations, and a study on the energetics of interacting dislocations. Then, to demonstrate the generality and potential of our treatment, we apply it to other complex dislocation configurations, including loops and low-angle grain boundaries.     

Improved Carroll's hyperelastic model considering compressibility and its finite element implementation

Acta Mechanica Sinica - In engineering component design, material models are increasingly used in finite element simulations for an expeditious and less costly analysis of the design prototypes. As...     

On strain localization analysis of elastoplastic materials at finite strains

The problem of strain localization into planar bands of rate-independent elastoplastic solids with smooth yield surface and plastic potential is analyzed reconsidering the work of Rice and Rudnicki in 1980. It is shown that strain localization with elastic unloading on one side of the band first becomes possible either at localization in the comparison solid corresponding to the loading branch of the constitutive equation or at the snap-back threshold. The elastic unloading is shown to start from the condition of neutral loading, occurring in fact at the onset of localization. The case of localization with elastic unloading into the band and plastic loading outside that was not considered by Rice and Rudnicki is taken into account.     

An unstructured-grid numerical model for interfacial multiphase fluids based on multi-moment finite volume formulation and THINC method

We present a practical numerical framework for incompressible interfacial multiphase flows on unstructured grids with arbitrary and hybrid elements. The numerical framework is constructed by combining VPM (volume-average/point-value multi-moment) and UMTHINC (unstructured multi-dimensional tangent of hyperbola interface capturing) schemes. To facilitate accurate and reliable simulations for interfacial multiphase flows on arbitrary and hybrid unstructured grids, we have made the following major new efforts in this work. (1) UMTHINC scheme on prismatic and pyramidal elements to facilitate computations on hybrid arbitrary unstructured grids; (2) Consistent numerical formulation for mass and momentum transports to simulate multiphase flows of large density ratio; (3) Combined FVM-FEM for accurate solution to diffusion equation; (4) Pressure-projection formulation in consistent with the balanced-force model. Integrating all these numerical techniques effectively enhances the accuracy and robustness in interface capturing and numerical solution of multiphase fluid dynamics, which results in a numerical framework of great significance for practical applications. Numerical verifications have been carried out through benchmark tests ranging from surface tension dominant flows of small scale to large scale flows with violently-changing interfaces. Numerical results demonstrate that the present framework is robust with adequate accuracy for simulating multiphase flows in complex geometries.     

Closed set of the uniqueness conditions and bifurcation criteria in generalized coupled thermoplasticity for small deformations

This paper reports the results of a study into global and local conditions of uniqueness and the criteria excluding the possibility of bifurcation of the equilibrium state for small strains. The conditions and criteria are derived on the basis of an analysis of the problem of uniqueness of a solution involving the basic incremental boundary problem of coupled generalized thermo-elasto-plasticity. This work forms a follow-up of previous research (?loderbach in Bifurcations criteria for equilibrium states in generalized thermoplasticity, IFTR Reports, 1980, Arch Mech 3(35):337–349, 351–367, 1983), but contains a new derivation of global and local criteria excluding a possibility of bifurcation of an equilibrium state regarding a comparison body dependent on the admissible fields of stress rate. The thermal elasto-plastic coupling effects, non-associated laws of plastic flow and influence of plastic strains on thermoplastic properties of a body were taken into account in this work. Thus, the mathematical problem considered here is not a self-conjugated problem.     

A thermodynamics framework and numerical aspects for transformation-induced plasticity at large strains

In this work, we present a macroscopic material model for simulation transformation-induced plasticity, which is an important phenomenon in metal forming processes. The model is formulated within a thermodynamic framework at large strains. In order to account for both, phase transformation and plasticity, yield functions are related to these effects. Then, applying the concept of maximum dissipation evolution equations are obtained for the inelastic strains, the transformation strains, a hardening variable and the volume fraction of martensite. Furthermore the numerical implementation of the constitutive equations into a finite element program is described. In a numerical example we investigate the austenite-to-martensite phase transformation in a shaft subjected to thermo-mechanical loading in a hybrid-forming process.     

An incremental theory of large strain and large displacement problems and its finite element formulation

Summary This paper is concerned with the development of an incremental finite'element theory for the large strain and the large displacement problems, referred to the current configuration of the body. Using the convected coordinate system which is embedded in the body, the incremental representations of strain and stress tensors and the energy relations are presented, and then the general procedure to construct the so-called element stiffness matrix in incremental form is considered. The finite element formulation is developed for a typical constitutive relation and it is shown that some correction matrices, some of which have been omitted in the previous works, are to be added to the element stiffness matrix. Finally the method to assemble the equations of the element to the global system is discussed and a simple finite element model satisfying the compatibility condition is presented.The finite element theory developed in this paper is able to be extended to the problems for the general thermodynamical process of a broad class of nonlinear materials.

---

Übersicht Mit Hilfe der Methode der finiten Elemente wird eineZuwachstheorie zurBehandlung von Problemen mit endlicher Verformung abgeleitet. Dabei wird ein im Körper eingebettetes, der momentanen Form angepaßtes Bezugssystem verwendet. Es werden Ausdrücke für die Energie sowie für die Änderungen der Spannungs- und Verformungs-Tensoren abgeleitet und es wird ein Verfahren zur Konstruktion der Steifigkeitsmatrix für ein Element angegeben. Ein typisches Stoffgesetz wird dabei zugrundegelegt. Dabei zeigt es sich, daß einige in früheren Arbeiten vernachlässigte Korrektur-Matrizen zu der Steifigkeits-Matrix des Elementes hinzugefügt werden müssen. Die Möglichkeiten der Zusammenfassung der für die Elemente geltenden Gleichungen zu einem globalen Gleichungssystem werden diskutiert und es wird ein den Verträglichkeitsbedingungen genügendes Elemente-Modell angegeben.Das angegebene Verfahren kann für allgemeine thermodynamische Prozesse in einer breiten Klasse nichtlinearer Materialien erweitert werden.