Formability of TWIP (twinning induced plasticity) automotive sheets

The main objective of this work is to experimentally and numerically evaluate the macro-performance of the automotive TWIP (twinning induced plasticity) sheet in conjunction with formability. In order to characterize the mechanical properties, the simple tension and compression tests were performed for anisotropic properties, while the strain rate test was carried out to evaluate strain rate sensitivity. The forming limit diagram was measured and incorporated into the simulation program, while the theoretical prediction of the diffuse and localized necking was also carried out utilizing Hill’s and the M-K theories as well as Dorn’s and Swift’s diffuse theories. Note that the generalized criteria of Hill’s, Dorn’s and Swift’s theories were derived for general anisotropic sheets as well in this work. For numerical simulations, the anisotropic yield functions Yld2000-2d and Hill48 as well as the isotropic Mises yield function were selectively applied along with the isotropic hardening law. Formability verification was performed, utilizing Yld2000-2d, for the hemispherical dome stretching, notch and simple tension tests with specimens selectively prepared by milling and punching, while anisotropic properties were verified through the three point bending and cylindrical cup drawing tests, comparing the performance of the three yield functions.     

Finite element and experimental method for analyzing the effects of martensite morphologies on the formability of DP steels

Abstract

In this article, we investigated the effect of martensite morphology on the mechanical properties and formability of dual phase steels. At first, three heat treatment cycles were subjected to a low-carbon steel to produce ferrite–martensite microstructure with martensite morphology of blocky-shaped, continuous, and fibrous. Tensile tests were then carried out so as to study mechanical properties, particularly the strength and strain hardening behavior of dual phase steels. In order to study the formability of dual phase samples, Forming Limit Diagram was obtained experimentally and numerically. Experimental forming limit diagram was obtained using Nakazima forming test, while Finite Element Method was utilized to numerically predict the forming limit diagram. The results indicated that the dual phase samples with fibrous martensite morphology had the highest tensile properties and strain rate hardening out of the three different microstructures. Blocky-shaped martensite morphology, on the other hand, had the worst mechanical properties. The study of the strain hardening behavior of dual phase sample by Kocks–Mecking-type plots, evinced two stages of strain hardening for all specimens with different microstructures: stages III and IV. The forming limit diagram of dual phase steels also proved that samples with fibrous martensite morphology had the best formability compared to other two microstructures. The simulated forming limit diagram manifested that there is a good agreement between experimental results and those obtained by FEM.     

应变硬化薄板夹持型滚弯成形过程分析

现存文献对薄板的连续滚弯成形过程的数学模拟多限于讨论三辊轮解锥型弯板过程,且基于较多的假定和板的小变形情况。而开发有更多的功能和更高的生产效率的四辊轮连续弯板机需研究夹持型滚弯成形过程。本文将「１」的工作,理想塑性材料平面应变薄板的夹持型连续滚弯成形过程的弹塑性大变形弯曲的力学分析,推广至一般应变硬化材料的情形。得出了控制微分方程的解析解。采用固支边界条件模拟刚度很大的平持系统对板的夹持。通过对线     

带接头管道在内压与循环载荷作用下的棘轮行为研究

液压管道在服役过程中受内压和循环弯曲载荷的共同作用．管道经常处于非比例循环加载状态,尤其是在管道接头位置处,容易产生棘轮行为,对管道的服役寿命有不利影响．因此,本文采用充液管道悬臂弯曲加载方式,对管道在接头位置处的棘轮响应进行研究。首先通过管材实验确定了材料的非线性等向/随动强化模型参数,并通过应变的实验测量结果与数值仿真结果的比较,验证了本构模型的有效性,然后建立了悬臂管道的有限元模型,模拟分析内压水平,内压小幅脉动,管道壁厚等因素对管道棘轮行为的影响．通过对带接头管道棘轮行为的研究分析,为进一步完善液压管道的设计,提高液压管道的可靠性,提供一定的理论基础．     

The shear fracture of dual-phase steel

