Uniaxial and non-proportionally multiaxial ratcheting of SS304 stainless steel at room temperature: experiments and simulations

The uniaxial and non-proportionally multiaxial ratcheting behaviors of SS304 stainless steel at room temperature were initially researched by experiment and then were theoretically described by a cyclic constitutive model in the framework of unified visco-plasticity. The effects of cyclic stress amplitude, mean stress, and their histories on the ratcheting were experimentally investigated under uniaxial and different multiaxial loading paths. The shapes of non-proportional loading paths were linear, circular, elliptical and rhombic, respectively. In the constitutive model, the rate-dependent behavior of the material was reflected by a viscous term; the cyclic flow and cyclic hardening behaviors of the material under asymmetrical stress-controlled cycling were reflected by the evolution rules of kinematic hardening back stress and isotropic deforming resistance, respectively. The effect of loading history on the ratcheting was also considered by introducing two fading memorization functions for maximum inelastic strain amplitude and isotropic deformation resistance, respectively, into the constitutive model. The effect of multiaxial loading path on the ratcheting was reflected by a non-proportional factor defined in this work. The predicting ability of the developed model was proved to be good by comparing the simulations with corresponding experiments.     

On the experimental determination of yield surfaces and some results of annealed 304 stainless steel

Factors affecting the experimental determination of yield surfaces are discussed. They include the elastic moduli and the zero offset strain, the strain domain used to determine the yield stress, the probing path, and the strain rate of probing. To obtain yield surfaces consistently, it is necessary to account for these factors. The initial and subsequent yield surfaces of annealed AISI type 304 stainless steel have been experimentally determined in the axial-torsional stress space. Three loading paths have been studied. They are a pure axial path, a pure torsional path, and a proportional axial-torsional path. Each path includes loading, unloading, reloading, and the cyclically steady state.     

Time-dependent behavior of passive skeletal muscle

An isotropic three-dimensional nonlinear viscoelastic model is developed to simulate the time-dependent behavior of passive skeletal muscle. The development of the model is stimulated by experimental data that characterize the response during simple uniaxial stress cyclic loading and unloading. Of particular interest is the rate-dependent response, the recovery of muscle properties from the preconditioned to the unconditioned state and stress relaxation at constant stretch during loading and unloading. The model considers the material to be a composite of a nonlinear hyperelastic component in parallel with a nonlinear dissipative component. The strain energy and the corresponding stress measures are separated additively into hyperelastic and dissipative parts. In contrast to standard nonlinear inelastic models, here the dissipative component is modeled using an evolution equation that combines rate-independent and rate-dependent responses smoothly with no finite elastic range. Large deformation evolution equations for the distortional deformations in the elastic and in the dissipative component are presented. A robust, strongly objective numerical integration algorithm is used to model rate-dependent and rate-independent inelastic responses. The constitutive formulation is specialized to simulate the experimental data. The nonlinear viscoelastic model accurately represents the time-dependent passive response of skeletal muscle.     

Experimental Investigation on Mechanical Behavior of Coarse Marble Under Six Different Loading Paths

A series of triaxial compression experiments were preformed for the coarse marble samples under different loading paths by the rock mechanics servo-controlled testing system. Based on the experimental results of complete stress-strain curves, the influence of loading path on the strength and deformation failure behavior of coarse marble is made a detailed analysis. Three loading paths (Paths I–III) are put forward to confirm the strength parameters (cohesion and internal friction angle) of coarse marble in accordance with linear Mohr-Coulomb criterion. Compared among the strength parameters, two loading paths (i.e. Path II by stepping up the confining pressure and Path III by reducing the confining pressure after peak strength) are suggested to confirm the triaxial strengths of rock under different confining pressures by only one sample, which is very applicable for a kind of rock that has obvious plastic and ductile deformation behavior (e.g. marble, chalk, mudstone, etc.). In order to investigate re-fracture mechanical behavior of rock material, three loading paths (Paths IV–VI) are also put forward for flawed coarse marble. The peak strength and deformation failure mode of flawed coarse marble are found depending on the loading paths (Paths IV–VI). Under lower confining pressures, the peak strength and Young’s modulus of damage sample (compressed until post-peak stress under higher confining pressure) are all lower compared with that of flawed sample; moreover mechanical parameter of damage sample is lower for the larger compressed post-peak plastic deformation of coarse marble. However under higher confining pressures (e.g. σ ₃?=?30 MPa), the axial supporting capacity and elastic modulus of damage coarse marble (compressed until post-peak stress under lower confining pressure) is not related to the loading path, while the deformation modulus and peak strain of damage sample depend on the difference of initial confining pressure and post-peak plastic deformation. The friction among crystal grains determines the strength behavior of flawed coarse marble under various loading paths. In the end, the effect of loading path on failure mode of intact and flawed coarse marble is also investigated. The present research provides increased understanding of the fundamental nature of rock failure under different loading paths.     

