Modelling Metals At Finite Deformation

During metal forming processes, substantial microstructural changes occur in the material due to large plastic deformations leading to different mechanical properties. It is of great interest to predict the behaviour of these materials at different fabriction stages and of the final product. At first glance, the behaviour of metals can be approached by an elastoplastic isotropic material model with a volumetric-deviatoric split and isotropic hardening. In order to perform the calculations, a logarithmic strain is considered in the principal directions of stress and strain space, allowing to make predictions even at finite deformations. Because of the actual nature of metals, the crystalline structure, the deformation at the microstructural level is much more complex. Due to the mathematically algorithmic form of an elastic predictor and a plastic corrector, the elastoplastic model can be extended to crystal plasticity which is similarly handled in terms of a critical resolved shear stress on defined slip planes in the crystal. Hardening can be modelled through a viscoplastic power law. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation model for anisotropic fibrous materials

Paper and paper–based materials such as cardboard are used in a wide variety of applications and in the development of new applications such as boxes an accurate simulation model is of major importance. Industrially made paper material typically has an orthotropic fibrous structure, due to the manufacturing process, where the fibers tend to align in the direction of motion in the machine. In this work a plasticity–based material model allowing for finite strains is developed. The model is suitable for materials with an anisotropic fibrous structure such as paper. The general framework is based on separate mappings describing the deformations of the continuum and the substructure and a multiplicative split of these mappings into elastic and plastic parts. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Deformation of composites with arbitrarily oriented orthotropic fibers under matrix microdamages

In the present work, a model of nonlinear deformation of stochastic composites under microdamaging is developed for the case of a composite with orthotropic inclusions, when microdefects are accumulated in the matrix. The composite is treated as an isotropic matrix strengthened by triaxial arbitrarily oriented ellipsoidal inclusions with orthotropic symmetry of the elastic properties. It is assumed that the process of loading leads to accumulation of damage in the matrix. Fractured microvolumes are modeled by a system of randomly distributed quasispherical pores. The porosity balance equation and relations for determining the effective elastic modules in the case of orthotropic components are taken as basic relations. The fracture criterion is specified as the limiting value of the intensity of average shear stresses acting in the intact part of the material. On the basis of the analytic and numerical approach, we propose an algorithm for the determination of nonlinear deformation properties of the investigated material. The nonlinearity of composite deformations is caused by the finiteness of deformations. By using the numerical solution, the nonlinear stress–strain diagrams are predicted and discussed for an orthotropic composite material for various cases of orientation of inclusions in the matrix.     

Application of the method of conditional moments to investigating the deformation properties of orthotropic composites with fiber microdamages

A model of deformation of stochastic composites subjected to microdamage is developed for the case of orthotropic materials with microdamages accumulating in the fibers. The composite is treated as a matrix strengthened with elliptic fibers with orthotropic elastic properties. The fractured microvolumes are modeled by a system of randomly distributed quasi-spherical pores. The porosity balance equation and relations for determining the effective elastic moduli for the case of a fibrous composite with orthotropic components are used as the fundamental relations. The fracture criterion is given as a limit value of the intensity of average shear stresses occurring in the undamaged part of the material, which is assumed to be a random function of coordinates and is described by the Weibull distribution. Based on an analytical and numerical approach, the algorithm for determining the nonlinear deformation properties of such a material is constructed. The nonlinearity of composite deformations is caused by the accumulation of microdamages in the fibers. By using a numerical solution, the nonlinear stress–strain diagrams for an orthotropic composite in uniaxial tension are obtained. Translated from Mekhanika Kompozitnykh Materialov, Vol. 45, No. 1, pp. 17–30, January–February, 2009.     

Simulation of micro-cutting considering finite deformation crystal plasticity

