Multiaxial and non-proportional loading responses,anisotropy and modeling of Ti–6Al–4V titanium alloy over wide ranges of strain rates and temperatures

Results from a series of multiaxial loading experiments on the Ti–6Al–4V titanium alloy are presented. Different loading conditions are applied in order to get the comprehensive response of the alloy. The strain rates are varied from the quasi-static to dynamic regimes and the corresponding material responses are obtained. The specimen is deformed to large strains in order to study the material behavior under finite deformation at various strain rates. Torsional Kolsky bar is used to achieve shear strain rates up to 1000 s⁻¹. Experiments are performed under non-proportional loading conditions as well as dynamic torsion followed by dynamic compression at various temperatures. The non-proportional loading experiments comprise of an initial uniaxial loading to a certain level of strain followed by biaxial loading, using a channel-type die at various rates of loadings. All the non-proportional experiments are carried out at room temperature. Experiments are also performed to investigate the anisotropic behavior of the alloy. An orthotropic yield criterion [proposed by Cazacu, O., Plunkett, B., Barlat, F., 2005. Orthotropic yield criterion for hexagonal closed packed metals. International Journal of Plasticity 22, 1171–1194.] for anisotropic hexagonal closed packed materials with strength differential is used to generate the yield surface. Based on the definition of the effective stress of this yield criterion, the observed material response for the different loading conditions under large deformation is modeled using the Khan–Huang–Liang (KHL) equation assuming isotropic hardening. The model constants used in the present study, were pre-determined from the extensive uniaxial experiments presented in the earlier paper [Khan, A.S., Suh, Y.S., Kazmi R., 2004. Quasi-static and dynamic loading responses and constitutive modeling of titanium alloys. International Journal of Plasticity 20, 2233–2248]. The model predictions are found to be extremely close to the observed material response.     

Constitutive modeling of cyclic plasticity deformation of a pure polycrystalline copper

Cyclic plasticity experiments were conducted on a pure polycrystalline copper and the material was found to display significant cyclic hardening and nonproportional hardening. An effort was made to describe the cyclic plasticity behavior of the material. The model is based on the framework using a yield surface together with the Armstrong–Frederick type kinematic hardening rule. No isotropic hardening is considered and the yield stress is assumed to be a constant. The backstress is decomposed into additive parts with each part following the Armstrong–Frederick type hardening rule. A memory surface in the plastic strain space is used to account for the strain range effect. The Tanaka fourth order tensor is used to characterize nonproportional loading. A set of material parameters in the hardening rules are related to the strain memory surface size and they are used to capture the strain range effect and the dependence of cyclic hardening and nonproportional hardening on the loading magnitude. The constitutive model can describe well the transient behavior during cyclic hardening and nonproportional hardening of the polycrystalline copper. Modeling of long-term ratcheting deformation is a difficult task and further investigations are required.     

变温下橡胶材料力学性能的实验分析

利用Zwick020材料试验机研究了天然橡胶材料在不同温度条件下单向拉伸的大变形力学行为。得到了不同温度和不同加载速率条件下材料的应力-应变关系曲线和材料的破坏条件,由此分析了高温和低温的温度变化条件及加载速率对材料力学性能的影响,得出天然橡胶材料在大变形条件下的温度和加载速率敏感特性。同时利用实验结果拟合出材料的应变能函数,并应用热超弹性模型对材料进行了理论分析。结果表明,理论分析结果与实验结果仅在一定温度范围内吻合得较好,因为材料软化或硬化,在很高或很低的温度情况下二者有一定差别。     

A viscoplastic constitutive model for nickel-base superalloy,part 2: modeling under anisothermal conditions

In Part 2 of this study, extensive deformation tests were carried out on the nickel-base polycrystalline superalloy IN738LC under isothermal and anisothermal conditions between 450 and 950 °C. Under the isothermal conditions, the material showed almost no rate/time-dependency below 700 °C, while it showed distinct rate/time-dependency above 800 °C. Regarding the cyclic deformation, slight cyclic hardening behavior was observed when the temperature was below 700 °C and the imposed strain rate was fast, whereas in the case of the temperature above 800 °C or under slower strain rate conditions, the cyclic hardening behavior was scarcely observed. Unique inelastic behavior was observed under in-phase and out-of-phase anisothermal conditions: with an increase in the number of cycles, the stress at higher temperatures became smaller and the stress at lower temperatures became larger in absolute value although the stress range was approximately constant during the cyclic loading. In other words, the mean stress continues to evolve cycle-by-cycle in the direction of the stress at lower temperatures. Based on the experimental results, it was assumed that evolution of the variable Y that had been incorporated into a kinematic hardening rule in Part 1 of this study is active under higher temperatures and is negligible under lower temperatures. The material constants used in the constitutive equations were determined with the isothermal data, and were expressed as functions of temperature empirically. The extended viscoplastic constitutive equations were applied to the anisothermal cyclic loading as well as the monotonic tension, stress relaxation, creep and cyclic loading under the isothermal conditions. It was demonstrated that the present viscoplastic constitutive model was successful in describing the inelastic behavior of the material adequately, including the anomalous inelastic behavior observed under the anisothermal conditions, owing to the consideration of the variable Y.     