Unexpected fractures at high-curvature die radii in sheet forming operations limit the adoption of advanced high strength steels (AHSS) that otherwise offer remarkable combinations of high strength and tensile ductility. Identified as “shear fractures” or “shear failures,” these often show little sign of through-thickness localization and are not predicted by standard industrial simulations and forming limit diagrams. To understand the origins of shear failure and improve its prediction, a new displacement-controlled draw-bending test was developed, carried out, and simulated using a coupled thermo-mechanical finite element model. The model incorporates 3D solid elements and a novel constitutive law taking into account the effects of strain, strain rate, and temperature on flow stress. The simulation results were compared with companion draw-bend tests for three grades of dual-phase (DP) steel over a range of process conditions. Shear failures were accurately predicted without resorting to damage mechanics, but less satisfactorily for DP 980 steel. Deformation-induced heating has a dominant effect on the occurrence of shear failure in these alloys because of the large energy dissipated and the sensitivity of strain hardening to temperature increases of the order of 75 °C. Isothermal simulations greatly overestimated the formability and the critical bending ratio for shear failures, thus accounting for the dominant effect leading to the inability of current industrial methods to predict forming performance accurately. Use of shell elements (similar to industrial practice) contributes to the prediction error, and fracture (as opposed to strain localization) contributes for higher-strength alloys, particularly for transverse direction tests. The results illustrate the pitfall of using low-rate, isothermal, small-curvature forming limit measurements and simulations to predict the failure of high-rate, quasi-adiabatic, large-curvature industrial forming operations of AHSS.     

Analysis of sheet metal formability through isotropic and kinematic hardening models

The present paper aims at analysing the sheet metal formability through several isotropic and kinematic hardening models. Specifically, a special attention is paid to the physically-based hardening model of Teodosiu and Hu (1995), which accounts for the anisotropic work-hardening induced by the microstructural evolution at large strains, as well as to some more conventional hardening models, including the isotropic Swift strain-hardening power law, and the Voce saturation strain-hardening law, combined with a non-linear kinematic hardening described by the Armstrong–Frederick law. The onset of localized necking is simulated by an advanced sheet metal forming limit model which connects, through the Marciniak–Kuczinsky analysis, the hardening models with the anisotropic yield criterion Yld2000-2d (Barlat et al., 2003). Both linear and complex strain paths are taken into account. The selected material is a DC06 steel sheet. The validity of each model is assessed by comparing the predicted forming limits with experimental results carefully obtained on this steel. The origin of discrepancy in the predicted results using different hardening models is thoroughly analyzed.     

A numerical method for rapid estimation of drawbead restraining force based on non-linear,anisotropic constitutive equations

Numerical procedures to predict drawbead restraining forces (DBRF) were developed based on the semi-analytical (non-finite-element) hybrid membrane/bending method. The section forces were derived by equating the work to pull sheet material through the drawbead to the work required to bend and unbend the sheet along with frictional forces on drawbead radii. As a semi-analytical method, the new approach was especially useful to analyze the effects of various constitutive parameters with less computational cost. The present model could accommodate general non-quadratic anisotropic yield function and non-linear anisotropic hardening under the plane strain condition. Several numerical sensitivity analyses for examining the effects of process parameters and material properties including the Bauschinger effect and the shape of yield surface on DBRF were presented. Finally, the DBRFs of SPCC steel sheet passing a single circular drawbead were predicted and compared with the measurements.     

泡沫铝本构行为研究进展

作为泡沫金属的典型代表, 泡沫铝是一种集结构性和功能性于一体的新型材料, 随着生产工艺的提高及国民经济的发展, 泡沫铝在航空航天、交通运输、建筑工程、机械制造等领域的应用日趋深入和广泛. 工程中复杂的载荷工况, 对泡沫铝本构模型的准确性及实用性提出了更高的要求. 国内外学者对泡沫铝进行了大量的试验研究及有限元数值分析工作, 随着各种研究和数值分析本构模型的提出、验证及修正等, 人们对泡沫铝的力学性能及其应用的认识逐渐提高. 本文首先简要概述了泡沫铝力学性能的试验研究及有限元数值分析的进展, 然后重点综述泡沫铝本构行为的国内外研究进展及现状. 最后, 针对目前存在的问题, 讨论并展望了泡沫铝本构行为的研究发展趋势. 完善本构模型表征所需的特征参量, 将各向异性或横观各项同性材料假定引入模型构建体系, 明确静水压缩响应及单轴压缩响应在材料强化中的权重, 建立例如随动强化般更加真实且准确反映泡沫铝强化过程的强化准则, 将应变率效应研究成果引入本构模型中等是现有泡沫铝本构模型体系中的重要研究方向.      

Spring-back evaluation of automotive sheets based on isotropic–kinematic hardening laws and non-quadratic anisotropic yield functions,part III: applications

