Nonisothermal instability of high-velocity elastoplastic flows

An important feature of the high-velocity deformation of solids is the localization of deformation, one of the causes of which may be the nonisothermal instability of plastic flow [1–6]. In connection with the intensive development of high-velocity technology in the treatment of materials, the investigation of the criteria for nonisothermal stability of processes of plastic deformation is of fundamental interest, since in certain cases they determine the optimum technological regimes [5]. The critical values of deformation velocities, above which the effects of thermal instability becomes decisive in the process of deformation of solids, are estimated by semiempirical methods in [1]. The non-boundary-value problem of the criteria for nonisothermal instability is analyzed in [2] for the point of view of flow stability in the so-called coupled formulation. The latter means that the heat-conduction equation is added to the basic equations determining the dynamics of an elastoplastic medium. The problem is solved in [6] in an analogous formulation, but for flow averaged over the spatial coordinate. The solution of the boundary-value problem for one-dimensional flow in this formulation is given in the present paper.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 133–138, May–June, 1986.     

Some discussions on finite deformation of continuous media

Three problems are discussed in this paper:1.The physical significance of the characteristic tensorθ~ii of finitedeformation is discussed as a complement to the literature.2.The four characterietic tensors of finite deformation,introduced inthe.literature,ard analyzed and discussed more thoroughly.3.The representation of the general finite deformation through simple.loading process is not always possible,the condition for its realization beingthat the given finite deformation satisfies the equations of compatibility.It is pointed out in this paper that for the finite deformation set forth byL.I.Sedov in an illustrative example in[9],the equations of compatibilitycannot be satisfied even for k=1,henceforth for whatever value of k,thefinite deformation set forth by Sedov can not be represented through simpleloading process.     

ON UNILATERAL CONTACT LARGE DEFORMATION PROBLEM WITH FRICTION(Ⅱ)-NONLINEAR FINITE ELEMENT TECHNIQUE AND ITS APPLICATION

Based on the variational equation derived in ref. [1], a nonlinear incremental F.E. equation is formulated for unilateral contact elastic and plastic large deformation problems. A new technique-co-moving coordinate finite element method is introduced, and a practical mathematical model for large deformation contact problem is described. To show the effectiveness of the method, problems of contact large deformation of cantilever beam, circular plate, as well as metal ring are computed. Compared with experiments, the results show good agreements.     

In-plane buckling of prismatic funicular arches with shear deformations

The in-plane buckling behavior of funicular arches is investigated numerically in this paper. A finite strain Timoshenko beam-type formulation that incorporates shear deformations is developed for generic funicular arches. The elastic constitutive relationships for the internal beam actions are based on a hyperelastic constitutive model, and the funicular arch equilibrium equations are derived. The problems of a submerged arch under hydrostatic pressure, a parabolic arch under gravity load and a catenary arch loaded by overburden are investigated. Buckling solutions are derived for the parabolic and catenary arch. Subsequent investigation addresses the effects of axial deformation prior to buckling and shear deformation during buckling. An approximate buckling solution is then obtained based on the maximum axial force in the arch. The obtained buckling solutions are compared with the numerical solutions of Dinnik (Stability of arches, 1946) [1] and the finite element package ANSYS. The effects of shear deformation are also evaluated.     

Deformation of liquid line with the collision of jets

With the study of the collision of metallic plates accelerated by an explosion, investigation of plastic deformations in the collision zone is of great importance. In [1] a method is expounded for investigations of deformation with explosion welding under conditions of wave formation; the method consists of the pressing of a thin wire into the metallic plate. With collision of the plates, the flowing metal carries the wire with it; the change in the form of the latter makes it possible to evaluate the character of the plastic deformation in the collision zone. From an investigation of the deformation of the wires, such important characteristics as the viscosity of the metals can be found. The method set forth in [1] is one of the few means of investigating the deformation of metals with high-speed collisions. The difficulties in the investigation are bound up with the small times of the process and with the high pressures, developing with an explosion and demolishing the experimental unit. In [1], from an analysis of viscous flow during the collision of plates, a dependence of the shift on the distance to the interface between the materials, described by a parabola, is derived. It is noted that, near the interface, the theoretical and experimental results differ considerably. It is important to make an attempt at a theoretical analysis of the wire deformation if, with the collision of metallic plates, a jet is formed and the interface between the materials is even. In this case, the flow differs strongly from the flow under conditions of wave formation [2], and near the interface the flow cannot be described by a parabolic dependence; this dependence will probably be exponential. To describe flow with the formation of a jet, in [3] a model of an ideal liquid is employed. Since the analysis of a collision between viscous jets with a free boundary is bound up with serious difficulties, in the first approximation it is advisable to consider the deformation of the liquid line with the deformation of ideal jets.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 65–71, November–December, 1977.The author wishes to express his thanks to N. S. Kozin and V. V. Efremov for a number of valuable observations with an evaluation of the work, and N. M. Maksimova for her aid in making the numerical calculation.     