Constitutive modeling of biodegradable polymers: Hydrolytic degradation and time-dependent behavior

A large range of biodegradable polymers has been used to produce implantable medical devices. Apart from biological compatibility, these devices shall be also functional compatible and perform adequate mechanical temporary support during the healing process. However, the mechanical behavior of biodegradable materials during its degradation, which is an important aspect of the design of these biodegradable devices, is still an unexplored subject. Based on the literature, the mechanical behavior of biodegradable polymers is strain rate dependent and exhibits hysteresis upon cyclic loading. On the other hand, ductility, toughness and strength of the material decay during hydrolytic degradation. In this work, it is considered a three-dimensional time-dependent model adapted from the one developed by Bergström and Boyce to simulate the performance of biodegradable structures undergoing large deformations incorporating the hydrolysis degradation. Since this model assumes that the mechanical behavior is divided into a time independent network and a non-linear time-dependent network, it enables to simulate the monotonic tests of a biodegradable structure loaded under different strain rates. The hysteresis effects during unloading–reloading cycles at different strain levels can be predicted by the model. A parametric study of the material model parameters evolution during the hydrolytic degradation was conducted to identify which parameters are more sensible to this degradation process. The investigated model could predict very well the experimental results of a blend of polylactic acid and polycaprolactone (PLA–PCL) in the full range of strains until rupture during hydrolytic degradation. From these results and analyses, a method is proposed to simulate the three-dimensional mechanical behavior during hydrolytic degradation.     

45钢的J-C损伤失效参量研究

为了在结构碰撞效应的有限元分析中描述材料行为,通过开展45钢在不同应力状态和温度下的准静态材料力学性能实验及拉伸SHB实验,考察了应力状态三轴度、温度和应变率对材料失效应变的影响。由实验数据得到了Johnson-Cook失效模型参量,并通过出现失效的Taylor撞击实验和数值模拟进行了一定的验证,表明模型描述与实验结果的趋势一致。     

Time-dependent ratchetting experiments of SS304 stainless steel

The time-dependent strain cyclic characteristics and ratchetting behaviours of SS304 stainless steel were investigated by uniaxial/multiaxial cyclic loading tests at room and elevated temperatures (350 and 700 °C). The effects of loading rate, peak/valley strain or stress holds, ambient temperature and non-proportional loading path on the cyclic softening/hardening and ratchetting behaviours of the material were discussed. It is shown that: the cyclic deformation of the material presents remarkable time-dependence at room temperature and 700 °C; the cyclic hardening feature and ratchetting strain depend significantly on straining or stressing rate, hold-time, ambient temperature and the non-proportionality of loading path; the time-dependent ratchetting is resulted from the slight opening of hysteresis loop and visco-plasticity together, and the viscosity is a dominating factor at 700 °C; at 350 °C, abnormal rate-dependence and quick shakedown of ratchetting are observed due to the dynamic strain aging of the material at this temperature. Some significant conclusions are obtained, which are useful to construct a constitutive model to describe the time-dependent ratchetting behaviour of the material. It is also stated that the unified visco-plastic constitutive model discussed here cannot provide reasonable simulation to the time-dependent ratchetting at 700 °C, especially to that with certain peak/valley stress hold, since the effect of the high viscosity on time-dependent ratchetting cannot be properly described by using a unified visco-plastic flow rule.     

考虑塑性应变记忆恢复的循环棘轮本构模型研究

在统一粘塑性循环本构理论框架下,以Ohno-Abdel-Karim非线性随动硬化模型为基础,建立了一个循环本构模型。模型通过引入塑性应变幅值记忆效应,并在塑性应变记忆项中加入恢复系数,提高了对循环硬化材料单轴棘轮行为的预言能力。将模型应用于316L不锈钢单轴棘轮行为的描述中,模拟不同平均应力、应力幅值下的棘轮应变,均与实验数据吻合较好,证明本文改进的本构模型能合理地描述循环硬化材料的单轴棘轮行为。     

