A mixed hardening rule coupled with Hill48’ yielding function to predict the springback of sheet U-bending

A mixed hardening model has been used to model the Bauschinger effect. This hardening model is based on Lemaitre and Chaboche nonlinear kinematic hardening theory to consider cyclic behavior and the Bauschinger effect. Hill’48 yielding criterion is used because of the general stress state and relative ease of formulation. The backward Euler return mapping algorithm is applied to calculate the stress and strain increment. The mixed hardening model is implemented based on UMAT subroutine of FEA code ABAQUS. The NUMISHEET’93 benchmark shows that the mixed hardening model coupled with anisotropic yield criteria can give a favorable springback angle prediction.     

Isotropic‐kinematic hardening framework to assess ratcheting response of steel samples undergoing asymmetric loading cycles

The present study intends to characterize ratcheting response of several steel alloys subject to asymmetric loading cycles through coupling the Ahmadzadeh‐Varvani kinematic hardening rule with isotropic hardening rules of Lee and Zavrel, Chaboche, and Kang. The Ahmadzadeh‐Varvani kinematic hardening rule was developed to address ratcheting progress over asymmetric stress cycles with relatively a simple framework and less number of coefficients. Inclusion of isotropic hardening rules to the framework improved ratcheting response of materials mainly over the first stage of ratcheting. Lee and Zavrel model (ISO‐I) developed an exponential function to account for accumulated plastic strain as yield surface is expanded over stage I and early stage II of ratcheting. Isotropic models by Chaboche (ISO‐II) and Kang (ISO‐III) encountered yield surface evolution in the framework by introducing an internal variable that takes into account the prior maximum plastic strain range. The choice of isotropic hardening model coupled to the kinematic hardening model is highly influenced by material softening/hardening response.     

An analytical model for predicting springback and side wall curl of sheet after U-bending

An analytical model for predicting sheet springback after U-bending is proposed in this paper based on Hill48 yielding criterion and plane strain condition. The model takes into account of the effects of deformation history, thickness thinning and neutral surface shift on the sheet springback of U-bending. Three rules for material hardening – kinematic, isotropic and combined hardening – have been used to consider the effect of complex deformation history that has undergone stretching, bending, and unbending deformations on the sheet springback. The model is applied to the benchmark of NUMISHEET’93 2-D draw bending problem. It indicates that the springback is overestimated when isotropic hardening is applied, while is underestimated when kinematic hardening is applied. For reverse loading problem, the combined hardening is a good approach to the practical material. In addition to that, the effects of blank holding force, friction coefficient between sheet and tools, sheet thickness and anisotropy have been investigated. When the shifting distance of neutral surface exceeds one-fourth of sheet thickness, the springback can be reduced effectively by increasing the blank holding force and friction between sheet and die. And the springback increases with anisotropy and friction between sheet and punch, and decreases with the sheet thickness.     

Experimental and numerical analyses of mean stress relaxation in cold expanded plate of Al‐alloy 2024‐T3 in double shear lap joints

In this paper, the mean stress relaxation behavior of simple Al‐alloy 2024‐T3 specimens and also the mean stress relaxation around the hole of cold expanded specimen are studied. The analyses are performed through the combination of the nonlinear isotropic hardening and Chaboche nonlinear kinematic hardening model accompanied by the results of experimental tests. The strain‐controlled axial tests are performed at two different strain amplitudes, while the stress‐controlled tests of cold expanded specimens are performed for three different imposed load amplitudes. The constitutive equations of the hardening model are coded as a UMAT subroutine in FORTRAN programming language and implemented in the commercial finite element code of ABAQUS. The accuracy of the hardening model has been proved in two steps: first by simulations of mean stress relaxation during the uniaxial strain‐controlled cyclic tests and second by simulation of strain ratcheting during the stress‐controlled cyclic loading. The stress and strain distributions after cold expansion process are examined as well as the mean stress relaxation due to cyclic loading. The results show the influences of imposed stress amplitude on increasing mean stress relaxation and also the effect of cold expansion level on reducing the mean stress relaxation.     

A simple integration algorithm for a non‐associated anisotropic plasticity model for sheet metal forming

In this paper, an anisotropic material model based on a non‐associated flow rule and nonlinear mixed isotropic‐kinematic hardening is developed. The quadratic Hill48 yield criterion is considered in the non‐associated model for both yield function and plastic potential to account for anisotropic behavior. The developed model is integrated based on fully implicit backward Euler's method. The resulting problem is reduced to only two simple scalar equations. The consistent local tangent modulus is obtained by exact linearization of the algorithm. All numerical development was implemented into user‐defined material subroutine for the commercial finite element code ABAQUS/Standard. The performance of the present algorithm is demonstrated by numerical examples. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulations of stress distributions in crankshaft sections under fillet rolling and bending fatigue tests