In order to improve the prediction capability of spring-back in the computational analysis of automotive sheet forming processes, the modified Chaboche type combined isotropic–kinematic hardening law was formulated to account for the Bauschinger and transient behavior in Part I. As for the yield stress function, the non-quadratic anisotropic yield potential, Yld2000-2d, was utilized under the plane stress condition. Experimental procedures to obtain the material parameters of the combined hardening law and the yield potential were presented in Part II for three automotive sheets. For verification purposes, comparisons of simulations and experiments were performed here for the unconstrained cylindrical bending, the 2-D draw bending and the modified industrial part (the double-S rail). For all three applications, simulations showed good agreements with experiments. Simplified one-dimensional plane strain analytical and numerical methods were also developed here to better understand the spring-back in forming processes.     

A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming

In this paper an anisotropic material model based on non-associated flow rule and mixed isotropic–kinematic hardening was developed and implemented into a user-defined material (UMAT) subroutine for the commercial finite element code ABAQUS. Both yield function and plastic potential were defined in the form of Hill’s [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. R. Soc. Lond. A 193, 281–297] quadratic anisotropic function, where the coefficients for the yield function were determined from the yield stresses in different material orientations, and those of the plastic potential were determined from the r-values in different directions. Isotropic hardening follows a nonlinear behavior, generally in the power law form for most grades of steel and the exponential law form for aluminum alloys. Also, a kinematic hardening law was implemented to account for cyclic loading effects. The evolution of the backstress tensor was modeled based on the nonlinear kinematic hardening theory (Armstrong–Frederick formulation). Computational plasticity equations were then formulated by using a return-mapping algorithm to integrate the stress over each time increment. Either explicit or implicit time integration schemes can be used for this model. Finally, the implemented material model was utilized to simulate two sheet metal forming processes: the cup drawing of AA2090-T3, and the springback of the channel drawing of two sheet materials (DP600 and AA6022-T43). Experimental cyclic shear tests were carried out in order to determine the cyclic stress–strain behavior and the Bauschinger ratio. The in-plane anisotropy (r-value and yield stress directionalities) of these sheet materials was also compared with the results of numerical simulations using the non-associated model. These results showed that this non-associated, mixed hardening model significantly improves the prediction of earing in the cup drawing process and the prediction of springback in the sidewall of drawn channel sections, even when a simple quadratic constitutive model is used.     

Investigation of advanced strain-path dependent material models for sheet metal forming simulations

Sheet metal forming processes often involve complex loading sequences. To improve the prediction of some undesirable phenomena, such as springback, physical behavior models should be considered. This paper investigates springback behavior predicted by advanced elastoplastic hardening models which combine isotropic and kinematic hardening and take strain-path changes into account. A dislocation-based microstructural hardening model formulated from physical observations and the more classical cyclic model of Chaboche have been considered in this work. Numerical implementation was carried out in the ABAQUS software using a return mapping algorithm with a combined backward Euler and semi-analytical integration scheme of the constitutive equations. The capability of each model to reproduce transient hardening phenomena at abrupt strain-path changes has been shown via simulations of sequential rheological tests. A springback analysis of strip drawing tests was performed in order to emphasize the impact of several influential parameters, namely: process, numerical and behavior parameters. The effect of the two hardening models with respect to the process parameters has been specifically highlighted.     

Determination of material intrinsic length and strain gradient hardening in microbending process

Microbending experiments of pure aluminum show that the springback angles increase with the decrease of foil thickness, which indicates obvious size effects and attributes to plastic strain gradient hardening. Then a constitutive model, taking into accounts both plastic strain and plastic strain gradient hardening, is proposed to analyze the microbending process of thin foil. The model is based on the relationship between shear yield stress and dislocation density, in which the material intrinsic length is related to material properties and average grain numbers along the characteristic scale direction of part. It is adopted in analytical model to calculate the non-dimensional bending moment and predict the springback angle after microbending. It is confirmed that the predictions by the proposed hardening model agree well with the experimental data, while those predicted by the classical plasticity model cannot capture such size effects. The contribution of plastic strain gradient increases with the decrease of foil thickness and is independent on the bending angle.     

Characterization of mechanical properties by indentation tests and FE analysis – validation by application to a weld zone of DP590 steel

In this work, the mechanical properties of the metal active gas (MAG) weld zone and heat affected zone (HAZ) were characterized utilizing the continuous indentation method together with its finite element (FEM) analysis. To verify the measured properties, uni-axial tension and three point bending tests were performed for DP590 welded specimens. For numerical simulations, the isotropic hardening law was assumed along with the non-quadratic anisotropic yield function, Yld2000-2d. As for the failure criterion of the base material and weld zones particularly for the failure evaluation in the uni-axial tension test, Hill’s bifurcation theory and the M–K theory were applied to calculate the forming limit diagrams considering all measured properties including strain-rate sensitivity.     

