Shakedown theory for elastic plastic kinematic hardening bodies

Shakedown static and kinematic theorems for elastic–plastic (generally nonlinear) kinematic hardening solids are derived in classical (path-independence) spirit with new constructions. The generally plastic-deformation-history-dependent hardening curve is assumed to be limited by the initial yield stress and ultimate yield strength, and to obey a positive hysteresis postulate for closed plastic cycles, but else can be arbitrary and unspecified. The theorems reveal that the shakedown of structures is not affected by the particular form of the hardening curve, but just by the initial and ultimate yield stresses. While the ultimate yield strength is clearly defined macroscopically and attached to the incremental collapse mode with unbounded plastic deformations, the initial yield stress, which is responsible for the bounded cyclic plasticity collapse mode, should not be taken as the convenient one at a fixed amount of plastic deformation (0.2%), but is suggested to be taken as low as the fatigue limit to preserve the classical load-history-independence spirit of the shakedown theorems. Otherwise, for our pragmatic application purpose, it may be given empirical values between the low fatigue limit and high ultimate yield stresses according to particular loading processes considered, which may range anywhere between the high-cycle and low-cycle ones. The theorems appear as simple as those of Melan and Koiter for perfect plasticity but applied to the much larger class of more realistic kinematic hardening materials.     

Modified Burzynski criterion with non-associated flow rule for anisotropic asymmetric metals in plane stress problems

The Burzynski criterion is developed for anisotropic asymmetric metals with the non-associated flow rule (NAFR) for plane stress problems. The presented pressure depending on the yield criterion can be calibrated with ten experimental data, i.e., the tensile yield stresses at 0°, 45°, and 90°, the compressive yield stresses at 0°, 15°, 30°, 45°, 75°, and 90° from the rolling direction, and the biaxial tensile yield stress. The corresponding pressure independent plastic potential function can be calibrated with six experimental data, i.e., the tensile R-values at 0°, 15°, 45°, 75°, and 90° from the rolling direction and the tensile biaxial R-value. The downhill simplex method is used to solve these ten and six high nonlinear equations for the yield and plastic potential functions, re- spectively. The results show that the presented new criterion is appropriate for anisotropic asymmetric metals.     

NUMERICAL SIMULATION OF NON-PLANAR PLASTIC ANISOTROPY OF COLD ROLLING POLYCRYSTALLINE MATERIALS

The plastic anisotropy of sheet metal is usually caused by preferred orientation of grains, developed by mechanical deformation and thermal treatment. In the present study, a Taylor-like polycrystal model suggested by Asaro and Needleman is applied to investigate the evolution of the anisotropic behavior of a face centered cubic (FCC) polycrystalline metal, which is considered having {111} (110) slip systems, by stretching it along an arbitrary direction after it has undergonea plane-strata compression that rationally simulates the cold rolling process of FCC polycrystalline pure aluminium. By using the Taylor-like polycrystal model, pole figures are obtained to describe the texture development of polycrystalline aggregate after plane-strain compression, and then the plastic anisotropy of polycrystalline aggregate is evaluated by stretching the polycrystalline aggregate in different direction in term of yield stress. According to the results, the contours of longitudinal flow stress in three-dimensional orientation space are given and analyzed. Experiment results similar to the prediction of planar anisotropy can be found inthe literature written by Takahashi et al. that in directly show the correctness of the prediction of non-planar plastic anisotropy by this analysis.     

基于有限变形单晶塑性滑移的微观热-力耦合模型及其有限元计算

基于率相关的晶体塑性滑移理论,论文考虑晶体内部塑性变形产生的热以及快速热冲击作用下温度急剧变化产生热应力的热力双向耦合效应,建立起微观单晶的瞬态热-弹-塑性耦合模型,推导出与温度有关的剪应变率和弹塑性切模量公式.根据论文建立的模型,对ABAQUS软件进行二次开发[1],数值模拟出/{100}单晶Cu在单轴拉伸状态下的应力、应变与温度的关系和弹性模量的变化,结果如下:轴向应力随温度升高先呈线性增加再呈非线性减小,轴向应变随温度增加而增加;弹性模量随塑性变形的增加而降低,与分子动力学模拟的趋势[2]是一致的.数值实验表明,论文建立的模型和算法是正确合理的,且计算量远远小于分子动力学模拟.     