Deformation of lead disks under impact

The effect of relatively low-velocity (1–3 m/sec) impact on a thin disk of imcompressible viscoplastic material placed in the gap between parallel rough surfaces is considered. The state of stress of the interlayer is assumed nearly hydrostatic during impact, the duration of which is limited by the elastic deformation of the elements of the striker system. The mathematical problem of determining the distributions of stresses, velocities, and temperatures for the axisymmetric deformation of a disk is reduced to the integration of an ordinary second-order differential equation. Numerical calculations for certain cases of impact are compared with the results of experiments on lead samples. Plane strain of an interlayer of viscoplastic material between rigid plates moving with a constant velocity is discussed in [1]. The state of stress of the interlayer for the same conditions of motion of the plates was studied in [2] for axial symmetry. In the present paper we take account of the impact nature of the loading and the elastic compression of the elements of the striker system, factors on which the deformation and the pressure developed in the impact depend.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 153–158, September–October, 1975.     

有限弹塑性变形的三维组集式本构模型

本文将文[1]中提出的三维组集式弹塑性本构模型推广应用于有限变形分析,导出了全量型和增量型本构关系在初始构形上的拉格朗日(Total Lagrange)形式和瞬时构形上的拉格朗日(Updated Lagrange)形式。文中对晶体单轴拉伸中的宏观剪切带进行了分析。预测结果与实验吻合。从而说明这种本构模型能够模拟有限变形中的几何非线性效应和晶体材料塑性变形中的宏观力学行为。     

Deformation and fracture of crystalline materials under complex loading-program conditions

In predicting the deformation and conditions of fracture of materials under complex loading program conditions one must consider the entire duration of the deformation process.The problem becomes complicated in cases of a simultaneous operation of various mechanisms of deformation and fracture, e. g., when plastic deformation is superposed (once or repeatedly) on creep.A promising phenomenological approach to this problem may be based on concepts of the mechanical equation of state of materials. A hypothesis of the existence of the equation of state depending on a finite number of structural parameters was formulated by Kröner [1] for the case of the three-dimensional law of plasticity and by Rabotnov [2] for the case of creep and fracture under uniaxial stress state conditions.This article is concerned with the application of the hypothesis of the mechanical equation of state to the problem of deformation and fracture of materials (in the uniaxial case) under complex loading program conditions.     

A new uniqueness theorem of the linear thermo-elastic dynamics on unbounded domains

In the present paper, the uniqueness of the solution to the initial boundary value problem of the linear thermo-elastic dynamics on unbounded domains is obtained under less restrictive conditions, including abandoning the positive semi-definiteness of the elasticity tensor and boundness of the material tensor and restrictions on the acoustic tensor and the coupled tensor, and the results in [1] are refined. The conclusion here is valid for the case on bounded domains and the linear elastic dynamics on unbounded domains, hence the results in [2–4] are refined too. Abandoning the positive semi-definiteness of elasticity tensor permits that the uniqueness of the kinetic process is still valid for deformation of the wider materials, especially for the case that there are phase-transition during deformation process provided that the constitutive equations are unchanged in forms. The project is partially supported by The Youth Foundtion of Science of the Higher-Education of Shanghai and YFNSC (No. 19802012)     

On internal dissipation inequalities and finite strain inelastic constitutive laws: Theoretical and numerical comparisons