Micro-machining processes on metalic microstructures are influenced by the crystal structure, i. e. the grain orientation. Furthermore, the chip formation underlies large deformations. To perform finite element simulations of micro-cutting processes, a large deformation material model is necessary in order to model the hyperelastic and finite plastic material behaviour. In the case of cp-titanium material with hcp-crystal structure the anisotropic behaviour must be considered by an appropriate set of slip planes and slip directions. In the present work the impact of the grain orientation on the plastic deformation is demonstrated by means of finite element simulations of a finite deformation single slip crystal plasticity model. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Computational Homogenisation of Polycrystalline Elastoplastic Microstructures at Finite Deformation

During sheet bulk metal forming processes both, flat geometries and three-dimensional structures change their shape significantly while undergoing large plastic deformations. As for forming processes, FE-simulations are often done before in situ experiments, a very accurate material model is required, performing well for a huge variety of different geometrical characteristics. Because of the crystalline nature of metals, anisotropies have to be taken into account. Macroscopically observable plastic deformation is traced back to dislocations within considered slip systems in the crystals causing plastic anisotropy on the microscopic and the macroscopic level. A finite crystal plasticity model is used to model the behaviour of polycrystalline materials in representative volume elements (RVEs) of the microstructure. A multiplicative decomposition of the deformation gradient into elastic and plastic parts is performed, as well as a volumetric-deviatoric split of the elastic contribution. In order to circumvent singularities stemming from the linear dependency of the slip system vectors, a viscoplastic power-law is introduced providing the evolution of the plastic slips and slip resistances. The model is validated with experimental microstructural data under deformation. Through homogenisation and optimisation techniques, effective stress-strain curves are determined and can be compared to results from real forming processes leading to a suitable effective material model. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Deformation of orthotropic composites with unidirectional ellipsoidal inclusions under matrix microdamages

In the present paper, a model of deformation of stochastic composites under microdamaging is developed for the case of orthotropic composite, when the microdamages are accumulated in the matrix. The composite is treated as an isotropic matrix strengthened by three-axial ellipsoidal inclusions with orthotropic symmetry of elastic properties. It is assumed that the loading process leads to accumulation of damages in the matrix. Fractured microvolumes are modeled by a system of randomly distributed quasispherical pores. The porosity balance equation and relations for determining the effective elastic moduli for the case of a composite with orthotropic components are taken as the basic relations. The fracture criterion is assumed to be given as the limit value of the intensity of average shear stresses occurring in the undamaged part of the material. Based on the analytical and numerical approach, an algorithm for the determination of nonlinear deformation properties of such a material is constructed. The nonlinearity of composite deformations is caused by the accumulation of microdamages in the matrix. Using the numerical solution, nonlinear stress-strain diagrams for the orthotropic composite in the case of biaxial extension are obtained. Published in Matematychni Metody ta Fizyko-Mekhanichni Polya, Vol. 51, No. 1, pp. 121–130, January–March, 2008.     

A Formulation of Additive Finite Anisotropic Thermo-Plasticity in Logarithmic Lagrangean Strain-Entropy Space

A phenomenological model of rate-independent thermo-plasticity at finite strains is discussed. The formulation is based on an additive decomposition of the strain measure into an elastic and plastic part as proposed by Green and Naghdi. A constitutive model in the logarithmic Lagrangean strain-entropy space is developed capable of modelling isotropic elastic and anisotropic plastic material behaviour. The staggered solution scheme for coupled thermo-mechanical problems employs an isentropic phase for the deformation and an iso-geometrical phase for the thermal field. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical model of one possibility of determining the tensor of elastic constants of an anisotropic body

An algorithm for a digital computer is presented which permits calculating the tensor of the elastic constants of a material with any given symmetry on the basis of experimentally measured Young moduli and Poisson ratios for several specimens oriented in unequivalent crystallographic directions. The results of determining the components of the tensor of elastic constants by means of it for a modification of orthotropic reinforced plastic are also given.Translated from Mekhanika Polimerov, Vol. 9, No. 3, pp. 436–442, May–June, 1973.     