A plasticity model for pressure-dependent anisotropic cellular solids

The initial and subsequent yield surfaces for an anisotropic and pressure-dependent 2D stochastic cellular material, which represents solid foams, are investigated under biaxial loading using finite element analysis. Scalar measures of stress and strain, namely characteristic stress and characteristic strain, are used to describe the constitutive response of cellular material along various stress paths. The coupling between loading path and strain hardening is then investigated in characteristic stress–strain domain. The nature of the flow rule that best describes the plastic flow of cellular solid is also investigated. An incremental plasticity framework is proposed to describe the pressure-dependent plastic flow of 2D stochastic cellular solids. The proposed plasticity framework adopts the anisotropic and pressure-dependent yield function recently introduced by Alkhader and Vural [Alkhader M., Vural M., 2009a. An energy-based anisotropic yield criterion for cellular solids and validation by biaxial FE simulations. J. Mech. Phys. Solids 57(5), 871–890]. It has been shown that the proposed yield function can be simply calibrated using elastic constants and flow stresses under uniaixal loading. Comparison of stress fields predicted by continuum plasticity model to the ones obtained from FE analysis shows good agreement for the range of loading paths and strains investigated.     

Erratum to: Experimental Investigation on Mechanical Behavior of Coarse Marble Under Six Different Loading Paths

A series of triaxial compression experiments were preformed for the coarse marble samples under different loading paths by the rock mechanics servo-controlled testing system. Based on the experimental results of complete stress-strain curves, the influence of loading path on the strength and deformation failure behavior of coarse marble is made a detailed analysis. Three loading paths (Paths I–III) are put forward to confirm the strength parameters (cohesion and internal friction angle) of coarse marble in accordance with linear Mohr-Coulomb criterion. Compared among the strength parameters, two loading paths (i.e. Path II by stepping up the confining pressure and Path III by reducing the confining pressure after peak strength) are suggested to confirm the triaxial strengths of rock under different confining pressures by only one sample, which is very applicable for a kind of rock that has obvious plastic and ductile deformation behavior (e.g. marble, chalk, mudstone, etc.). In order to investigate re-fracture mechanical behavior of rock material, three loading paths (Paths IV–VI) are also put forward for flawed coarse marble. The peak strength and deformation failure mode of flawed coarse marble are found depending on the loading paths (Paths IV–VI). Under lower confining pressures, the peak strength and Young’s modulus of damage sample (compressed until post-peak stress under higher confining pressure) are all lower compared with that of flawed sample; moreover mechanical parameter of damage sample is lower for the larger compressed post-peak plastic deformation of coarse marble. However under higher confining pressures (e.g. σ ₃ = 30 MPa), the axial supporting capacity and elastic modulus of damage coarse marble (compressed until post-peak stress under lower confining pressure) is not related to the loading path, while the deformation modulus and peak strain of damage sample depend on the difference of initial confining pressure and post-peak plastic deformation. The friction among crystal grains determines the strength behavior of flawed coarse marble under various loading paths. In the end, the effect of loading path on failure mode of intact and flawed coarse marble is also investigated. The present research provides increased understanding of the fundamental nature of rock failure under different loading paths.     

Inelastic behavior of an AS4/PEEK composite under combined transverse compression and shear. Part I: experiments

The nonlinear behavior in shear and transverse compression of unidirectional AS4/PEEK and their interaction are investigated experimentally. The composite is rate dependent even at room temperature and its rate exponent is similar to that of neat PEEK. The material is tested under pure shear, pure compression and under biaxial loading histories. The biaxial tests are performed in a custom facility on thin strips of the material. The facility allows freedom to choose the loading path in the biaxial stress and strain spaces of interest. Tests are performed for three biaxial loading paths. In the first, the specimen is sheared then compressed while the shear stress is held constant; in the second, the specimen is compressed then sheared while the compressive stress is held constant; and in the third, the specimen is loaded simultaneously by proportional amounts of compression and shear. It was found that the induced deformation is influenced significantly by the loading history followed. Also, initial loading in shear or compression has only a modest effect on subsequent loading of the other type. An unorthodox yielding behavior for the composite results from this lack of interaction. Finally, the stresses at failure are found to trace an elliptical path in the shear–transverse compression plane, but the failure stress state is not significantly affected by the loading path followed.     