Constitutive modeling for anisotropic/asymmetric hardening behavior of magnesium alloy sheets: Application to sheet springback

The constitutive model for the unusual asymmetric hardening behavior of magnesium alloy sheet presented in a companion paper (Lee, M.G., Wagoner, R.H., Lee, J.K., Chung, K., Kim, H.Y., 2008. Constitutive modeling for anisotropic/asymmetric hardening behavior of magnesium alloy sheet, Int. J. Plasticity 24(4), 545–582) was applied to the springback prediction in sheet metal forming. The implicit finite element program ABAQUS was utilized to implement the developed constitutive equations via user material subroutine. For the verification purpose, the springback of AZ31B magnesium alloy sheet was measured using the unconstrained cylindrical bending test of Numisheet (Numisheet ’2002 Benchmark Problem, 2002. In: Yang, D.Y., Oh, S.I., Huh, H., Kim, Y.H. (Eds.), Proceedings of 5th International Conference and Workshop on Numerical Simulation of 3D Sheet Forming Processes, Jeju, Korea) and 2D draw bend test. With the specially designed draw bend test the direct restraining force and long drawn distance were attainable, thus the measurement of the springback could be made with improved accuracy comparable with conventional U channel draw bend test. Besides the developed constitutive models, other models based on isotropic constitutive equations and the Chaboche type kinematic hardening model were also considered. Comparisons were made between simulated results by the finite element analysis and corresponding experiments and the newly proposed model showed enhanced prediction capability, which was also supported by the simple bending analysis adopting asymmetric stress–strain response.     

Stress integration method for a nonlinear kinematic/isotropic hardening model and its characterization based on polycrystal plasticity

Sheet metal forming processes generally involve non-proportional strain paths including springback, leading to the Bauschinger effect, transient hardening, and permanent softening behavior, that can be possibly modeled by kinematic hardening laws. In this work, a stress integration procedure based on the backward-Euler method was newly derived for a nonlinear combined isotropic/kinematic hardening model based on the two-yield’s surfaces approach. The backward-Euler method can be combined with general non-quadratic anisotropic yield functions and thus it can predict accurately the behavior of aluminum alloy sheets for sheet metal forming processes. In order to characterize the material coefficients, including the Bauschinger ratio for the kinematic hardening model, one element tension–compression simulations were newly tried based on a polycrystal plasticity approach, which compensates extensive tension and compression experiments. The developed model was applied for a springback prediction of the NUMISHEET’93 2D draw bend benchmark example.     

Modeling multi-axial deformation of planar anisotropic elasto-plastic materials,part I: Theory

In sheet metal forming processes local material points can experience multi-axial and multi-path loadings. Under such loading conditions, conventional phenomenological material formulations are not capable to predict the deformation behavior within satisfying accuracy. While micro-mechanical models have significantly improved the understanding of the deformation processes under such conditions, these models require large sets of material data to describe the micromechanical evolution and quite enormous computation expenses for industrial applications. To reduce the drawbacks of phenomenological material models under the multi-path loadings a new anisotropic elasto-plastic material formulation is suggested. The model enables the anisotropic yield surface to grow (isotropic hardening), translate (kinematic hardening) and rotate (rotation of the anisotropy axes) with respect to the deformation, while the shape of the yield surface remains essentially unchanged.Essentially, the model is formulated on the basis of an Armstrong–Frederick type kinematic hardening, the plastic spin theory for the reorientation of the symmetry axes of the anisotropic yield function, and additional terms coupling these expressions. The capability of the model is illustrated with multi-path loading simulations in ‘tension-shear’ and ‘reverse-shear’ to assess its performance with ‘cross’ hardening and ‘Bauschinger’ effects.     

Effects of initial anisotropy on the finite strain deformation behavior of glassy polymers