Calculation of the Thermoelastoplastic Nonaxisymmetric Stress-Strain State of Layered Orthotropic Shells of Revolution

The methods for determining the nonaxisymmetric thermoelastoplastic stress-strain state of layered orthotropic shells of revolution are developed. It is assumed that the layered package deforms without mutual slippage or separation of layers. The problem is solved using the geometrically nonlinear theory of shells based on the Kirchhoff-Love hypotheses. In the isotropic layers, plastic deformations may appear, whereas the orthotropic layers deform in the elastic region. It is assumed that the mechanical properties of the materials are temperature-dependent. The thermoplasticity equations are presented in a form corresponding to the method of additional deformations. The order of the system of partial differential equations obtained is reduced with the help of trigonometric series in the circumferential coordinate. The resulting systems of ordinary differential equations are solved by the Godunov technique of discrete orthogonalization. The nonaxisymmetric thermoelastoplastic stress-strain states of layered shells of revolution are considered as examples.     

Coupling of a structural analysis and flow simulation for short-fiber-reinforced polymers: property prediction and transfer of results

The aim of this paper is to discuss the basic theories of interfaces able to transfer the results of an injection molding analyis of fiber-reinforced polymers, performed by using the commercial computer code Moldflow, to the structural analysis program ABAQUS. The elastic constants of the materials, such as Young’s modulus, shear modulus, and Poisson’s ratio, which depend on both the fiber content and the degree of fiber orientation, were calculated not by the usual method of “orientation averaging,” but with the help of linear functions fitted to experimental data. The calculation and transfer of all needed data, such as material properties, geometry, directions of anisotropy, and so on, is performed by an interface developed. The interface is suit able for midplane elements in Moldflow. It calculates and transfers to ABAQUS all data necessary for the use of shell elements. In addition, a method is described how a nonlinear orthotropic behavior can be modeled starting from the generalized Hooke’s law. It is also shown how such a model can be implemented in ABAQUS by means of a material subroutine. The results obtained according to this subroutine are compared with those based on an orthotropic, linear, elastic simulation. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 45, No. 3, pp. 367-378, May-June, 2009.     

Kinematics of finite elastoplastic deformations

In this paper we review various approaches to the decomposition of total strains into elastic and nonelastic (plastic) components in the multiplicative representation of the deformation gradient tensor. We briefly describe the kinematics of finite deformations and arbitrary plastic flows. We show that products of principal values of distortion tensors for elastic and plastic deformations define principal values of the distortion tensor for total deformations. We describe two groups of methods for decomposing deformations and their rates into elastic and nonelastic components. The methods of the first group additively decompose specially built tensors defined in a common basis (initial, current, or “intermediate”). The second group implies a certain relation connecting tensors that describe elastic and plastic deformations. We adduce an example of constructing constitutive relations for elastoplastic continuums at large deformations from thermodynamic equations.     

Computational modeling of cardiac tissue with strongly coupled electromechanics and orthotropic viscoelastic effects

Modeling of complex mechanisms leading to the functioning of the heart has been an active field of research since decades. Difficulties associated with in vivo experiments motivate the utilization of computational models in order to gain a better appreciation of heart electromechanics. Although rate dependent behaviour of the orthotropic passive heart tissue has been comprehensively studied in the literature [1], effects of this phenomenon on fully coupled cardiac electromechanics are unrevealed yet. Therefore, this contribution is concerned with the investigation of viscous effects on the electromechanical response of the myocardium. To this end, we adopt the fully implicit finite element framework which strongly couples the mechanical and electrophysiological problem of the myocardium in a mono- and bi-domain setting [2,3], respectively. Viscous effects, however, are consistently embedded into this framework by making use of the orthotropic viscoelastic material model for the passive myocardium, which considers different relaxation mechanisms for the different orientation directions [5]. The performance of the proposed model is assessed by comparing finite element simulations of spiral waves in heart tissue for elastic and viscoelastic formulations. We further investigate the influence of viscosity on the defibrillation phenomenon by means of the finite element formulation of bidomain electrophysiology. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Elastoplastic torsion of a cylindrical rod for finite deformations