Crystal plasticity-based forming limit prediction for non-cubic metals: Application to Mg alloy AZ31B

A viscoplastic crystal plasticity model is incorporated within the Marciniak–Kuczynski (M–K) approach for forming limit curve prediction. The approach allows for the incorporation of crystallographic texture-induced anisotropy and the evolution of the same. The effects of mechanical twinning on the plastic response and texture evolution are also incorporated. Grain-level constitutive parameters describing the temperature dependent behavior of hexagonal close packed Mg alloy, AZ31B, sheets at discrete temperatures are used as a first application of the model. A trade-off between significant strain hardening behavior at lower temperatures (∼150 °C), and significant strain rate hardening at higher temperatures (∼200 °C) lead to similarities in the predicted forming limits. The actual formability of this alloy depends strongly on temperature within this range, and this distinction with the current modeling is related to more localized instability-based failure mechanisms at the lower temperatures than is assumed in the M–K approach. It is shown that the strain path dependence in the strain hardening response is significant and that it influences the forming limits in a predictable way. For broader applicability, a means of incorporating dynamic recrystallization into the crystal plasticity model is required.     

Modeling and validation of the large deformation inelastic response of amorphous polymers over a wide range of temperatures and strain rates

A robust physically consistent three-dimensional constitutive model is developed to describe the finite mechanical response of amorphous polymers over a wide range of temperatures and strain rates, including the rubbery region and for impact loading rates. This thermomechanical model is based on an elastic–viscoplastic rheological approach, wherein the effects of temperature, strain rate, and hydrostatic pressure are accounted for. Intramolecular, as well as intermolecular, interactions under large elastic–inelastic behavior are considered for the mechanisms of deformation and hardening. For a wide range of temperatures and strain rates, our simulated results for poly(methyl methacrylate) (PMMA) and polycarbonate (PC) are in good agreement with experimental observations.     

Uniaxial ratcheting of SS304 stainless steel at high temperatures: visco-plastic constitutive model

The uniaxial ratcheting of SS304 stainless steel at high temperatures (300, 600 and 700 °C) were analyzed experimentally, and described by a cyclic constitutive visco-plasticity model. The rate dependence of the material was accounted for by introducing a viscous term. The cyclic hardening and cyclic flow behavior of the material under asymmetrical stress-controlled cycling were described by the evolution rules of kinematic hardening back stress and isotropic deforming resistance. Under the isothermal condition, temperature effect was included by terms involving temperature in the evolution equations of isotropic deforming resistance. The effect of load history on ratcheting was also considered by introducing a fading memory function of the maximum inelastic strain amplitude and isotropic deformation resistance. After the material constants were determined from the experimental data, the uniaxial ratcheting of SS304 stainless steel was numerically simulated and compared with the corresponding experimental results at high temperatures. The predicted results agree well with the experimental ones.     

金属铍的压缩变形行为

利用材料试验机及Hopkinson杆装置系统开展热等静压金属铍在不同温度下的静动态压缩力学行为研究,获得了温度、应变率对金属铍屈服强度和加工硬化行为的影响规律。结果表明:金属铍在压缩应力状态下呈现出良好的塑性,同时其力学性能具有显著的应变率敏感性与热软化效应,屈服强度和流动应力随应变率提高呈明显增大趋势,随着温度升高逐渐降低。同时,室温下其加工硬化行为随着应变增大表现为分段硬化特征,随温度升高则趋于理想塑性。最后,采用修正的Johnson-Cook本构模型对实验结果进行了拟合,模型计算结果与实验结果吻合较好。     

Cyclic plasticity modeling with the distribution of non-linear relaxations approach

Motivated by the distribution of non-linear relaxation (DNLR) approach, a phenomenological model is proposed in order to describe the cyclic plasticity behavior of metals under proportional and non-proportional loading paths with strain-controlled conditions. Such a model is based on the generalization of the Gibbs's relationship outside the equilibrium of uniform system and the use of the fluctuation theory to analyze the material dissipation due to its internal reorganization. The non-linear cyclic stress–strain behavior of metals notably under complex loading is of particular interest in this study. Since the hardening effects are described appropriately and implicitly by the model, thus, a host of inelastic behavior of metals under uniaxial and multiaxial cyclic loading paths are successfully predicted such as, Bauschinger, strain memory effects as well as additional hardening. After calibrating the model parameters for two metallic materials, the model has demonstrated obviously its ability to describe the cyclic elastic-inelastic behavior of the nickel base alloy Waspaloy and the stainless steel 316L. The model is then implemented in a commercial finite element code simulating the cyclic stress–strain response of a thin-walled tube specimen. The numerical responses are in good agreement with experimental results.     