Solid phase deformation processing of glassy polymers produces highly anisotropic polymer components as a result of the massive reorientation of molecular chains during the large strain forming operation. Indeed, the polymer preform used as the starting materials is usually anisotropic owing to its prior deformation history. The process end product has often been fashioned for a particular application, i.e. to possess an increased flow strength along a particular axis, thereby exploiting the orientation induced anisotropy effects. The fully three-dimensional issues involved in the use of glassy polymer components include anisotropic flow strenghts, limiting extensibilities, and deformation patterns. These characteristics have been altered by the initial forming operation but are obviously not expected to be enhanced in all directions. The presence of anisotropy in structural components may also lead to premature failure or unexpected shear localization. In this report the effects of initial deformation and the associated anisotropies are investigated through uniaxial compression tests on preoriented polycarbonate (PC) and polymethylmethacrylate (PMMA) specimens. The evolving anisotropy is monitored by testing materials preoriented by various amounts of strain and under different states of deformation. The tensorial nature of the anisotropic material is characterized by examining the preoriented material response in three orthogonal directions. A model for the large strain deformation response of glassy polymers has been shown by Arruda and Boyce [in press] to be well predictive of the evolution of anisotropy during deformation in initially isotropic materials. Here the authors evaluate the ability of the model developed in Arruda and Boyce [in press] to predict several aspects of the anisotropic response of preoriented materials. Using material properties determined from the characterization of the isotropic material response and a knowledge of the anisotropic state of the preoriented material, model simulations are shown to accurately capture all aspects of the large strain anisotropic response including flow strengths, strain hardening characteristics, cross-sectional deformation patterns, and limiting extensibilities. Although anisotropy has been shown to evolve with temperature and strain rate in Boyce, Arruda and Jayachandran [in press] and also state of deformation in Arruda and Boyce [in press], we submit an experimental observation that the subsequent deformation response of preoriented polymers may be predicted using only a measure of optical anisotropy, and not the prior strain or thermal history. Optical anisotropy, as measured for example by birefringence, therefore represents a true internal variable indicative of the evolution of anisotropy with inelastic strain, state of strain, and temperature.     

Modeling multi-axial deformation of planar anisotropic elasto-plastic materials,part II: Applications

In order to predict the deformations under multi-axial and multi-path loadings in a phenomenological framework, a new rotational-isotropic-kinematic (RIK) hardening model has been suggested in the theory part of the paper combining isotropic, kinematic and rotational hardening. Essential features of this material model are Armstrong–Frederick type backstress components for kinematic hardening and a plastic spin for the rotational hardening describing the evolution of the symmetry axes of the anisotropic yield function.The purpose of this article is to illustrate the significance of the RIK hardening model in sheet metal forming applications as well as in springback predictions. With the rotational hardening and a correction term related to the kinematic hardening, the flow stress in each orientation can be described with few material parameters. Several benchmark problems are considered to illustrate and assess the performance of the RIK hardening model in comparison with other hardening models and experimental results.     

K0超固结结构性黏土的本构模型

K₀固结黏土在自然界广泛分布, 其通常同时具有超固结性与天然结构性, 而K₀超固结性又与K₀正常固结性质存在很大差异. 为了有效的描述K₀超固结性质, 在结构性模型基础上, 做了如下三点改进, 使得原模型拓展为同时考虑K₀超固结特性与天然结构性影响的本构模型. (1)引入相对应力比来描述屈服面, 并引入初始各向异性转轴参量ξ来表达初始各向异性对屈服面在p-q空间的位置影响. (2)基于给定的屈服面方程, 推导得到变相应力比参量, 并将变相应力比引入到统一硬化参数中, 利用统一硬化参数可以有效描述初始各向异性固结黏土在剪切加载下的剪缩与剪胀, 应变硬化及软化现象. (3)引入反映结构性胶结强度性质的胶结参量p_e, 并给出p_e随塑性偏应变的衰减演化方程, 利用胶结参量可描述结构性黏土的剪胀特性. 预测与试验结果对比表明, 所提的K₀超固结结构性模型可有效描述K₀超固结黏土的刚度提高效应, 黏土的包辛格效应, 结构性黏土胶结强度的损失现象以及结构性黏土的应变软化现象. 证明了所提模型的适用性以及合理性.      

A three-dimensional model of anisotropic hardening in metals and its application to the analysis of sheet metal formability

The hardening model proposed by Z. Mróz based on the uniaxial fatigue behavior of many metals is adopted to derive an incremental constitutive equation for general three-dimensional problems. This constitutive law is then employed in the analysis of metal forming problems to assess the influence of loading cycles, of the types involved in standard forming processes, on the ultimate formability of sheet metals. The predicted forming limit curves differ quantitatively from results obtained via an isotropie hardening model and differ qualitatively from those obtained via a kinematic model. Also investigated are the effects of such loading cycles on material response to simple tensile loading, which is often used to characterize a material. Significant differences between the present model and the other two models considered are observed in such characterizers of simple tensile behavior as the stress-strain curve, the anisotropy parameter and the uniform elongation. These differences suggest a rather simple experiment to identify the proper material model to be used in analyses of problems which involve loading cycles. Comparisons with some experimental results reveal that the employment of an anisotropic hardening model, such as the generalized Mróz model derived herein, is indeed crucial in accurately predicting material response to complicated loading histories.