Path-dependent failure of inflated aluminum tubes

Our recent investigation on the formability of Al alloy tubes under combined internal pressure and axial load is expanded by examining the effect of the loading path traced. A set of Al-6260-T4 tubes were loaded along orthogonal stress paths to failure and the results are compared to those of the corresponding radial paths. It is confirmed that failure strains are path-dependent, but also is demonstrated that failure stresses become path-dependent if the prestrain is significant. The experiments are simulated using the previously developed finite element models and the calibration of the Yld2000-2D [Barlat, F., Brem, J.C., Yoon, J.W., Chung, K., Dick, R.E., Lege, D.J., Pourboghrat, F., Choi, S.-H., Chu, E., 2003. Plane stress yield function for aluminum alloy sheets-part I: theory. Int. J. Plasticity 19, 1297--1319] anisotropic yield function shown in [Korkolis, Y.P., Kyriakides, S., 2008b. Inflation and burst of anisotropic aluminum tubes. Part II: an advanced yield function including deformation-induced anisotropy. Int. J. Plasticity 24, 1625–1637] to yield accurate predictions of rupture for nine radial paths. The models are shown to reproduce the path dependence of the failure stresses and strains quite well. A group of additional radial and corner paths are subsequently examined numerically to enrich the existing data on path-dependence of failure. It is again shown that the amount of plastic prestraining in either of the two directions influences the difference of the failure stresses and strains between the radial and the corner stress paths.     

An experimental study of inhomogeneous cyclic plastic deformation of 1045 steel under multiaxial cyclic loading

An experimental study was conducted on the inhomogeneous cyclic plastic deformation of 1045 steel under multiaxial cyclic loading. Thin-walled tubular specimens were used and small strain gages were bonded on the specimen surface to characterize the local deformation. The controlled loading paths included cyclic tension–compression, cyclic torsion, proportional axial-torsion, 90°-out-of-phase axial-torsion, and fully reversed torsion with a constant axial stress. The maximum stress in each experiment was lower than the lower yield stress of the material. It was found that the cyclic plastic deformation within the gage section of the specimen under multiaxial stress state followed the three-stage process that was observed from uniaxial loading, namely, incubation, propagation, and saturation. The plastic deformation was significantly inhomogeneous during the propagation stage, and the inhomogeneity continued through the saturation stage. The duration of each stage and the saturated strains were dependent on the cyclic stress amplitude and the loading path. Multiaxial stress state reduced the incubation stage. With identical equivalent stress magnitude, the nonproportional loading path resulted in the shortest incubation and propagation stages, and the saturated equivalent plastic strain magnitude was the smallest. Although the deformation over the gage section was inhomogeneous, the plastic deformation in a given local area was found to be practically isotropic.     

花岗岩体中应力波传播计算的动态本构关系

在花岗岩体的弹性区域,对实测径向质点速度波形运用Lagrangian分析方法,得到了球面应力波传播的加载速率相关的应力应变关系曲线。由速率相关本构关系所作的计算结果表明,它能较正确地描述岩石中球面应力波传播过程中所体现的主要特征,即应力波峰值衰减指数大于1和波形剖面的展宽。     

复杂加载下混凝土的弹塑性本构模型

混凝土材料在不同应力路径下或复杂加载条件下会表现出差异性显著的应力应变关系,在小幅循环加载条件下,其应力应变关系会表现出类似于弹性变形的滞回曲线.在不同应力水平下,混凝土的应力应变关系以及破坏特性都具有静水压力相关特点,即随着静水压力增大,各向异性强度特性弱化.此外,混凝土受压及受拉破坏机理不同,因而对应于混凝土硬化损伤亦有不同,即可分为受压硬化损伤,受拉硬化损伤及两者的混合硬化损伤类型.基于Hsieh模型,对该模型进行了三点改进.（1）针对小幅循环加载下混凝土无塑性变形的试验规律,而模型中在应力水平较低的循环加载条件下始终存在塑性变形的预测问题,采用在边界面模型框架下,设置了应力空间的弹性域,初始屈服面与后续临界状态屈服面几何相似的假定.（2）基于广义非线性强度准则将原模型采用变换应力方法将其推广为三维弹塑性本构模型,采用变换后模型可合理的考虑不同应力路径对于子午面以及偏平面上静水压力效应形成的影响,并避免了边界面应力点奇异问题.（3）分别对拉压两种加载损伤模式建议了相应的硬化参数表达式,可分别用于描述上述加载中产生的应变软化及强度退化行为.基于多种加载路径模拟表明:所建立的三维弹塑性本构模型可合理地用于描述混凝土的一般应力应变关系特性.      