A constitutive model for orthotropic elasto-plasticity at large strains

Summary The paper presents a thermodynamically consistent constitutive model for elasto-plastic analysis of orthotropic materials at large strain. The elastic and plastic anisotropies are assumed to be persistent in the material but the anisotropy axes can undergo a rigid rotation due to large plastic deformations. The orthotropic yield function is formulated in terms of the generally nonsymmetric Mandel stress tensor such that its skew-symmetric part is additionally taken into account. Special attention is focused on the convexity of the yield surface resulting in the nine-dimensional stress space. Of particular interest are new convexity conditions which do not appear in the classical theory of anisotropic plasticity. They impose additional constraints on the material constants governing the plastic spin. The role of the plastic spin is further studied in simple shear accompanied by large elastic and large plastic deformations. If the plastic spin is neglected, the shear stress response is characterized by oscillations with an amplitude strictly dependent on the degree of the plastic anisotropy.accepted for publication 2 March 2004     

Mathematical model for cold rolling based on energy method

There is an increasing requirement for improved accuracy of the rolling models which are widely used in rolling plant in order to produce high-quality products, because this accuracy is important for rolling schedule setup and automatic control. A three-dimensional mathematical model for cold rolling using tangential velocity field and energy method is firstly proposed to investigate the deformation of the strip at the roll gap based on elastic and plastic mechanics. The field and geometrical approximation yield criterion are used to integrate the internal plastic deformation power. The friction power is analyzed using the co-line vector inner product method. Finally, the analytical expressions of roll separating force and roll torque are obtained quickly considering the effect of roll flattening on the roll separating force. The predicted roll separating forces are consistent with other researchers’ models, especially for on-line measured roll separating forces in a tandem cold rolling plant. The accurate predicted results provide valuable guidelines for optimization of rolling process.     

Laser heating of multilayer assembly and stress levels: elasto-plastic consideration

Laser heating of multilayer assembly results in temperature profiles, which differ in each layer. Depending upon the material and laser pulse properties, high temperature gradients can occur. This, in turn, results in excessive stress levels in the region irradiated by a laser beam. When the magnitude of stress level exceeds the yield stress of the material, plastic deformation is resulted. In this case, elasto-plastic analysis should be applied when modeling the thermal stresses due to a laser heating pulse. In the present study, time exponentially varying pulse laser heating of four and two layer assemblies is considered. The four layer assembly composes of gold, chromium, gold and silicon while two layer assembly is formed from gold and silicon. The resulting temperature field is obtained numerically using a control volume approach. The elasto-plastic analysis employed to compute the stress levels in the substrate material. It is found that stress levels higher than the yield stress of the substrate material occurs in the surface region. This, in turn, results in plastic zone in this region. The size of the plastic region extends towards the solid bulk as the heating progresses.     

Texture development and plastic anisotropy of B.C.C. strain hardening sheet metals

A Taylor-like polycrystal model is adopted here to investigate the plastic behavior of body centered cubic (b.c.c.) sheet metals under plane-strain compression and the subsequent in-plane biaxial stretching conditions. The <111> pencil glide system is chosen for the slip mechanism for b.c.c. sheet metals. The {110} <111> and {112} <111> slip systems are also considered. Plane-strain compression is used to simulate the cold rolling processes of a low-carbon steel sheet. Based on the polycrystal model, pole figures for the sheet metal after plane-strain compression are obtained and compared with the corresponding experimental results. Also, the simulated plane-strain stress—strain relations are compared with the corresponding experimental results. For the sheet metal subjected to the subsequent in-plane biaxial stretching and shear, plastic potential surfaces are determined at a given small amount of plastic work. With the assumption of the equivalence of the plastic potential and the yield function with normality flow, the yield surfaces based on the simulations for the sheet metal are compared with those based on several phenomenological planar anisotropic yield criteria. The effects of the slip system and the magnitude of plastic work on the shape and size of the yield surfaces are shown. The plastic anisotropy of the sheet metal is investigated in terms of the uniaxial yield stresses in different planar orientations and the corresponding values of the anisotropy parameter R, defined as the ratio of the width plastic strain rate to the through-thickness plastic strain rate under in-plane uniaxial tensile loading. The uniaxial yield stresses and the values of R at different planar orientations from the polycrystal model can be fitted well by a yield function recently proposed by Barlat et al. (1997b).     