Experimental investigation and modeling of the tension behavior of polycarbonate with temperature effects from low to high strain rates

The effects of strain rate and temperature on the tension stress–strain responses of polycarbonate are experimentally investigated over a wide range of strain rates (0.001–1700 s⁻¹) and temperatures (0–120 °C). A modified split Hopkinson tension bar is used for high-rate uniaxial tension tests. Experimental results indicate that the stress–strain responses of polycarbonate at high strain rates exhibit the nonlinear characteristics including the obvious yielding and strain softening. The tension behavior is strongly dependent on the strain rate and temperature. The values of yield stress and strain at yield present a dramatic increase at higher strain rates and decrease with the increase in temperature. Moreover, there exists a significant rate-sensitivity transition in the polycarbonate tension yield behavior. Based on the experimental investigation, a physically based three-dimensional elastoplastic constitutive model for the finite deformation of glassy polymers is used to characterize the rate-temperature dependent yield and post-yield behavior of polycarbonate when subjected to tension loading. The model results are shown close to the experimental data within the investigated strain-rate and temperature ranges.     

受轴向冲击薄壁圆管的几何畸变相似律研究

对于受轴向冲击载荷作用的薄壁圆管动态响应的相似律问题,由于圆管的薄壁特性导致厚度无法与高度和半径按相同的比例进行结构缩放,从而产生模型的几何畸变,此时传统的相似律已无法描述原型与畸变模型之间的动态响应规律。基于薄壁圆管轴向冲击问题的控制方程,通过能量守恒和量纲分析,推导了考虑几何畸变条件下轴向冲击载荷作用的理想弹塑性薄壁圆管动态响应的相似律。通过在给定应变与应变率区间上建立比例模型预测的流动屈服应力与原型流动屈服应力的最佳逼近关系,将几何畸变相似律进一步推广至包含应变率和应变硬化的材料。通过数值方法验证了提出的几何畸变模型相似律的适用性。分析结果表明,提出的考虑厚度畸变的受轴向冲击薄壁圆管的相似律可用于预测原型结构的冲击动态响应,并显著降低比例模型与原型结构平均载荷和能量的偏差。     

Modelling large deformation behaviour under loading–unloading of semicrystalline polymers: Application to a high density polyethylene

In this work, the large deformation behaviour under monotonic loading and unloading of a high density polyethylene (HDPE) is studied. To analyze the nonlinear time-dependent response of the material, mechanical tests were conducted at room temperature under constant true strain rates and stress relaxation conditions. A physically-based inelastic model written under finite strain formulation is proposed to describe the mechanical behaviour of HDPE. In the model, the inelastic mechanisms involve two parallel elements: a visco-hyperelastic network resistance acting in parallel with a viscoelastic–viscoplastic intermolecular resistance where the amorphous and crystalline phases are explicitly taken into consideration. The semicrystalline polymer is considered as a two-phase composite. The influence of the crystallinity on the loading and unloading behaviour is investigated. Numerical results are compared to experimental data. It is shown that the model is able to accurately reproduce the experimental observations corresponding to monotonic loading, unloading and stress relaxation behaviours at different strain levels. Finally, the model capabilities to capture cyclic loading–unloading behaviour up to large strains are discussed. To demonstrate the improved modelling capabilities, simulations are also performed using the original model of Boyce et al. [Boyce, M.C., Socrate, S., Llana, P.G., 2000. Constitutive model for the finite deformation stress–strain behavior of poly(ethylene terephthalate) above the glass transition. Polymer 41, 2183–2201] modified by Ahzi et al. [Ahzi, S., Makradi, A., Gregory, R.V., Edie, D.D., 2003. Modeling of deformation behavior and strain-induced crystallization in poly(ethylene terephthalate) above the glass transition temperature. Mechanics of Materials 35, 1139–1148].     

Dislocation-mediated strain hardening in tungsten: Thermo-mechanical plasticity theory and experimental validation

