Determination of Anisotropic Plastic Constitutive Parameters Using the Virtual Fields Method

The aim of the present study is to retrieve all the anisotropic plastic constitutive parameters from uniaxial loading. A complex geometry which can provide very heterogeneous stress states in a uniaxial tensile test was chosen for steel sheet specimens. A digital image correlation technique was used for the full-field heterogeneous strain measurement. The orthotropic Hill1948 yield criterion with Swift isotropic hardening was adopted as an elasto-plastic constitutive model. The virtual fields method (VFM) was employed as an inverse analytical tool to determine the constitutive parameters. All the parameters were successfully identified using the VFM by combining two tensile test results obtained in rolling and transverse directions.     

Identification of Post-Necking Hardening Phenomena in Ductile Sheet Metal

A combined theoretical/experimental approach accurately quantifying post-necking hardening phenomena in ductile sheet materials that initially exhibit diffuse necking in tension is presented. The method is based on the minimization of the discrepancy between the internal and the external work in the necking zone during a quasi-static tensile test. The main focus of this paper is on the experimental validation of the method using an independent material test. For this purpose, the uniaxial tube expansion test is used to obtain uniaxial strain hardening behavior beyond the point of maximum uniform strain in a tensile test. The proposed method is used to identify the post-necking hardening behavior of a cold rolled interstitial-free steel sheet. It is demonstrated that commonly adopted phenomenological hardening laws cannot accurately describe all hardening stages. An alternative phenomenological hardening model is presented which enables to disentangle pre- and post-necking hardening behavior. Additionally, the influence of the yield surface on the identified post-necking hardening behavior is scrutinized. The results of the proposed method are compared with the hydraulic bulge test. Unlike the hydraulic bulge test, the proposed method predicts a decreased hardening rate in the post-necking regime which might be associated with probing stage IV hardening. While inconclusive, the discrepancy with the hydraulic bulge test suggests differential work hardening at large plastic strains.     

On the identification of a high-resolution multi-linear post-necking strain hardening model

The Finite Element Model Updating (FEMU) technique is an inverse method that enables to arrive at a complete solution to the problem of diffuse necking of a thick tensile specimen. Conventionally, FEMU relies on the identification of a phenomenological strain hardening law that inherently limits the accuracy of the method due to the predefined character of the adopted strain hardening law. A high-resolution multi-linear post-necking strain hardening model enables to describe more generically the actual strain hardening behaviour. A numerical concept study is used to scrutinise the identification of such a model using FEMU. It is shown that, unlike progressive identification strategies, a global identification strategy followed by a smoothing operation based on area conservation yields sufficiently accurate results. To study the experimental feasibility, the latter strategy is used to identify the post-necking strain hardening behaviour of a thick S690QL high-strength steel. To this purpose, a notched tensile specimen was loaded up to fracture, while the elongation was measured using Digital Image Correlation (DIC). It is shown that the global identification strategy suffers from experimental noise associated with DIC and the load signal.     

Nonproportional loading effects on elastic-plastic behavior based on stress resultants for thin plates of strain hardening materials

A stress resultant constitutive law in rate form is constructed for power-law hardening materials. The change of plate thickness is considered in the constitutive law. The elastic-plastic behavior of a plate element based on the stress resultant constitutive law under uniaxial combined tension and bending is determined under a limited number of nonproportional and unloading paths. The results based on the stress resultant constitutive law and the through-the-thickness integration method are compared within the context of both the small-strain and finite deformation approaches. The results indicate that the selection of the normalized equivalent stress resultant and the corresponding work-conjugate normalized equivalent generalized strain is appropriate for describing the hardening behavior in the stress resultant space. However, the hardening rule in a power law form must be modified for low hardening materials at large plastic deformation when finite deformation effects are considered.     

Inferring Post-Necking Strain Hardening Behavior of Sheets by a Combination of Continuous Bending Under Tension Testing and Finite Element Modeling

Experimental Mechanics - This paper presents a combined experimental and simulation approach to identify post-necking hardening behavior of ductile sheet metal. The method is based on matching a...     

一般加载规律的弹塑性本构关系

将有关文献给出一般加载规律一维全量理论的简单模型推广到一般加载规律的一维增量理论,进而推广到一般加载规律的多维增量理论,在此基础上,建立了推导一般加载规律的多维增量理论的本构关系的一种途径。应用这种途径,从应力空间的加载函数和应变空间的加载函数出发,推导了等向强化材料和被加热的等向强化材料的一般加载规律的弹塑性本构关系的两种表示形式。理论和实例均表明,这种途径对等向强化材料、随动强化材料和理想弹塑性材料均适用。     

Anisotropic plasticity model coupled with Lode angle dependent strain-induced transformation kinetics law