A solution of the problem of the torsion of a cylindrical rod was obtained in /1/ for a general, isotropic, incompressible elastic material. The present paper gives an analytical solution of the elastoplastic torsion problem for finite deformations, written in terms of quadratures of elliptic functions. The non-linear kinematics of elastoplastic deformation is introduced into the defining equations with the help of a multiplicative decomposition of the deformation gradient into elastic and plastic components /2, 3/. The elastic deformation and rate of plastic deformation are related to the state of stress of the body, in accordance with the defining Mooney-Rivlin equations /4/ and the law of flow for finite deformations associated with the Tresca yield condition /5/. A non-linear first-order partial differential equation and the initial data at the elastoplastic boundary are obtained in order to determine the angle of rotation within the plastic zone of the basis formed from the eigenvectors of the stress tensor, relative to the radial direction. The integration of the resulting equation is reduced to determining the general integral of the Ricatti equation with right-hand side determined from the angular velocity of flow of the material within the plastic zone. It is shown that neglecting the finiteness of the deformation leads to too high an estimate of the rigidity of the rod.     

Steel arches subjected to blast loading: A non-discretisation analysis approach

An innovative numerical non-discretisation semi-analytical methodology for the non-linear dynamic analysis of circular steel arch members subjected to blast loading is presented in this paper. The steel arch has a singly-symmetric cross-section with both elastic and plastic domains to account for the spread of yielding, while it is restrained at its two ends by translational springs in both the horizontal and vertical directions as well as counterpart rotational springs which simulate semi-rigid connections. The rate-dependent effects of the steel material due to rapid dynamic loading, as well as geometric non-linearities and material non-linearities are taken into account in the analysis. The effect of the included angle on the dynamic behaviour of a circular steel arch is investigated comprehensively, while the proposed methodology is validated against ABAQUS finite elements, for which the results show that the developed formulation is accurate in capturing the behaviour of the steel arch member subjected to blast loading, providing an efficacious formulation for further structural design and evaluation.     

正交各向异性材料的混合硬化弹塑性本构方程

建立了混合硬化正交各向异性材料的屈服准则,进而推导了与之相关的塑性流动法则.根据简单应力状态的实验曲线,可得到广义等效应力-应变关系.初始屈服曲面与材料的弹性常数有关,材料退化为各向同性且只考虑各向同性硬化时,屈服函数退化为Huber-Mises屈服函数,相关的本构方程退化为Prandtl-Reuss方程.     

Torsional stability of a glass-reinforced plastic cylindrical shell with an elastic core

The problem of the stability of a glass-reinforced plastic cylindrical shell with an elastic core subjected to twisting moments applied to the edges of the shell is considered. As in various other studies [4–6], the glass-reinforced plastic is treated as an elastically orthotropic material. The core is treated as an isotropic elastic cylinder, whose outer surface is bonded to the shell. Expressions for the critical stresses are obtained for an infinitely long shell and a shell of finite length.Moscow. Translated from Mekhanika Polimerov, No. 6, pp. 1082–1086, November–December, 1970.     

From Dislocation Motion to an Additive Velocity Gradient Decomposition, and Some Simple Models of Dislocation Dynamics

A mathematical theory of time-dependent dislocation mechanics of unrestricted geometric and material nonlinearity is reviewed. Within a ``small deformation" setting, a suite of simplified and interesting models consisting of a nonlocal Ginzburg Landau equation, a nonlocal level set equation, and a nonlocal generalized Burgers equation is derived. In the finite deformation setting, it is shown that an additive decomposition of the total velocity gradient into elastic and plastic parts emerges naturally from a micromechanical starting point that involves no notion of plastic deformation but only the elastic distortion, material velocity, dislocation density and the dislocation velocity. Moreover, a plastic spin tensor emerges naturally as well.