Effect of nonlinear strain paths on forming limits under isotropic and anisotropic hardening

The path-dependence of the conventional Forming Limit Diagram (FLD) is an important issue for its applications in industry. Great efforts have been made to understand the nature of the path-dependence with both experimental and theoretical approaches, many of them attempting to find a path-independent way for the application of forming limits. In this paper, we focus on the nonlinear strain path effect on forming limit predictions using both isotropic and anisotropic hardening models. The Forming Limit Diagram (FLD), Forming Limit Stress Diagram (FLSD) and Forming Limit Effective Strain Diagram (epFLD) of sheet metals subject to linear and nonlinear strain paths are analyzed and compared using the Marciniak–Kuczynski approach. An anisotropic hardening model based on Yoshida and Uemori development is adopted in this study, and it is coupled with the traditional Hill’48 yield surface. This model is capable of describing the complex Bauschinger phenomenon after the material undergoes the reverse loading process such as the early re-yielding, work-hardening stagnation and permanent softening. Two different scenarios for the change of strain paths are also investigated. In the first scenario, the sheet material is initially loaded with a fixed strain increment ratio, unloaded to the free stress state, and then reloaded with a different strain increment ratio until the forming limit is reached. In the second scenario, the material does not undergo elastic unloading. Instead, the strain path is abruptly changed to a different strain increment ratio and the material undergoes continuous loading until the forming limit is reached. It is found that the work-hardening behavior after the pre-straining and the loading scenario plays an important role in the path dependent behavior of forming limits. Detailed analysis reveals that the M–K approach may have contributed to the significance of path-dependence observed in this study, especially at high pre-strain levels.     

A New Experimental Technique for the Multi-axial Testing of Advanced High Strength Steel Sheets

This paper deals with the development of a new experimental technique for the multi-axial testing of flat sheets and its application to advanced high strength steels. In close analogy with the traditional tension-torsion test for bulk materials, the sheet material is subject to combined tension and shear loading. Using a custom-made dual actuator hydraulic testing machine, combinations of normal and tangential loading are applied to the boundaries of a flat sheet metal specimen. The specimen shape is optimized to provide uniform stress and strain fields within its gage section. Finite element simulations are carried out to verify the approximate formulas for the shear and normal stress components at the specimen center. The corresponding strain fields are determined from digital image correlation. Two test series are performed on a TRIP-assisted steel sheet. The experimental results demonstrate that this new experimental technique can be used to investigate the large deformation behavior of advanced high strength steel sheets. The evolution of the yield surface of the TRIP700 steel is determined for both radial and non-proportional loading paths.     

砂土的应力路径本构模型

将微元应力路径线性逼近,转变成与其充分接近且易于计算应变的等平均应力微元和等应力比微元,计算任意加荷应力路径所产生的塑性应变,建立了双屈服面的砂土应力路径本构模型．模型体现了岩土塑性理论分量屈服和非关联流动法则的要求,在p,q平面内根据双线性的屈服线确定了加卸载准则．结合广义非线性强度理论采用变换应力三维化方法简单、合理地使模型实现三维化．通过试验数据的验证表明,砂土应力路径本构模型可以合理地描述各种应力路径下砂土的变形和强度特性。     

Study on dynamic constitutive relations for concrete with finite deformation

The differences between finite deformation and infinitesimal deformation are discussed. They are exercised on elasto-viscoplastic constitutive relations of concrete. Then, a rate-dependent mechanics model was presented on the basis of Ottosen‘ s fourparameter yield criterion, where different loading surface transferring laws were taken into account, when material was in hardening stage or in softening stage, respectively. The model is well established, so that it can be applied to simulate the response of concrete subject to impact loading. Green-Naghdi stress rate was introduced as objective stress rate. Appropriate hypothesis was postulated in accordance with many experimental results, which could reflect the mechanical behaviour of concrete with large deformation. Available thoughts as well as effective methods are also provided for the research on related engineering problems.     

Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline polymers