A phenomenological macroscopic plasticity model is developed for steels that exhibit strain-induced austenite-to-martensite transformation. The model makes use of a stress-state dependent transformation kinetics law that accounts for both the effects of the stress triaxiality and the Lode angle on the rate of transformation. The macroscopic strain hardening is due to nonlinear kinematic hardening as well as isotropic hardening. The latter contribution is assumed to depend on the dislocation density as well as the current martensite volume fraction. The constitutive equations are embedded in the framework of finite strain isothermal rate-independent anisotropic plasticity. Experimental data for an anisotropic austenitic stainless steel 301LN is presented for uniaxial tension, uniaxial compression, transverse plane strain tension and pure shear. The model parameters are identified using a combined analytical–numerical approach. Numerical simulations are performed of all calibration experiments and excellent agreement is observed. Moreover, we make use of experimental data from ten combined tension and shear experiments to validate the proposed constitutive model. In addition, punch and notched tension tests are performed to evaluate the model performance in structural applications with heterogeneous stress and strain fields.     

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.     

The influence of particle shape,volume fraction and distribution on post-necking deformation and fracture in uniaxial tension of AA5754 sheet materials

Finite element analysis was performed over a small particle field, edge constraint plane strain post-necking model. The aim is to understand the roles of particle shape, volume fraction and distribution over the post-necking deformation and fracture of AA5754-O sheet materials. For models containing one single particle, the post-necking deformation decreases when the particle varies from circular to elliptical. The inter-particle spacing, the major parameter of distribution to determine whether a pair of particles belongs to a stringer or not, was varied for models with two particles of circular or elliptical shape. The general trend is that the post-necking deformation and fracture strains decrease with decreasing spacing between particles. There is considerable difference in terms of both fracture topographies and strains for models containing 16 particles when distributions varied from random/uniform to stringer distributions. The post-necking deformation and fracture strains monotonically decrease with particle volume fractions for models with 4–64 particles of random or stringer distribution. This indicates that the post-necking behavior for AA5754-O alloys where the matrix material is rather ductile is not solely controlled by a single or pair of particles although they may become initiation places of damage. Multiple damaging sources such as stringers or large particles can act cooperatively and speed up the damaging propagation of the material, and therefore produce small post-necking deformation and early fracture. The center clustering of particles can be beneficial for post-necking behavior and bendability of sheet materials.     

Identification of Mechanical Material Behavior Through Inverse Modeling and DIC

Inverse methods offer a powerful tool for the identification of the elasto-plastic material parameters. One of the advantages with respect to classical material testing is the fact that those inverse methods are able to deal with heterogeneous deformation fields. The basic principle of the inverse method that is presented in this paper, is the comparison between experimentally measured strain fields and those computed by the finite element (FE) method. The unknown material parameters in the FE model are iteratively tuned so as to match the experimentally measured and the numerically computed strain fields as closely as possible. This paper describes the application of an inverse method for the identification of the hardening behavior and the yield locus of DC06 steel, based on a biaxial tensile test on a perforated cruciform specimen. The hardening behavior is described by a Swift type hardening law and the yield locus is modeled with a Hill 1948 yield surface.     

Spring-back evaluation of automotive sheets based on isotropic-kinematic hardening laws and non-quadratic anisotropic yield functions: Part II: characterization of material properties

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 are presented here in Part II for three automotive sheets: AA5754-O, AA6111-T4 and DP-Steel. The modified Chaboche model was confirmed to well represent the measured hardening behavior including the Bauschinger and transient behavior. While the theoretical and numerical formulations of the constitutive law are discussed in Part I, experimental verifications for spring-back of formed parts are further discussed in Part III.     

Strain gradient crystal plasticity analysis of a single crystal containing a cylindrical void

The effects of void size and hardening in a hexagonal close-packed single crystal containing a cylindrical void loaded by a far-field equibiaxial tensile stress under plane strain conditions are studied. The crystal has three in-plane slip systems oriented at the angle 60° with respect to one another. Finite element simulations are performed using a strain gradient crystal plasticity formulation with an intrinsic length scale parameter in a non-local strain gradient constitutive framework. For a vanishing length scale parameter the non-local formulation reduces to a local crystal plasticity formulation. The stress and deformation fields obtained with a local non-hardening constitutive formulation are compared to those obtained from a local hardening formulation and to those from a non-local formulation. Compared to the case of the non-hardening local constitutive formulation, it is shown that a local theory with hardening has only minor effects on the deformation field around the void, whereas a significant difference is obtained with the non-local constitutive relation. Finally, it is shown that the applied stress state required to activate plastic deformation at the void is up to three times higher for smaller void sizes than for larger void sizes in the non-local material.     

Spring-back evaluation of automotive sheets based on isotropic-kinematic hardening laws and non-quadratic anisotropic yield functions: Part I: theory and formulation

In order to improve the prediction capability of spring-back in automotive sheet forming processes, the modified Chaboche type combined isotropic-kinematic hardening law was formulated based on the modified equivalent plastic work principle to account for the Bauschinger effect and transient behavior. As for the yield stress function, the non-quadratic anisotropic yield potential, Yld2000-2d, was utilized under the plane stress condition. Besides the theoretical aspect of the constitutive law including the general plastic work principle for monotonously proportional loading, the method to determine hardening parameters as well as numerical formulations to update stresses were developed based on the incremental deformation theory and the consistency requirement as summarized in Part I, while the characterization of material properties and verifications with experiments are discussed in Part II and III, respectively.     

