Modified Burzynski criterion with non-associated flow rule for anisotropic asymmetric metals in plane stress problems

The Burzynski criterion is developed for anisotropic asymmetric metals with the non-associated flow rule (NAFR) for plane stress problems. The presented pressure depending on the yield criterion can be calibrated with ten experimental data, i.e., the tensile yield stresses at 0°, 45°, and 90°, the compressive yield stresses at 0°, 15°, 30°, 45°, 75°, and 90° from the rolling direction, and the biaxial tensile yield stress. The corresponding pressure independent plastic potential function can be calibrated with six experimental data, i.e., the tensile R-values at 0°, 15°, 45°, 75°, and 90° from the rolling direction and the tensile biaxial R-value. The downhill simplex method is used to solve these ten and six high nonlinear equations for the yield and plastic potential functions, re- spectively. The results show that the presented new criterion is appropriate for anisotropic asymmetric metals.     

A novel quadratic yield model to describe the feature of multi-yield-surface of rolled sheet metals

Since Hill’s quadratic yield model [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. Roy. Soc. Lond. A193, 281–297] cannot address enough experimental results for fairly describing the “anomalous” yield behavior as observed in some of rolled sheet metals, a new quadratic yield model is proposed. As the concept of multiple yielding systems is introduced into the new quadratic yield model, seven commonly used experimental results, three uniaxial tension stresses, one equibiaxial tension stress and three strain ratios, can all be taken into account for characterizing the anisotropy of rolled sheet metals. If more experimental results are extra needed for further improving the prediction, this yield model is still workable. As the experimental parameters are defined as functions of loading direction of corresponding test separately from the major part of yield model, the increase of experimental results regarding the same test does not vary the quadratic form of yield model. The representation of this yield model with axes of principal stresses demonstrates the similar form to Hill’s quadratic model. Therefore, many previous studies developed from Hill’s quadratic yield model can be directly upgraded by the new model to reach a higher accurate level.     

A non-associated flow rule for sheet metal forming

An improved model of material behavior is proposed that shows good agreement with experimental data for both yield and plastic strain ratios in uniaxial, equi-biaxial, and plane-strain tension under proportional loading for steel, aluminum and possibly other alloys. This model is based on a non-associated flow rule in which the plastic potential and yield surface functions are defined by quadratic functions of the stress tensor. The plastic potential aspect of the model is identical to that proposed by Hill for a quadratic anisotropic plastic potential defined in terms of measured r values. The new model differs in that the yield surface, although also defined by a quadratic function of the stress tensor, is defined independently of the plastic potential in terms of measured yield stresses. The model is developed and implemented in an FEM code that is based on a convected coordinate system. Since the associated flow rule, which assumes equivalency between the plastic potential and yield functions, is commonly accepted as a valid law in the theory of plastic deformation of most metals, the arguments for the associated flow rule are also discussed.     

Texture development and plastic anisotropy of B.C.C. strain hardening sheet metals

A Taylor-like polycrystal model is adopted here to investigate the plastic behavior of body centered cubic (b.c.c.) sheet metals under plane-strain compression and the subsequent in-plane biaxial stretching conditions. The <111> pencil glide system is chosen for the slip mechanism for b.c.c. sheet metals. The {110} <111> and {112} <111> slip systems are also considered. Plane-strain compression is used to simulate the cold rolling processes of a low-carbon steel sheet. Based on the polycrystal model, pole figures for the sheet metal after plane-strain compression are obtained and compared with the corresponding experimental results. Also, the simulated plane-strain stress—strain relations are compared with the corresponding experimental results. For the sheet metal subjected to the subsequent in-plane biaxial stretching and shear, plastic potential surfaces are determined at a given small amount of plastic work. With the assumption of the equivalence of the plastic potential and the yield function with normality flow, the yield surfaces based on the simulations for the sheet metal are compared with those based on several phenomenological planar anisotropic yield criteria. The effects of the slip system and the magnitude of plastic work on the shape and size of the yield surfaces are shown. The plastic anisotropy of the sheet metal is investigated in terms of the uniaxial yield stresses in different planar orientations and the corresponding values of the anisotropy parameter R, defined as the ratio of the width plastic strain rate to the through-thickness plastic strain rate under in-plane uniaxial tensile loading. The uniaxial yield stresses and the values of R at different planar orientations from the polycrystal model can be fitted well by a yield function recently proposed by Barlat et al. (1997b).     