A micromechanically based constitutive model for the elasto-viscoplastic deformation and texture evolution of semi-crystalline polymers is developed. The model idealizes the microstructure to consist of an aggregate of two-phase layered composite inclusions. A new framework for the composite inclusion model is formulated to facilitate the use of finite deformation elasto-viscoplastic constitutive models for each constituent phase. The crystalline lamellae are modeled as anisotropic elastic with plastic flow occurring via crystallographic slip. The amorphous phase is modeled as isotropic elastic with plastic flow being a rate-dependent process with strain hardening resulting from molecular orientation. The volume-averaged deformation and stress within the inclusions are related to the macroscopic fields by a hybrid interaction model. The uniaxial compression of initially isotropic high density polyethylene (HDPE) is taken as a case study. The ability of the model to capture the elasto-plastic stress-strain behavior of HDPE during monotonic and cyclic loading, the evolution of anisotropy, and the effect of crystallinity on initial modulus, yield stress, post-yield behavior and unloading-reloading cycles are presented.     

A new crystal plasticity scheme for explicit time integration codes to simulate deformation in 3D microstructures: Effects of strain path,strain rate and thermal softening on localized deformation in the aluminum alloy 5754 during simple shear

The predominant deformation mode during material failure is shear. In this paper, a crystal plasticity scheme for explicit time integration codes is developed based on a forward Euler algorithm. The numerical model is incorporated in the UMAT subroutine for implementing rate-dependent crystal plasticity model in LS-DYNA/Explicit. The sheet is modeled as a face centered cubic (FCC) polycrystalline aggregate, and a finite element analysis based on rate-dependent crystal plasticity is implemented to analyze the effects of three different strain paths consisting predominantly of shear. Finite element meshes containing texture data are created with solid elements. The material model can incorporate information obtained from electron backscatter diffraction (EBSD) and apply crystal orientation to each element as well as account for texture evolution. Single elements or multiple elements are used to represent each grain within a microstructure. The three dimensional (3D) polycrystalline microstructure of the aluminum alloy AA5754 is modeled and subjected to three different strain rates for each strain path. The effects of strain paths, strain rates and thermal softening on the formation of localized deformation are investigated. Simulations show that strain path is the most dominant factor in localized deformation and texture evolution.     

On shakedown theory for elastic-plastic materials and extensions

The idea that an elastic-plastic structure under given loading history may shake down to some purely elastic state (and hence to a safe state) after a finite amount of initial plastic deformation, can apply to many sophisticated material models with possible allowable changes of additional material characteristics, as has been done in the literature. Despite some claims to the contrary, it is shown; however, that the shakedown theorems in a Melan-Koiter path-independent sense have been extended successfully only for certain elastic-plastic hardening materials of practical significance. Shakedown of kinematic hardening material is determined by the ultimate and initial yield stresses, not the generally plastic deformation history-dependent hardening curve between. The initial yield stress is no longer the convenient one (corresponding to the plastic deformation at the level of 0.2%) as in usual elastic-plastic analysis but to be related to the shakedown safety requirement of the structure and should be as small as the fatigue limit for arbitrary high-cycle loading. Though the ultimate yield strength is well defined in the standard monotonic loading experiment, it also should be reduced to the so-called “high-cycle ratchetting” stress for the path-independent shakedown analysis. A reduced simple form of the shakedown kinematic theorem without time integrals is conjectured for general practical uses. Application of the theorem is illustrated by examples for a hollow cylinder, sphere, and a clamped disk, under variable (including quasiperiodic dynamic) pressure.     

拉─扭复合加载下不锈钢的弹塑性本构关系─I．实验

讨论了１Ｃｒ１８Ｎｉ９Ｔｉ不锈钢薄壁圆管试件沿三段折线、不同曲率的圆形和椭圆应变路径承受拉－扭复合载荷的实验．在塑性应变空间中，观察加载路径的内蕴几何参数对应力矢量大小、方向影响的规律．结果表明：响应的延迟角、瞬时软化和重新强化性质与路径的内蕴几何学密切相关；Ｌｅｎｓｋｙ的“局部确定性”假设不完全符合事实；变形历史和应变分量相互间的耦合效应对响应存在显著的影响．初步的电镜实验表明，材料中的位错组态和塑性应变历史密切相关