A plane stress anisotropic plastic flow theory for orthotropic sheet metals

A phenomenological theory is presented for describing the anisotropic plastic flow of orthotropic polycrystalline aluminum sheet metals under plane stress. The theory uses a stress exponent, a rate-dependent effective flow strength function, and five anisotropic material functions to specify a flow potential, an associated flow rule of plastic strain rates, a flow rule of plastic spin, and an evolution law of isotropic hardening of a sheet metal. Each of the five anisotropic material functions may be represented by a truncated Fourier series based on the orthotropic symmetry of the sheet metal and their Fourier coefficients can be determined using experimental data obtained from uniaxial tension and equal biaxial tension tests. Depending on the number of uniaxial tension tests conducted, three models with various degrees of planar anisotropy are constructed based on the proposed plasticity theory for power-law strain hardening sheet metals. These models are applied successfully to describe the anisotropic plastic flow behavior of 10 commercial aluminum alloy sheet metals reported in the literature.     

Elastic–plastic behavior of steel sheets under in-plane cyclic tension–compression at large strain

Elastic–plastic behavior of two types of steel sheets for press-forming (an aluminum-killed mild steel and a dual-phase high strength steel of 590 MPa ultimate tensile strength) under in-plane cyclic tension–compression at large strain (up to 25% strain for mild steel and 13% for high strength steel) have been investigated. From the experiments, it was found that the cyclic hardening is strongly influenced by cyclic strain range and mean strain. Transient softening and workhardening stagnation due to the Bauschinger effect, as well as the decrease in Young's moduli with increasing prestrain, were also observed during stress reversals. Some important points in constitutive modeling for such large-strain cyclic elasto-plasticity are discussed by comparing the stress–strain responses calculated by typical constitutive models of mixed isotropic–kinematic hardening with the corresponding experimental observations.     

Modeling the effects of material non-linearity using moving window micromechanics

In the analysis of materials with random heterogeneous microstructure the assumption is often made that material behavior can be represented by homogenized or effective properties. While this assumption yields accurate results for the bulk behavior of composite materials, it ignores the effects of the random microstructure. The spatial variations in these microstructures can focus, initiate and propagate localized non-linear behavior, subsequent damage and failure. In previous work a computational method, moving window micromechanics (MW), was used to capture microstructural detail and characterize the variability of the local and global elastic response. Digital images of material microstructure described the microstructure and a local micromechanical analysis was used to generate spatially varying material property fields. The strengths of this approach are that the material property fields can be consistently developed from digital images of real microstructures, they are easy to import into finite element models (FE) using regular grids, and their statistical characterizations can provide the basis for simulations further characterizing stochastic response. In this work, the moving window micromechanics technique was used to generate material property fields characterizing the non-linear behavior of random materials under plastic yielding; specifically yield stress and hardening slope, post yield. The complete set of material property fields were input into FE models of uniaxial loading. Global stress strain curves from the FE–MW model were compared to a more traditional micromechanics model, the generalized method of cells. Local plastic strain and local stress fields were produced which correlate well to the microstructure. The FE–MW method qualitatively captures the inelastic behavior, based on a non-linear flow rule, of the sample continuous fiber composites in transverse uniaxial loading.     

A Finite Element Analysis of a Tapered Flat Sheet Tensile Specimen

A uniaxial tension sheet metal coupon with a tapered instead of a straight gage section has been used for centering the location of diffuse neck and for measuring sheet stretchability in a non-uniform strain field. A finite element analysis of such a tensile coupon made of automotive steel sheet metals has been carried out to assess the effect of the tapered gage section geometry and material plastic strain hardening characteristics on the development of local plastic deformation pattern and local stress state, especially beyond the onset of diffuse necking but before localized necking. In particular, the finite element analysis was used in this study to evaluate the accuracy and reliability of an experimental data analysis method for estimating the post-necking effective plastic stress-strain curve based on the direct local surface axial plastic strain measurements for base metal, heat-affected zone, and weld metals of a dual-phase steel DP600. It is concluded that the estimated lower and upper bounds of the effective stress-strain curve at large strains are not satisfactory for low strain-hardening materials such as heat-affected zone and weld metals with the tapered tension coupons. A simple correction method utilizing only the additional local surface strain measurement in the transverse direction is proposed and it is shown to be effective in correcting the estimated effective stress-strain curve of dual-phase steel weld metals obtained for two tapered gage section geometries.     

考虑材料循环塑性的疲劳裂纹扩展模拟

提出了一种考虑材料循环塑性性能的研究疲劳裂纹扩展与闭合行为的有限元模拟方法．对所选用的循环塑性本构关系进行了基本实验检验．探讨了在疲劳裂纹扩展有限元分析中网格尺寸的影响，给出了网格优化准则．研究了在循环硬化条件下考虑裂纹合效应时裂纹面张开廓形、裂纹尖端应力、应变场和正反向塑性区的演变规律．对于循环硬化和不同循环应力比Ｒ等因素对裂纹张开应力水平的影响也作了考察