Advances in experiments on metal sheets and tubes in support of constitutive modeling and forming simulations

This work is a review of experimental methods for observing and modeling the anisotropic plastic behavior of metal sheets and tubes under a variety of loading paths, such as biaxial compression tests; biaxial tension tests on metal sheets and tubes using closed-loop electrohydraulic testing machines; the abrupt strain path change method for detecting a yield vertex and subsequent yield loci without unloading; in-plane stress reversal tests on metal sheets; and multistage tension tests. Observed material responses are compared with the predictions of phenomenological plasticity models. Special attention is paid to the plastic deformation behavior of materials commonly used in industry, and to verifying the validity of conventional anisotropic yield criteria for those materials and associated flow rules at large plastic strains. The effects of using appropriate anisotropic yield criteria on the accuracy of simulations of forming defects, such as large springback and fracture, are also presented to highlight the importance of accurate material testing and modeling.     

Plane stress yield function for aluminum alloy sheets—part 1: theory

A new plane stress yield function that well describes the anisotropic behavior of sheet metals, in particular, aluminum alloy sheets, was proposed. The anisotropy of the function was introduced in the formulation using two linear transformations on the Cauchy stress tensor. It was shown that the accuracy of this new function was similar to that of other recently proposed non-quadratic yield functions. Moreover, it was proved that the function is convex in stress space. A new experiment was proposed to obtain one of the anisotropy coefficients. This new formulation is expected to be particularly suitable for finite element (FE) modeling simulations of sheet forming processes for aluminum alloy sheets.     

Experimental and numerical study on formability of friction stir welded TWB sheets based on hemispherical dome stretch tests

In order to investigate formability performance and also to obtain guidelines for the stamping process design of friction stir welded TWB (tailor welded blank) sheets, the hemispherical dome stretching test was experimentally performed and the results of the base and friction stir welded samples were compared. Also, in order to better understand the experimental results, numerical analysis was performed. In this work, five automotive sheets, 6111-T4, 5083-H18, 5083-O aluminum alloy, dual-phase steel (DP590) and AZ31 magnesium alloy sheets were considered by (friction stir) welding the same materials. To represent mechanical properties for the numerical analysis, the non-quadratic orthotropic yield function, Yld2000-2d, was utilized for the aluminum alloy and DP590 sheets, while the Cazacu anisotropic/asymmetric yield function was applied for the AZ31 sheet considering different hardening behavior in tension and compression.     

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

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

Endochronic description of plastic anisotropy in sheet metal

It is shown in this paper that an extended form of Hills quadratic yield criterion for anisotropic sheet metal can be derived from an endochronic theory of plasticity. The extended form considers the combined isotropic–kinematic hardening and the anomalous behavior observed in the anisotropic plastic behavior of sheet metals can be accounted for by the concept of kinematic hardening.This form of anisotropic endochronic theory can accommodate the usual requirement of normality between the plastic strain rate and the yield function. In addition, the theory leads naturally to the expressions for back stresses. This work provides an additional example to show that the form of the intrinsic time is directly related to the form of the yield function.It is suggested that the coefficients of the quadratic yield function be determined from the yield stresses obtained from a set of tension tests.     

Forming simulation of aluminum sheets using an anisotropic yield function coupled with crystal plasticity theory

This paper describes the application of a coupled crystal plasticity based microstructural model with an anisotropic yield criterion to compute a 3D yield surface of a textured aluminum sheet (continuous cast AA5754 aluminum sheet). Both the in-plane and out-of-plane deformation characteristics of the sheet material have been generated from the measured initial texture and the uniaxial tensile curve along the rolling direction of the sheet by employing a rate-dependent crystal plasticity model. It is shown that the stress–strain curves and R-value distribution in all orientations of the sheet surface can be modeled accurately by crystal plasticity if a “finite element per grain” unit cell model is used that accounts for non-uniform deformation as well as grain interactions. In particular, the polycrystal calculation using the Bassani and Wu (1991) single crystal hardening law and experimental electron backscatter data as input has been shown to be accurate enough to substitute experimental data by crystal plasticity data for calibration of macroscopic yield functions. The macroscopic anisotropic yield criterion CPB06ex2 (Plunkett et al., 2008) has been calibrated using the results of the polycrystal calculations and the experimental data from mechanical tests. The coupled model is validated by comparing its predictions with the anisotropy in the experimental yield stress ratio and strain ratios at 15% tensile deformation. The biaxial section of the 3D yield surface calculated directly by crystal plasticity model and that predicted by the phenomenological model calibrated with experimental and crystal plasticity data are also compared. The good agreement shows the strength of the approach. Although in this paper, the Plunkett et al. (2008) yield function is used, the proposed methodology is general and can be applied to any yield function. The results presented here represent a robust demonstration of implementing microscale crystal plasticity simulation with measured texture data and hardening laws in macroscale yield criterion simulations in an accurate manner.     