A self-consistent thermo-mechanical model to study the strain-hardening behavior of polycrystalline tungsten was developed and validated by a dedicated experimental route. Dislocation–dislocation multiplication and storage, as well dislocation-grain boundary (GB) pinning were the major mechanisms underlying the evolution of plastic deformation, thus providing a link between the strain hardening behavior and material's microstructure. The microstructure of the polycrystalline tungsten samples has been thoroughly investigated by scanning and electron microscopy. The model was applied to compute stress–strain loading curves of commercial tungsten grades, in the as-received and as-annealed states, in the temperature range of 500–1000 °C. Fitting the model to the independent experimental results obtained using a single crystal and as-received polycrystalline tungsten, the model demonstrated its capability to predict the deformation behavior of as-annealed samples in a wide temperature range and applied strain. The relevance of the dislocation-mediated plasticity mechanisms used in the model have been validated using transmission electron microscopy examination of the samples deformed up to different amounts of strain. On the basis of the experimental validation, the limitations of the model are determined and discussed.     

Cyclic multiaxial and shear finite deformation responses of OFHC Cu. Part II: An extension to the KHL model and simulations

In this part, the Khan–Huang–Liang (KHL) constitutive model was extended to account for kinematic hardening characteristic behavior of materials. The extended model is then generalized and used to simulate experimental response of oxygen free high conductivity (OFHC) copper under cyclic shear straining and biaxial tension–torsion (multiaxial ratchetting) experiments presented in Part I (Khan et al., 2007). In addition, a new modification for the non-linear kinematic hardening rule of Karim–Ohno (Abdel-Karim and Ohno, 2000) is proposed to simulate multiaxial ratchetting behaviors. Although, the kinematic hardening contributes the most to the response, it is shown that, the loading rate effect, and a coupled isotropic and kinematic hardening effect should also be considered while simulating the multiaxial ratchetting behavior of OFHC copper. Furthermore, the newly modified kinematic hardening rules is able to fairly well simulate the multiaxial ratchetting experiments under different loading conditions, irrespective of the value of applied axial tensile stress, shear strain amplitude, pre-cyclic hardening and/or loading sequence.     

Multiaxial cyclic deformation and non-proportional hardening employing discriminating load paths

Some novel discriminating multiaxial cyclic strain paths with incremental and random sequences were used to investigate cyclic deformation behavior of materials with low and high sensitivity to non-proportional loadings. Tubular specimens made of 1050 QT steel with no non-proportional hardening and 304L stainless steel with significant non-proportional hardening were used. 1050 QT steel was found to exhibit very similar behavior under various multiaxial loading paths, whereas significant effects of loading sequence were observed for 304L stainless steel. In-phase cycles with a random sequence of axial-torsion cycles on an equivalent strain circle were found to cause cyclic hardening levels similar to 90° out-of-phase loading of 304L stainless steel. In contrast, straining with a small increment of axial-torsion on an equivalent strain circle results in higher stress than for in-phase loading of 304L stainless steel, but the level of hardening is lower than for 90° out-of-phase loading. Tanaka’s non-proportionality parameter coupled with a Armstrong–Fredrick incremental plasticity model, and Kanazawa et al.’s empirical formulation as a representative of such empirical models were used to predict the stabilized stress response of the two materials under variable amplitude axial-torsion strain paths. Consistent results between experimental observations and predictions were obtained by employing the Tanaka’s non-proportionality parameter. In contrast, the empirical model resulted in significant over-prediction of stresses for 304L stainless steel.     

Dynamic limit analysis formulation for impact simulation of structural members

This paper introduces an extended concept of limit analysis to deal with the dynamic equilibrium condition considering the inertia and strain-rate effect for dynamic behavior of structures. The conventional limit analysis method has been applied to only static collapse analysis of structures without consideration of dynamic effects in the structural behavior. A dynamic formulation for the limit analysis has been derived for incremental analysis dealing with time integration, strain and stress evaluation, strain hardening, strain-rate hardening and thermal softening. The time dependent term in the governing equation is integrated with the WBZ-α method. The dynamic material behavior is described by the Johnson–Cook model in order to consider strain-rate hardening and thermal softening as well as strain hardening. Simulations have been carried out for impact analysis of a Taylor bar and an S-rail and their numerical results are compared with elasto-plastic explicit analysis results by LS-DYNA3D. Comparison demonstrates that the dynamic finite element limit analysis can predict the crashworthiness of structural members effectively with less effort and computing time than the commercial code compared. The crashworthiness of a structure with the rate-dependent constitutive model is also compared to that with the quasi-static constitutive relation in order to investigate the dynamic effect on deformation of structures.