The residual stresses due to fillet rolling and the bending stresses near the fillets of crankshaft sections under bending fatigue tests are important driving forces to determine the bending fatigue limits of crankshafts. In this paper, the residual stresses and the bending stresses near the fillet of a crankshaft section under fillet rolling and subsequent bending fatigue tests are investigated by a two-dimensional plane strain finite element analysis based on the anisotropic hardening rule of Choi and Pan [Choi KS, Pan J. A generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials (in preparation)]. The evolution equation for the active yield surface during the unloading/reloading process is first presented based on the anisotropic hardening rule of Choi and Pan (in preparation) and the Mises yield function. The tangent modulus procedure of Peirce et al. [Peirce D, Shih CF, Needleman A. A tangent modulus method for rate dependent solids. Comput Struct 1984;18:875–87] for rate-sensitive materials is adopted to derive the constitutive relation. A user material subroutine based on the anisotropic hardening rule and the constitutive relation was written and implemented into ABAQUS. Computations were first conducted for a simple plane strain finite element model under uniaxial monotonic and cyclic loading conditions based on the anisotropic hardening rule, the isotropic and nonlinear kinematic hardening rules of ABAQUS. The results indicate that the plastic response of the material follows the intended input stress–strain data for the anisotropic hardening rule whereas the plastic response depends upon the input strain ranges of the stress–strain data for the nonlinear kinematic hardening rule. Then, a two-dimensional plane-strain finite element analysis of a crankshaft section under fillet rolling and subsequent bending was conducted based on the anisotropic hardening rule of Choi and Pan (in preparation) and the nonlinear kinematic hardening rule of ABAQUS. In general, the trends of the stress distributions based on the two hardening rules are quite similar after the release of roller and under bending. However, the compressive hoop stress based on the anisotropic hardening rule is larger than that based on the nonlinear kinematic hardening rule within the depth of 2 mm from the fillet surface under bending with consideration of the residual stresses of fillet rolling. The critical locations for fatigue crack initiation according to the stress distributions based on the anisotropic hardening rule appear to agree with the experimental observations in bending fatigue tests of crankshaft sections.     

Simulation of low cycle fatigue stress‐strain response in 316LN stainless steel using non‐linear isotropic kinematic hardening model—A comparison of different approaches

Cyclic stress‐strain response of 316LN stainless steel subjected to low cycle fatigue at strain amplitude of ±0.4% and at 873 K is simulated using finite element analysis with non‐linear isotropic‐kinematic hardening Chaboche model. Four different approaches have been used in simulating cyclic stress response and hysteresis loops: 3 based on Chaboche model‐parameters and the fourth on direct experimental data (stabilized loop and cyclic stress‐strain curve [CSSC]). Among them, simulations performed with direct experimental data have not yielded expected initial cyclic response. The source of data used for evaluation of kinematic‐hardening (KH) parameters determined the extent of closeness between experimental results and Chaboche‐model predictions. KH parameters determined from first‐cycle loop and modified‐CSSC predicted the overall stress‐strain response (from initial to stabilized condition) with reasonable fit, compared with other approaches. All 4 approaches though predicted stabilized response, simulations based on “KH‐parameters from stabilized‐cycle” accurately described stabilized response with coefficient of determination (r²) 0.995.     

Effect of hardening models on different ductile fracture criteria in sheet metal forming

Prediction of the fracture is one of the challenging issues which gains attention in sheet metal forming as numerical analyses are being extensively used to simulate the process. To have better results in predicting the sheet metal fracture, appropriate ductile fracture criterion (DFC), yield criterion and hardening rule should be chosen. In this study, the effects of different hardening models namely isotropic, kinematic and combined hardening rules on the various uncoupled ductile fracture criteria are investigated using experimental and numerical methods. Five different ductile fracture criteria are implemented to a finite element code by the user subroutines. The criterion constants of DFCs are obtained by the related experimental tests. The in-plane principle strains obtained by the finite element analyses for different DFCs are compared with the experimental results. Also, the experimental results are used to evaluate the principle strain values calculated by the finite element analysis for different combinations of DFCs and hardening rules. It is shown that some DFCs give better predictions if the appropriate hardening model is employed.     

Twist springback characteristics of dual-phase steel sheet after non-axisymmetric deep drawing