Internal dissipation always occurs in irreversible inelastic deformation processes of materials. The internal dissipation inequalities (specific mathematical forms of the second law of thermodynamics) determine the evolution direction of inelastic processes. Based on different internal dissipation inequalities several finite strain inelastic constitutive laws have been formulated for instance by Simo [Simo, J.C., 1992. Algorithms for static and dynamic multiplicative plasticity that preserve the classical return mapping schemes of the infinitesimal theory. Computer Methods in Applied Mechanics and Engineering 99, 61–112]; Simo and Miehe [Simo, J.C., Miehe, C., 1992. Associative coupled thermoplasticity at finite strains: formulation, numerical analysis and implementation. Computer Methods in Applied Mechanics and Engineering 98, 41–104]; Lion [Lion, A., 1997. A physically based method to represent the thermo-mechanical behavior of elastomers. Acta Mechanica 123, 1–25]; Reese and Govindjee [Reese, S., Govindjee, S., 1998. A theory of finite viscoelasticity and numerical aspects. International Journal of Solids and Structures 35, 3455–3482]; Lin and Schomburg [Lin, R.C., Schomburg, U., 2003. A finite elastic–viscoelastic–elastoplastic material law with damage: theoretical and numerical aspects. Computer Methods in Applied Mechanics and Engineering 192, 1591–1627]; Lin and Brocks [Lin, R.C., Brocks, W., 2004. On a finite strain viscoplastic theory based on a new internal dissipation inequality. International Journal of Plasticity 20, 1281–1311]; and Lin and Brocks [Lin, R.C., Brocks, W., 2005. An extended Chaboche’s viscoplastic law at finite strains: theoretical and numerical aspects. Journal of Materials Science and Technology 21, 145–147]. These constitutive laws are consistent with the second law of thermodynamics. As the internal dissipation inequalities are described in different configurations or coordinate systems, the related constitutive laws are also formulated in the corresponding configurations or coordinate systems. Mathematically, these constitutive laws have very different formulations. Now, a question is whether the constitutive laws provide identical constitutive responses for the same inelastic constitutive problems. In the present work, four types of finite strain viscoelastic and endochronically plastic laws as well as three types of J₂-plasticity laws are formulated based on four types of dissipation inequalities. Then, they are numerically compared for several problems of homogeneous or complex finite deformations. It is demonstrated that for the same inelastic constitutive problem the stress responses are identical for deformation processes without rotations. In the simple shear deformation process with large rotation, the presented viscoelastic and endochronically plastic laws also show almost identical stress responses up to a shear strain of about 100%. The three laws of J₂-plasticity also produce the same shear stresses up to a shear strain of 100%, while different normal stresses are generated even at infinitesimal shear strains. The three J₂-plasticity laws are also compared at three complex finite deformation processes: billet upsetting, cylinder necking and channel forming. For the first two deformation processes similar constitutive responses are obtained, whereas for the third deformation process (with large global rotations) significant differences of constitutive responses can be observed.     

Well production indicators in a deformable reservoir interacting with rocks

Steady-state flow towards a well through a thin porous deformable two-layer reservoir with allowance for deformation of the surrounding rocks is investigated. The permeability of the reservoir is considered to be a function of the displacements of its top and bottom. The effect of deformation on the well production indicators is studied. The results obtained agree qualitatively with the data of full-scale experiments. Earlier, in [1–5], in considering the self-consistent processes of flow through porous media and their deformation attention was concentrated on the analysis of the stress-strain state of the rocks and reservoir and on unsteady problems within the framework of the nonlocally elastic flow regime.Kazan'. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 86–93, January–February, 1995.     

On Non-Symmetric Deformations of Neo-Hookean Solids

Deformations possible (i.e., those satisfying the governing three-dimensional equations of equilibrium and the incompressibility constraint) within a class of non-symmetric deformations for a neo-Hookean nonlinearly elastic body were determined in [1], where it was found that only three special cases of the class of deformation fields considered could be solutions. One of these is the trivial solution, one the solution describing radially symmetric deformation, and the other a (non-symmetric, non-homogeneous) deformation contained within a family of universal deformations. In this paper, the results reported in [1] are shown to hold for a substantially broadened deformation field. This revised version was published online in July 2006 with corrections to the Cover Date.     

Obtaining the load-extension curves of samples subjected to an explosion load

The energy of explosives has been increasingly used in recent times for the mechanical testing of construction materials subjected to high rates of deformation [1–8]. We describe in the present paper a small apparatus for the uniaxial expansion testing of samples with the aid of explosives at deformation rates of about 10³ sec^–1; we report on some results of testing steel and aluminum alloy samples.Translated from Zhurnal Prikladnoi Mekhanika i Tekhnicheskoi Fizika, No. 1, pp. 184–186, January–February. 1975.     

Detwinning of High-Purity Zirconium: <Emphasis Type="Italic">In-Situ</Emphasis> Neutron Diffraction Experiments