Numerical analysis of diffuse and localized necking in orthotropic sheet metals

In the present paper the diffuse and localized necking models according to Swift [Swift, H.W., 1952. Plastic instability under plane stress, Journal of the Mechanics and Physics of Solids, 11–18], Hill [Hill, R., 1952. On discontinuous plastic states, with special reference to localized necking in thin sheets. Journal of the Mechanics and Physics of Solids 1, 19–30] and Marciniak and Kuczyński [Marciniak, Z., Kuczyński, K., 1967. Limit strains in the process of stretch-forming sheet metal. International Journal of Mechanical Sciences 9, 609–620], respectively, are considered. A theoretical framework for the mentioned models is proposed that covers rigid–plastic as well as elastic–plastic constitutive modelling using various advanced phenomenological yield functions that are able to account very accurately for plastic anisotropy. The mentioned necking models are applied to different orthotropic sheet metals in order to assess their predictive capabilities and to stress out some potential sources for discrepancies between simulations and experiments. In particular, the impact of the applied hardening curve and the equibiaxial r-value, which was recently introduced by Barlat [Barlat, F., Brem, J.C., Yoon, J.W., Chung, K., Dick, R.E., Choi, S.-H., Pourboghrat, F., Chu, E., Lege, D.J., 2003. Plane stress yield function for aluminium alloy sheets – part 1: theory. International Journal of Plasticity 19, 297–1319], on formability prediction is investigated. Furthermore, the impact of the Portevin–LeChatelier effect on the formability of aluminum sheet metals is discussed.     

Plastic strip model of cracked weld joint with residual stresses

The effect of residual stresses on the fracture behavior of a cracked weld joint is studied by making use of the continuous dislocation formulation. Considered are the plastic zone length of the strip model zone and the opening displacement of a crack that is normal to both weld line and base metal boundary; they depend on the character of the yield stresses for the base metal (BM), weld material (WM), and heat affected zone (HAZ). The crack driving force is found to increase with the tensile residual stress while crack initiation and growth are suppressed if the residual stress is compressive. Moreover, the plastic zone and crack opening displacement are found to decrease linearly with the HAZ yield strength as the HAZ width is increased for HAZ yield strength greater than that of BM.     

Plasticity of formable all-metal sandwich sheets: Virtual experiments and constitutive modeling

The mechanical behavior of a metallic sandwich sheet material composed of two flat face sheets and two bi-directionally corrugated core layers is analyzed in detail. The manufacturing of the sandwich material is simulated to obtain a detailed unit cell model which accounts for the non-uniform thickness distribution and residual stresses associated with the stamping of the core layers. Virtual experiments are performed by subjecting the unit cell model to various combinations of bi-axial in-plane loading including the special cases of uniaxial tension, uniaxial compression, equi-biaxial tension and shear. The results demonstrate that the core structure’s contribution to the in-plane load carrying capacity of the sandwich sheet material is similar to that of the face sheets. The numerical results are also used to identify the effective yield surface and hardening response of both the core layer and the face sheets. An anisotropic yield function with linear pressure dependency is proposed to approximate the equal-plastic work surfaces for the core structure and face sheets. Furthermore, a new two-surface model with non-linear interpolation based on plastic work density is presented to describe the observed combined isotropic-distortional hardening of the core structure.     

Constitutive modeling for anisotropic/asymmetric hardening behavior of magnesium alloy sheets