塑性材料的统一屈服条件

从一个新的角度对经典屈服条件进行研究，并以最大剪应力是使塑性屈服发生的主要是否 是唯一因素为出发点，提出了包含材料特性以参数的统一屈服条件。     

硬化系数对界面端弹塑性奇异应力场的影响

本文利用弹塑性边界元分析方法,对具有不同硬化系数的线性硬化结合材料界面端进行了计算,分析结果表明,当硬化系数较大时,界面附近的弹塑怀应力与将弹塑性本构关系简化为线性后得到的理论结果相接近,而当硬化系数相对较少时,理论分析的奇异应力场的主控区变得非常小,在屈服域的绝大部分区间,应力奇异性与理论解有较大区别,本文的结果还表明,硬化系数越小,过渡区（弹塑性厅异应力场支配区到屈服边界）越大,屈服区域应力分布变得平坦,在小规模屈服条件异次数一致）,即可用弹性厅异应力场来近似地描述小规模屈服时的弹塑性界面端,但应力强度系数则比弹性时略大,且随硬化系数的减小而增大。     

Creep transition in torsion

Stresses for a circular cylinder of compressible material subjected to torsion are derived in closed form for steady state creep. It is shown that the asymptotic solution through stress leads from elastic state to plastic and then to creep and through stress difference leads to the creep state. The effect of compressibility is presented graphically. The results indicate that the value of maximum shear stress for a cylinder of compressible material is greater than that for an incompressible material and increases with the increase in a measure index n. For an incompressible material, as a particular case, the results obtained are the same as given by Marin [9].     

Yield behavior of photoplastic materials

A mixture of flexible and rigid polyester resins has been used previously by Morris and Riley¹⁶ and by Zachary and Riley⁹ to model plastic deformations. The last of these papers furnished mechanical and optical properties under uniaxial tension and compression for several different mixture ratios of the polyester resins and also presents some results under multiaxial-stress conditions from thin cylinders under internal pressure. In a recent paper, Burger, Gomide and Scott¹⁴ used the rigid polyester resin at elevated temperature to model plastic deformations in upset rings; the behavior of the rigid polyester was verified with diametrically compressed disks and uniaxial-compression specimens. A very important similitude requirement for model to prototype scaling in photoplasticity work is that the macroscopic yield behavior of model and prototype materials must be the same. Thus, not only uniaxial tension and compression properties must be examined, but also yield properties under multiaxial-stress states have to be determined. The purpose of this paper is to provide additional information on the yield behavior of polyester mixtures which appear suitable for model studies of manufacturing methods such as rolling and extruding. For these processes, mixture ratio, test temperature and strain rate can be used to control the shape of the stress-strain curve and the yield behavior. The experimental procedure used to determine the initial yield locus of the photoplastic materials employed a new specimen geometry proposed by Arcan, Hashin and Voloshin¹⁸ which produces uniform biaxial-stress fields of opposite sign in one section of the specimen. Both polycarbonate and polyester materials were evaluated using this procedure and results are compared with those available in the technical literature.     

A generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials

In this paper, a generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials is derived. The evolution equation for the active yield surface with reference to the memory yield surface is obtained by considering the continuous expansion of the active yield surface during the unloading/reloading process. The incremental constitutive relation based on the associated flow rule is then derived for a general yield function for pressure insensitive and sensitive materials. Detailed incremental constitutive relations for materials based on the Mises yield function, the Hill quadratic anisotropic yield function and the Drucker–Prager yield function are derived as the special cases. The closed-form solutions for one-dimensional stress–plastic strain curves are also derived and plotted for materials under cyclic loading conditions based on the three yield functions. In addition, the closed-form solutions for one-dimensional stress–plastic strain curves for materials based on the isotropic Cazacu–Barlat yield function under cyclic loading conditions are summarized and presented. For materials based on the Mises and the Hill anisotropic yield functions, the stress–plastic strain curves show closed hysteresis loops under uniaxial cyclic loading conditions and the Masing hypothesis is applicable. For materials based on the Drucker–Prager and Cazacu–Barlat yield functions, the stress–plastic strain curves do not close and show the ratcheting effect under uniaxial cyclic loading conditions. The ratcheting effect is due to different strain ranges for a given stress range for the unloading and reloading processes. With these closed-form solutions, the important effects of the yield surface geometry on the cyclic plastic behavior due to the pressure-sensitive yielding or the unsymmetric behavior in tension and compression can be shown unambiguously. The closed form solutions for the Drucker–Prager and Cazacu–Barlat yield functions with the associated flow rule also suggest that a more general anisotropic hardening theory needs to be developed to address the ratcheting effects for a given stress range.     

Correlation between swift effects and tension–compression asymmetry in various polycrystalline materials

The Swift phenomenon, which refers to the occurrence of permanent axial deformation during monotonic free-end torsion, has been known for a very long time. While plastic anisotropy is considered to be its main cause, there is no explanation as to why in certain materials irreversible elongation occurs while in others permanent shortening is observed.In this paper, a correlation between Swift effects and the stress–strain behavior in uniaxial tension and compression is established. It is based on an elastic–plastic model that accounts for the combined influence of anisotropy and tension–compression asymmetry. It is shown that, if for a given orientation the uniaxial yield stress in tension is larger than that in compression, the specimen will shorten when twisted about that direction; however, if the yield stress in uniaxial compression is larger than that in uniaxial tension, axial elongation will occur. Furthermore, it is shown that on the basis of a few simple mechanical tests it is possible to predict the particularities of the plastic response in torsion for both isotropic and initially anisotropic materials. Unlike other previous interpretations of the Swift effects, which were mainly based on crystal plasticity and/or texture evolution, it is explained the occurrence of Swift effects at small to moderate plastic strains. In particular, the very good quantitative agreement between model and data for a strongly anisotropic AZ31–Mg alloy confirm the correlation established in this work between tension–compression asymmetry and Swift effects. Furthermore, it is explained why the sign of the axial plastic strains that develop depends on the twisting direction.     

An experimental study of uniaxial creep,cyclic creep and relaxation of aisi type 304 stainless steel at room temperature

Following previous work (Krempl, 1979), a servocontrolled testing machine and strain measurement at the gage length were used to study the uniaxial rate(time)-dependent behavior of AISI Type 304 stainless steel at room temperature. The test results show that the creep strain accumulated in a given period of time depends strongly on the stress-rate preceding the creep test. In constant stress-rate zero-to-tension loading the creep strain accumulated in a fixed time-period at a given stress level is always higher during loading than during unloading. Continued cycling causes an exhaustion of creep ratchetting which depends on the stress-rate. Periods of creep and relaxation introduced during completely reversed plastic cycling show that the curved portions of the hysteresis loop exhibit most of the inelasticity. In the straight portions, creep and relaxation are small and there exists a region commencing after unloading where the behavior is similar to that at the origin for virgin materials. This region does not extend to zero stress.The results are at variance with creep theory and with viscoplasticity theories which assume that the yield surface expands with the stress. They support the theory of viscoplasticity based on total strain and overstress.     

Basic equations for the microscopic theory of plasticity of polycrystalline materials with given texture

An expression for the yield stress of anisotropic materials is applied to the anisotropic strength of hard rolled copper foils whose crystallographic texture is known. We assume that this crystallographic texture is the only cause of the anisotropic plastic behaviour of the material. The model used for the yield stress is also used to deduce:

1. Stress-strain relations for isotropic polycrystalline materials;
2. A formula for the fully plastic strain tensor, applied to anisotropic hard rolled copper foils.