Twinning is an important deformation mode in hexagonal metals to accommodate deformation along the c-axis. It differs from slip in that it accommodates shear by means of crystallographic reorientation of domains within the grain. Such reorientation has been shown to be reversible (detwinning) in magnesium alloy aggregates. In this paper we perform in-situ neutron diffraction reversal experiments on high-purity Zr at room temperature and liquid nitrogen temperature, and follow the evolution of twin fraction. The experiments were motivated by previous studies done on clock-rolled Zr, subjected to deformation history changes (direction and temperature), in the quasi-static regime, for temperatures ranging from 76 K to 450 K. We demonstrate here for the first time that detwinning of { 10[`1] 2 } á 10[`1][`1] ñ\left\{ {10\overline 1 2} \right\}\left\langle {10\overline 1 \overline 1 } \right\rangle tensile twins is favored over the activation of a different twin variant in grains of high-purity polycrystalline Zr. A visco-plastic self-consistent (VPSC) model developed previously, which includes combined slip and twin deformation, was used here to simulate the reversal behavior of the material and to interpret the experimental results in terms of slip and twinning activities.     

Initial stresses in elastic solids: Constitutive laws and acoustoelasticity

On the basis of the nonlinear theory of elasticity, the general constitutive equation for an isotropic hyperelastic solid in the presence of initial stress is derived. This derivation involves invariants that couple the deformation with the initial stress and in general, for a compressible material, it requires 10 invariants, reducing to 9 for an incompressible material. Expressions for the Cauchy and nominal stress tensors in a finitely deformed configuration are given along with the elasticity tensor and its specialization to the initially stressed undeformed configuration. The equations governing infinitesimal motions superimposed on a finite deformation are then used to study the combined effects of initial stress and finite deformation on the propagation of homogeneous plane waves in a homogeneously deformed and initially stressed solid of infinite extent. This general framework allows for various different specializations, which make contact with earlier works. In particular, connections with results derived within Biot's classical theory are highlighted. The general results are also specialized to the case of a small initial stress and a small pre-deformation, i.e. to the evaluation of the acoustoelastic effect. Here the formulas derived for the wave speeds cover the case of a second-order elastic solid without initial stress and subject to a uniaxial tension [Hughes and Kelly, Phys. Rev. 92 (1953) 1145] and are consistent with results for an undeformed solid subject to a residual stress [Man and Lu, J. Elasticity 17 (1987) 159]. These formulas provide a basis for acoustic evaluation of the second- and third-order elasticity constants and of the residual stresses. The results are further illustrated in respect of a prototype model of nonlinear elasticity with initial stress, allowing for both finite deformation and nonlinear dependence on the initial stress.     

Local and overall extremal properties of time-independent materials and non-linear elastic comparison materials

For time-independent materials which undergo non-linear deformations from some given reference configuration two (dual) hypotheses are considered. Firstly it is supposed that the work done to a given state of deformation is bounded below and that the bound is attainable on a physically possible path; secondly that the complementary work to a given state of stress is bounded above and that this bound too is attainable on a physically possible path. The consequences of these assumptions are analysed, and the results of Ponter and Martin [1] in the linear theory are generalized to account for non-linear deformations, due attention being paid to questions of stability.A non-linear elastic comparison material is defined whose strain energy is equal to the work done on a minimum path for the time-independent material. Extremum principles for non-linear elastic materials given in [2] are then applied to the comparison elastic material, and bounds are thereby placed on the work and complementary-work functional of the time-independent material. Corresponding overall properties of the time-independent and elastic materials are then compared by defining respective overall constitutive laws and overall stress and deformation variables.Following the definition of strengthening (weakening) of a non-linear elastic solid given by Ogden[2] a time-independent material is said to be strengthened (weakened) when its comparison elastic material is strengthened (weakened). Local and overall aspects of this definition are examined.     

Crystallographic orientation and size dependence of tension-compression asymmetry in molybdenum nano-pillars

Uniaxial tension and compression experiments on [0 0 1] and [0 1 1] oriented molybdenum nano-pillars exhibit tension-compression asymmetry, a difference in attained stresses in compression vs. tension, which is found to depend on crystallographic orientation and sample size. We find that (1) flow stresses become higher at smaller diameters in both orientations and both loading directions, (2) compressive flow stresses are higher than tensile ones in [0 0 1] orientation, and visa versa in [0 1 1] orientation, and (3) this tension-compression asymmetry is in itself size dependent. We attribute these phenomena to the dependence of twinning vs. antitwinning deformation on loading direction, to the non-planarity of screw dislocation cores in Mo crystals, and to the possibly lesser role of screw dislocations in governing nano-scale plasticity compared with bulk Mo.     

有限变形下多晶晶体塑性模型算法及应用

用Sanna和Zacharia^[1]所提出的延性单晶本构模型的积分算法和Taylor多晶模型假设研究了时间步长和硬化模型的选取对多晶集合体的应力应变响应和织构演化的影响。该算法是利用变形梯度乘法分解获得弹性变形梯度演化方程，用增量迭代法积分该方程，显式更新各滑移上的临界分切剪应力。算例的结果表明该算法具有时间步大，计算效率高的特点，另外，不同硬化模型的选取对多晶集合体应力应变响应的预测有明显的影响但对织构演化的预测影响不大。