This paper aims to investigate the twist springback characteristics of advanced high strength steel sheet subjected to deep drawing. A C-rail benchmark, which leads to a particularly pronounced twist springback characteristics, was developed. For an accurate numerical modeling of the process, a non-quadratic anisotropic yield criterion integrated with combined isotropic and kinematic hardening model was used to describe the strain-stress behavior including anisotropy and Bauschinger effects. The corresponding mechanical experiments, namely uniaxial tension and forward-reverse simple shear tests were performed to determine the material parameters. The digital image correlation technique was applied for component tests as well as the deformation and stress-strain analysis. The experimental validation of the elastic-plastic finite element model was assessed by comparing maximum in-plane strain, thickness reduction distribution and twist springback of the drawn rail. To explore the source of twist springback, the deformation associated with in-plane stress and bending moment was analyzed. The results indicate that the bending moment before springback caused by non-symmetric stress states play an important role in twist springback and control. Certain regions of the die radius were varied in a numerical analysis to control the bending moment for the minimization of twist springback as well as the preliminary results of the relationship between the ratio of variable die radius and twist springback.     

屈服准则和硬化模型对5052铝板回弹仿真的影响

为选择适合于5052铝合金回弹仿真的材料模型,对LS-Dyna软件中4个材料模型MAT_36、MAT_122、MAT_125和MAT_226所采用的屈服准则和硬化模型进行了分析,采用这4个模型对5052铝板U形件的回弹进行了仿真,对回弹过程中圆角区的应力释放进行了讨论.同时,进行了U形件的回弹试验,并与仿真结果进行了比较.结果表明,4个材料模型中,基于Yoshida-Uemori随动硬化模型和Barlat’89屈服准则的材料模型MAT_226具有最好的回弹预测精度,由各向同性硬化模型和Hill’48屈服准则组合的材料模型MAT_122的回弹预测结果与试验结果的偏差最大.硬化模型对回弹预测精度的影响大于屈服准则的影响.     

A finite element analysis of plane strain dynamic crack growth in materials displaying the Bauschinger effect

In this paper dynamic crack growth in an elastic-plastic material is analyzed under mode I plane strain small-scale yielding conditions using a finite element procedure. The main objective of this paper is to investigate the influence of anisotropic strain hardening on the material resistance to rapid crack growth. To this end, materials that obey an incremental plasticity theory with linear isotropic or kinematic hardening are considered. A detailed study of the near-tip stress and deformation fields is conducted for various crack speeds. The results demonstrate that kinematic hardening does not oppose the role of inertia in decreasing the plastic strains and stresses near the crack tip with increase in crack speed to the same extent as isotropic strain hardening. A ductile crack growth criterion based on the attainment of a critical crack opening displacement at a small micro-structural distance behind the tip is used to obtain the dependence of the theoretical dynamic fracture toughness with crack speed. It is found that for any given level of strain hardening, the dynamic fracture toughness displays a much more steep increase with crack speed over the quasi-static toughness for the kinematic hardening material as compared to the isotropic hardening case.     

Cyclic uniaxial and multiaxial loading with yield surface distortion consideration on prediction of ratcheting

In this study, the yield surface distortion was incorporated in the cyclic plasticity modeling as well as its center movement. The combination of Chaboche’s model and the yield surface distortion model of Baltov was used in a set of uniaxial and multiaxial loadings. The variation of the stress amplitude and the mean stress and different multiaxial loadings such as tension-torsion tests were studied. It was shown that the consideration of the distortion of the yield surface via the distortion parameter and its sign in modeling has an important effect on the plastic strain increment determination and so on the ratcheting rate. The combined model was applied to the experimental results. It was shown that the combination of the nonlinear kinematic hardening model of Chaboche and the yield surface distortion leads to a good estimation of the ratcheting strain increment in different uniaxial and multiaxial tests.     

Study of hardening and cyclic plastic behaviour around the crack tip of C(T) specimen made by as‐received and annealed copper

The purpose of this research is determining experimentally the characteristics of tension and cyclic plastic behaviours of as‐received and annealed coppers and studying distribution of stress/strain field near the crack tip. Samples made by pure copper were annealed at 420°C for 40 minutes in electric furnace. To determine the properties of the cyclic plastic behaviour, proper tests with symmetric strain‐controlled conditions were performed on standard samples. Chaboche nonlinear hardening model was used to determine the cyclic plastic behaviour of both materials. According to results, annealing process creates isotropic hardening in the copper and also changes its initial kinematic hardening behaviour. Effects of the annealing and hardening on the variations of the stresses and strains around the crack tip were investigated. Also, ratcheting and mean stress relaxations versus number of cycles, inside the plastic region, were studied.     