Magnesium alloy sheets have been extending their field of applications to automotive and electronic industries taking advantage of their excellent light weight property. In addition to well-known lower formability, magnesium alloys have unique mechanical properties which have not been thoroughly studied: high in-plane anisotropy/asymmetry of yield stress and hardening response. The reason of the unusual mechanical behavior of magnesium alloys has been understood by the limited symmetry crystal structure of HCP metals and thus by deformation twinning. In this paper, the phenomenological continuum plasticity models considering the unusual plastic behavior of magnesium alloy sheet were developed for a finite element analysis. A hardening law based on two-surface model was further extended to consider the general stress–strain response of metal sheets including Bauschinger effect, transient behavior and the unusual asymmetry. Three deformation modes observed during the continuous in-plane tension/compression tests were mathematically formulated with simplified relations between the state of deformation and their histories. In terms of the anisotropy and asymmetry of the initial yield stress, the Drucker–Prager’s pressure dependent yield surface was modified to include the anisotropy of magnesium alloy. The numerical formulations and characterization procedures were also presented and finally the correlation of simulation with measurements was performed to validate the proposed theory.     

QUASI－FLOW CORNER THEORY ON LARGE PLASTIC DEFORMATION OF DUCTILE METALS AND ITS APPLICATIONS

ＱＵＡＳＩ－ＦＬＯＷＣＯＲＮＥＲＴＨＥＯＲＹＯＮＬＡＲＧＥＰＬＡＳＴＩＣＤＥＦＯＲＭＡＴＩＯＮＯＦＤＵＣＴＩＬＥＭＥＴＡＬＳＡＮＤＩＴＳＡＰＰＬＩＣＡＴＩＯＮＳＨｕＰｉｎｇ（胡平）ＬｉｕＹｕｑｉ（柳玉启）ＧｕｏＷｅｉ（郭威）ＴａｉＦｅｎｇ（台风）（Ｒ...     

Analytical derivation of earing in circular cup drawing based on simple tension properties

In the circular cylindrical cup drawing process of sheet materials, an earing profile develops, incurred by the planar anisotropic properties of sheets. Therefore, proper analysis of earing in cup drawing is important to evaluate anisotropic properties and also to control the development of earing. Even though anisotropic properties are commonly measured in the simple tension test, deformation in circular cylindrical cup drawing is in a near plane strain mode (at the flange) so that numerical simulations utilizing yield functions are common practices to analyze earing. In this work, simplified analytical derivation of earing development in circular cylindrical cup drawing is proposed, based on two simple tension anisotropic properties: the yield stress and the r-value. Good performance of the analytical derivation was verified for AA2090-T3, which has strong anisotropy and six ears in cup drawing. Since the current approach directly utilizes measured simple tension data without involving any yield functions, computational cost is significantly lower. Besides, the current derivation can handle any set of detailed anisotropic measurements in the simple tension test, unlike numerical approaches involving yield functions, which need the development of sophisticated yield functions in the first place.     

Stretch-induced stress patterns and wrinkles in hyperelastic thin sheets

Wrinkles are commonly observed in stretched thin sheets and membranes. This paper presents a numerical study on stretch-induced wrinkling of hyperelastic thin sheets based on nonlinear finite element analyses. The model problem is set up for uniaxial stretching of a rectangular sheet with two clamped ends and two free edges. A two-dimensional stress analysis is performed first under the plane-stress condition to determine stretch-induced stress distribution patterns in the elastic sheets, assuming no wrinkles. As a prerequisite for wrinkling, development of compressive stresses in the transverse direction is found to depend on both the length-to-width aspect ratio of the sheet and the applied tensile strain in the longitudinal direction. A phase diagram is constructed with four different distribution patterns of the stretch-induced compressive stresses, spanning a wide range of aspect ratio and tensile strain. Next, an eigenvalue analysis is performed to find the potential buckling modes of the elastic sheet under the prescribed boundary conditions. Finally, a nonlinear post-buckling analysis is performed to show evolution of stretch-induced wrinkles. In addition to the aspect ratio and tensile strain, it is found that the critical condition for wrinkling and the post-buckling behavior both depend sensitively on the sheet thickness. In general, wrinkles form only when both the magnitude and the distribution area of the compressive stresses are sufficiently large. The wrinkle wavelength decreases with increasing strain, in good agreement with the prediction by a scaling analysis. However, as the tensile strain increases, the wrinkle amplitude first increases and then decreases, eventually flattened beyond a moderately large critical strain, in contrast to the scaling analysis.