For the anisotropic copper foils considered the calculated curves of the yield stress and of the strain tensor as a function of the angle x between rolling and tensile direction agree qualitatively with the measured values. However, the theory is not complete, since the yield stress and the plastic strain tensor are both a function of a parameter Q, the fraction of the number of available crystallographic slip planes on which the maximum shear stress has reached the critical value τ_a. We assume that for “fully” plastic deformation a certain critical fraction Q _e of the total number of slip planes has to be active. The fraction Q _e is called the critical active quantity. With the parameter Q _e we adjust the calculated curves to the measured ones. The dependence of Q _e on the properties of the material (e.g. the crystallographic texture) is discussed in Appendix I.     

Variational methods for the steady state response of elastic–plastic solids subjected to cyclic loads

Solids (or structures) of elastic–plastic internal variable material models and subjected to cyclic loads are considered. A minimum net resistant power theorem, direct consequence of the classical maximum intrinsic dissipation theorem of plasticity theory, is envisioned which describes the material behavior by determining the plastic flow mechanism (if any) corresponding to a given stress/hardening state. A maximum principle is provided which characterizes the optimal initial stress/hardening state of a cyclically loaded structure as the one such that the plastic strain and kinematic internal variable increments produced over a cycle are kinematically admissible. A steady cycle minimum principle, integrated form of the aforementioned minimum net resistant power theorem, is provided, which characterizes the structure’s steady state response (steady cycle) and proves to be an extension to the present context of known principles of perfect plasticity. The optimality equations of this minimum principle are studied and two particular cases are considered: (i) loads not exceeding the shakedown limit (so recovering known results of shakedown theory) and (ii) specimen under uniform cyclic stress (or strain). Criteria to assess the structure’s ratchet limit loads are given. These, together with some insensitivity features of the structure’s alternating plasticity state, provide the basis to the ratchet limit load analysis problem, for which solution procedures are discussed.     

Polar plasticity of thin-walled composites made of ideal fibre-reinforced material

The present study investigates the influence that polar material response has on the plastic behaviour of thin-walled structures made of ideal fibre-reinforced materials (Spencer, 1972); or, equivalently, on the response of thin-walled fibrous composites within the first branch of the matrix dominated form (MDM) of the bimodal theory of plasticity (Soldatos, 2011, Dvorak and Bahei-El-Din, 1987). The plasticity studies mentioned above assume that fibres are infinitely thin and, therefore, perfectly flexible. They possess no bending stiffness and, hence, their negligible bending resistance cannot influence the developed stress state, which is accordingly described by a symmetric stress tensor. In contrast, the present study considers that if fibres resistant in bending are embedded in a material at high volume concentrations, their flexure produces couple-stress and, as a result of this kind of polar material behaviour, the stress tensor becomes non-symmetric. Under plane stress conditions that dominate behaviour of thin-walled structures, the stress-space and, therefore, conditions of plastic yield and relevant yield surfaces are thus four-dimensional. However, shapes and properties of initial yield surfaces relevant to the f₁-branch of MDM are studied comprehensively by considering their projection on particular planes of such a four-dimensional stress-space. It then becomes easier understood that, in the regime of polar material response, a thin-walled structure made of ideal fibre-reinforced material deforms plastically when suitable combinations of shear stress values are reached simultaneously, rather than when only one of two unequal shear stress components reaches some maximum absolute value. Thus, polar material plasticity dismisses the conventional concept of material yield stress in shear and replaces it with a pair of two independent yield moduli. Existence of the latter is perceived as a theoretical justification of the expectation that, due to the presence of fibres, two rather than one shear yield parameters of the composite should be present and accountable for. The non-zero values of those parameters are shown to exert paramount influence on the form of the yield surface of the ideal fibre-reinforced material of interest.     

Technological thermal stresses in the shrink fitting of cylindrical bodies with consideration of plastic flows

This paper presents a solution of a sequence of one-dimensional boundary-value problems of thermal stresses defining the elastic–plastic deformation processes used in the shrink fitting of cylindrical bodies. The initiation and development of plastic flow in the materials of the assembly elements are studied taking into account the temperature dependence of the yield stress of these materials. During temperature equalization, the flow can slow down, followed by unloading and formation of a residual stress field providing tension. The conditions of formation and motion of the boundaries of the elastic and plastic states in plastic flow and during unloading are determined.