Cyclic loading of beams based on the Chaboche model

The time independent non-unified version of the Chaboche constitutive model for the cyclic loading which includes the kinematic and isotropic hardening is discussed in detail. The performance of the Chaboche constitutive model in predicting ratcheting response of CS1026 steel for a broad set of mechanical uniaxial and biaxial loading histories is considered. A numerical iterative method is used to calculate the stresses and strains in beams due to cyclic loading. The reported experimental data of the stainless steel, available in the literature, is used for the verification of the results. It is concluded that the Chaboche model performs quite well in predicting the uniaxial ratcheting or shakedown responses. In addition, imposing the isotropic hardening effect to the constitutive equations results to lower ratcheting rate at initial cycles. While the kinematic hardening effect remains the major factor in prediction of the ratcheting response.     

Constitutive Modeling and Algorithmic Implementation of a Plasticity-like Model for Trabecular Bone Structures

Trabecular bone is a highly porous orthotropic cellular solid material present inside human bones such as the femur (hip bone) and vertebra (spine). In this study, an infinitesimal plasticity-like model with isotropic/kinematic hardening is developed to describe yielding of trabecular bone at the continuum level. One of the unique features of this formulation is the development of the plasticity-like model in strain space for a yield envelope expressed in terms of principal strains having asymmetric yield behavior. An implicit return-mapping approach is adopted to obtain a symmetric algorithmic tangent modulus and a step-by-step procedure of algorithmic implementation is derived. To investigate the performance of this approach in a full-scale finite element simulation, the model is implemented in a non-linear finite element analysis program and several test problems including the simulation of loading of the human femur structures are analyzed. The results show good agreement with the experimental data.     

Predicting damage and failure under low cycle fatigue in a 9Cr steel

Experimental data have been generated and finite element models developed to examine the low cycle fatigue (LCF) life of a 9Cr (FB2) steel. A novel approach, employing a local ductile damage initiation and failure model, using the hysteresis total stress–strain energy concept combined with element removal, has been employed to predict the failure in the experimental tests. The 9Cr steel was found to exhibit both cyclic softening and nonlinear kinematic hardening behaviour. The finite element analysis of the material's cyclic loading was based on a nonlinear kinematic hardening criterion using the Chaboche constitutive equations. The models’ parameters were calibrated using the experimental test data available. The cyclic softening model in conjunction with the progressive damage evolution model successfully predicted the deformation behaviour and failure times of the experimental tests for the 9Cr steels performed.     

Application of a semianalytical strain assessment and multiaxial fatigue analysis to compare rolling contact fatigue in twin-disk and full-scale wheel/rail contact conditions

A semianalytical model is introduced to assess rolling contact fatigue problems in railway applications. The constitutive law is based on the nonlinear kinematic and isotropic hardening model of Chaboche–Lemaitre, which allows the cyclic elastoplastic strain under the contact surface to be evaluated. The much higher computational effectiveness in comparison with finite element (FE) analyses is quantified. The Dang Van multiaxial fatigue criterion is implemented to evaluate the rolling contact fatigue in the subsurface elastic region where cracking is relatively rare but more dangerous than surface cracks. The influence of the presence of sulfides in the wheel matrix in decreasing fatigue strength is also assessed by means of Murakami's approach. The model is used to compare conditions under small-scale twin-disk experiments to full-scale wheel/rail contact conditions. It is found that, for the same Hertzian pressure, the small-scale contact is more conservative in that it causes a deeper plasticized layer as compared with the elliptical full-scale contact. In the investigated cases, crack initiation is also not expected according to Dang Van criterion in neither of the studied contact conditions.     

Service-type creep-fatigue experiments with cruciform specimens and modelling of deformation

Advanced material models for the application to component life prediction require multiaxial experiments. A biaxial testing system for cruciform test pieces has been established in order to provide data for creep, creep-fatigue and thermomechanical fatigue (TMF) experiments. For this purpose a cruciform specimen was developed with the aid of Finite element calculation and the specimen design was optimised for tension and compression load. The testing system is suitable for strain (displacement) and load control mode. A key feature deals with the opportunity to perform thermomechanical experiments. Further, a constitutive material model is introduced which is implemented as a user subroutine for Finite element applications. The constitutive material model of type Chaboche considers both isotropic as well as kinematic hardening and isotropic damage. Identification of material parameters is achieved by a combination of Neural networks and subsequent Nelder–Mead Method.