Algorithmic Formulations of Evolutionary Anisotropic Plasticity Models Based on Non-Associated Flow Rule

In the present paper, orthotropic elasto-plastic constitutive formulations for sheet metal forming based on non-associated flow rule that assume distortion of yield function/plastic potential with ongoing deformation process are analyzed. The yield function/plastic potential are considered as two different functions with functional form as orthotropic quadratic Hill or non-quadratic Karafillis-Boyce stress function. Based on the principle of plastic work equivalence, anisotropy parameters of the utilized yield function/plastic potential are set as functions of the equivalent plastic strain. In the constitutive formulation, for this internal variable, evolution equation consistent with the same principle of plastic work equivalence is introduced. For DC06 sheet sample with reported significant variation of the incremental r-values with straining, predictions of the evolution of the yield stress and r-value directional dependences with straining obtained by the analyzed models are presented. The algorithmic formulations of the analyzed constitutive models are derived by application of the implicit return mapping algorithm. For the derived stress integration procedures the accuracy is investigated by calculating iso-error maps. The maps are compared according to the flow rule and involved orthotropic stress functions. It has been revealed that although there is a difference in maps configuration there is no prominent difference in error magnitudes.


INTRODUCTION
Sheet metals are produced by complex thermo-mechanical processing routes among which the production rolling steps introduce a preferential orientation to the grains called texture.This preferential orientation results in anisotropic material behavior that is closely related to the formability of the material.The finite element codes that are widely used to design sheet forming processes most Latin American Journal of Solids and Structures 14 (2017) 1853-1871 frequently utilize the phenomenological plasticity constitutive theories that are based on yield surface concept and plastic potential theory.The basis of these phenomenological constitutive models are: a yield criterion representing surface that separates the elastic and plastic regions in the stress space, a plastic potential whose gradient represents direction of the plastic strain rate and a hardening rule by which evolution of the yield surface is described.
The simplest and the most utilized approach in phenomenological plasticity theories that are intended for the metallic materials is the concept of isotropic hardening and use of the assumption known as associated flow rule or normality condition.
Associated flow rule assumption implies that the yield function and plastic potential function are identical and it is often regarded as the necessary condition for the metal plasticity fundamentals (stability of plastic flow and uniqueness).However, according to some theoretical observations (Mroz, 1963;Runesson and Mroz, 1989;Stoughton, 2002;Stoughton andYoon, 2004, 2006;Cvitanić et al., 2008), an appropriate plastic flow description can be achieved by using less restrictive constrains over the constitutive equations than the associated flow rule.By non-associated flow rule, yielding and plastic flow are described by two different functions.In the last few years, several metal plasticity non-associated formulations were developed and this approach became particularly engaging in sheet metal plasticity.An acceptable constitutive model for a sheet material based on associated flow rule requires a yield function that can simultaneously predict anisotropy of yield stresses and anisotropy of plastic flow.Under the assumption of associated flow rule, various anisotropic phenomenological yield functions have been proposed for metallic sheets (Hill, 1948(Hill, , 1979(Hill, , 1990(Hill, , 1993;;Barlat and Lian, 1989;Karafillis and Boyce, 1993;Barlat et al., 1991Barlat et al., , 1997Barlat et al., , 2003Barlat et al., , 2005;;Banabic et al., 2003Banabic et al., , 2005;;Aretz, 2005;and others).Anisotropy is introduced into these stress functions by the parameters determined by the data related to the directional dependences of the yield stresses and Lankford parameters as well as data related to the biaxial tensile tests.Lankford parameter, also known as r-value or plastic strain ratio, is used as the measure of plastic flow.It is defined as the ratio of the sheet specimen transverse and thickness true plastic strain increments in uniaxial tensile testing.According to the standards, r-value is calculated by linear regression of the transverse versus longitudinal plastic strain plot between certain limits of the measured strains.Hence, it is considered as a constant value.In orthotropic sheet metal models based on nonassociated flow rule, yield function and plastic potential function might have an identical functional form but their parameters can be defined by different experimental data.That means that the yield function can be adjusted to the yield stresses and the plastic potential can be adjusted to r-values.Such formulation can address sufficient number of experimental data and can result in acceptable predictions of the uniaxial material behavior even if yield/plastic potential stress functions with a low number of anisotropic parameters are utilized (Lademo et al., 1999;Stougthon, 2002;Stoughton and Yoon, 2004;Cvitanić et al., 2008;Safaei et al., 2014).
Isotropic hardening concept implies a proportional expansion of the yield surface with ongoing deformation process, without any change of its shape and position.In constitutive formulations intended for sheet metals, by application of isotropic hardening and standard calculation procedure for r-value, anisotropy parameters are calculated using the initial yield stress ratios and/or constant r-values.Such approach results in fixed anisotropy parameters and accordingly any possible distortion of yield function/plastic potential with continuation of the plastic deformation process is ne-Latin American Journal of Solids and Structures 14 (2017) 1853-1871 glected.In addition, uniform size change of the orthotropic yield function in the stress space is controlled by single scalar hardening variable and mostly flow stress curve for the rolling direction is adopted as the expansion rule.However, some recent experimental studies report alternation of the yield stress ratios and/or instantaneous r-values with evolution of sheet texture during deformation process (Zamiri and Pourboghrat, 2007;An et al., 2013;Safaei et al., 2014).Furthermore, numerous studies related to the application of the orthotropic plasticity formulations with constant anisotropy parameters in predicting complex forming processes indicate that possible model improvements could be achieved by incorporating the evolution of yield stress ratios and r-values into the model.This is particularly evident in the simulations of the cylindrical cup drawing problem.In this forming process, the plastic anisotropy has important effect on the formation of undulative top edge of the finally drawn cup often termed as earing.Numerous phenomenological orthotropic yield functions with constant (fixed) anisotropy parameters under associated or non-associated flow rule have been tested in this forming process (Yoon et al., 1998(Yoon et al., , 1999(Yoon et al., , 2000(Yoon et al., , 2004(Yoon et al., , 2006;;Cvitanić et al., 2008;Taherizadeh et al., 2010Taherizadeh et al., , 2011;;Park and Chung, 2012;Safaei et al., 2013;Vrh et al., 2014;and others).These extensive numerical studies clearly indicate that the predicted earing trend (location of peaks and valleys) is the mirror image of the predicted r-value trend with the respect to the transverse direction, and that the prediction of the earing profile amplitude (the maximum difference in cup heights) is correlated to the prediction of the yield stress anisotropy amplitude.Considering predictions of the complex phenomenological yield functions that are capable to accurately predict directional dependence of the r-values, it can be observed that there is still discrepancy between their predictions of earing profile and experimental cup heights.Although the description of the material behavior is not the only factor that influences the finite element predictions, due to the expressed correlation between predicted r-values and predicted earing profile, it can be concluded that reliable material model should take into account not only directional dependences of the uniaxial material data determined at the initial yield state or their averaged values but also their possible evolution with continuation of the deformation process.That means that the anisotropy parameters of the yield function/plastic potential should be altered by a certain measure of plastic flow.
In this paper, constitutive formulations based on non-associated flow rule and distortional evolution of the orthotropic yield function/plastic potential with ongoing deformation process are considered.As yield function/plastic potential orthotropic four parametric quadratic stress function proposed by Hill (1948) or non-quadratic stress function proposed by Karafillis-Boyce (1993) are utilized.Anisotropy parameters of the utilized yield function/plastic potential are set as polynomial functions of the equivalent plastic strain and evolution equation for this hardening parameter consistent with the principle of the plastic work equivalence is introduced.Previously developed algorithmic formulations of the evolutionary anisotropic model based on the considered stress functions and associated flow rule (Cvitanić et al., 2016) are now extended for the case of non-associated flow rule.Considered algorithmic formulations are derived by application of the implicit return mapping algorithm.Numerical analysis of accuracy of the proposed stress integration procedures is performed and compared according to the flow rule and the functional form of the involved stress functions by calculating iso-error maps.The constitutive descriptions are developed and analyzed considering experimental data for DC06 steel sheet reported by Safaei et al. (2014).Therefore, for the hardening material the plastic multiplier obeys the complementary conditions 0, 0, Latin American Journal of Solids and Structures 14 (2017) 1853-1871 3 ANISOTROPIC STRESS FUNCTIONS WITH ANISOTROPY EVOLUTION The analyzed constitutive formulations utilize yield function and plastic potential with functional form of orthotropic quadratic Hill (1948) stress function or non-quadratic Karafillis-Boyce (1993) stress function.In sheet metal forming it is common practice to assume that the sheet is approximately subjected to plane stress conditions and that material exhibits orthotropic symmetry in plastic properties.Accordingly, the utilized yield functions/plastic potentials are stated as functions of in-plane stress components xx s , yy s and zz s , where x-axis denotes the original sheet rolling direction and y-axis denotes the direction in sheet plane transverse to the rolling direction.The z-axis denotes the sheet normal direction.
The utilized Hill stress function and Karafillis-Boyce stress function inclose four parameters by which material plastic anisotropy can be described.Under the associated flow rule, the utilized stress function acts as yield function as well as plastic potential, therefore, its anisotropy parameters can be adjusted to the yield stresses or experimental data indicating plastic flow.In the nonassociated formulation, anisotropy parameters of the Hill or Karafillis-Boyce yield function are defined in terms of three directional yield stresses obtained in the uniaxial tension of the specimens oriented at 0, 45 and 90 to the rolling direction and equibiaxial yield stress.The associated yield stresses are denoted as 0 s , 45 s , 90 d e are width and thickness plastic strain increments, respectively, obtained in uniaxial sheet specimen tension.In calculating parameters of the plastic potential, three plastic strain ratios obtained in the uniaxial tensions along 0, 45 and 90 to the rolling direction and the yield stress for the referent direction 0 s are used.The associated plastic strain ratios are denoted as 0 r , 45 r and 90 r .

Hill Stress Function and Karafillis-Boyce Stress Function
The orthotropic quadratic Hill (1948) If the Hill stress function is adjusted to the plastic strain ratios, the calculation procedure results in following explicit expressions The orthotropic Karafillis-Boyce (1993) stress function is a linear combination of two convex non-quadratic functions where c is a weighting parameter and exponent m is an even number.In the above equation, 1 s  , 2 s  and 3 s  are the principal values of the so-called isotropic plasticity equivalent stress tensor.For plane stress conditions, these values can be calculated as where In Karafillis-Boyce stress function material anisotropy is defined by the parameters C , 1 a , 2 a and 3 g .For isotropic material, these parameters have the values and stress components defined by Eq. ( 12) reduce to the components of the stress deviator tensor.
Latin American Journal of Solids and Structures 14 (2017) 1853-1871 For Karafillis-Boyce stress function, calibration procedure i.e. calculation of anisotropy parameters, leads to the system of non-linear equations that can be solved using a numerical iterative procedure.

Evolution of Hill and Karafillis-Boyce Yield Function/Plastic Potential
In this paper, experimentally determined directional dependences of the uniaxial material properties and their evolution with ongoing plastic deformation process for 0.8 mm thick DC06 steel sheet reported by Safaei et al. (2014) are used to develop and analyze evolutionary anisotropic constitutive models.The utilized data are related to the uniaxial straining tests of the seven specimens with orientations 0°, 15°, 30°, 45°, 60°, 75° and 90° to the rolling direction.The reported data utilized in model development are: 1) parameters of the combined Swift-Voce hardening law by which the experimental true stress and longitudinal true plastic strain are approximated; 2) parameters of the 3 rd order polynomial fit by which the experimental transverse and longitudinal true plastic strains are approximated.In the previous study conducted by the authors (Cvitanić et al., 2016), using these data, parameters of the Hill and Karafillis-Boyce yield function/plastic potential are introduced as polynomial functions of the equivalent plastic strain.In the first step, following procedure based on the principle of plastic work equivalence presented in Safaei et al. (2014), for each specimen orientation yield stresses y s and instantaneous r-values corresponding to the certain amount of the equivalent true plastic strain are calculated.In the adopted approach, longitudinal true plastic strain in the rolling direction is used as the equivalent true plastic strain.In calculating instantaneous r-values incompressibility hypothesis is applied, thus Eq. ( 6) is considered in the following form e + are the increments of true longitudinal and transverse plastic strains corresponding to the loading direction q and direction 90 q +  , respectively.According to Eq. ( 14), instantaneous r-values corresponding to the certain longitudinal true plastic strain can be calculated using the slope m q of the available 3 rd order polynomial fit by which the experimental transverse and longitudinal true plastic strains are approximated.
According to the adopted calculation procedure, yield stresses and r-values corresponding to the seven orientations (0°, 15°, 30°, 45°, 60°, 75° and 90°) and amounts of the equivalent plastic strain starting from 0.001 to 0.301 at each 0.002 increment are calculated.The calculated yield stress ratios (yield stresses normalized with yield stress for the rolling direction) and r-values corresponding to the selected values of the equivalent plastic strain are presented in Figure 1 and Figure 2, respectively.From Figure 1 it can be observed that the directional dependence trend of the yield stress ratios at the start of plastic deformation process ( ) is rather distorted with further straining.Furthermore, as presented in Figure 2, for the considered material, there is a significant decrease of the r-values and certain alternation of the r-value directional dependence with ongoing deformation.These observations clearly indicate the suitability of the non-associated flow rule approach, by which the evolution of yield stress and r-value directional dependences can be independently considered and described.In Figures 1 and 2, predictions of the yield stress and r-value directional dependences corresponding to the selected equivalent plastic strains obtained by the Hill yield function/plastic potential are also presented.The predictions of the Karafillis-Boyce yield function/plastic potential are almost identical to the predictions obtained by the Hill functions and therefore they are not separately presented.

Latin
The anisotropy parameters of Hill and Karafillis-Boyce yield functions are calculated using the yield stresses corresponding to the orientations 0, 45 and 90 ( 0 45 90 , , s s s ) and assuming that the and 90 ( 0 45 90 , , r r r ).The parameters of the analyzed Hill and Karafillis-Boyce yield function/plastic potential are calculated considering calculated yield stresses and r-values corresponding to the equivalent plastic strain values starting from 0.001 to 0.301 at each 0.002 increment.From Figures 1 and 2, it can be observed that the analyzed yield functions poorly predict the pronounced directional dependence at the onset of the deformation process while the analyzed plastic potentials result in good predictions of the r-value directional dependences, particularly for greater strain levels.
In Figure 3 predictions of yield stress ratios obtained by the models based on the associated flow rule and Hill or Karafillis-Boyce stress function are presented.In these models, utilized orthotropic stress function adjusted to r-values acts as plastic potential and as yield function.It can be observed that the analyzed associated models fail in predicting directional dependence of the yield stress ratios for the considered deformation levels.These results in addition support the use of the non-associated flow rule approach coupled with analyzed four parametric stress functions in describing plastic anisotropy and its evolution for the considered material.
In order to relate the anisotropy parameters with the equivalent plastic strain, the fourth order polynomial fit is used The calculated polynomial parameters for the analyzed yield functions/plastic potentials are obtained using the least square method and are presented in Table 1.Figures 4 and 5 show values of each anisotropy parameter corresponding to the selected values of the equivalent plastic strain and related polynomial fit.From these figures it can be seen that there is a good correlation between adopted fit and plastic potential parameters, while there is a certain discrepancy for yield function parameters for the lower plastic levels.Contours of the analyzed yield functions and plastic potentials in normalized stress space for zero shear stress corresponding to the selected values of the equivalent plastic strain are presented in Figure 6 and Figure 7, respectively.For material that obeys the isotropic hardening concept, yield contours that correspond to the different amount of plastic strain should coincide if presented in normalized stress space.Considering obtained contours, it can be observed that there is a significant discrepancy between initial and subsequent yield contours and notable evolution of plastic potential contour particularly for the quadratic potential function.

STRESS INTEGRATION PROCEDURE
In this section, for the analyzed evolutionary anisotropic plasticity models based on non-associated flow rule, computational procedure for calculating state variables at time t n+1 ( ) and known increment of total deformation Dε is derived.The derived procedure is based on the implicit return mapping (Simo and Hughes, 1988;Yoon et al., 1999;Cvitanić et al., 2008) and presents an extension of the procedure previously derived for evolutionary anisotropic model based on the associated flow rule (Cvitanić et al., 2016).By the application of the implicit return mapping procedure, the stress solution is obtained in two steps.In the elastic predictor step, the strain increment Dε is assumed to be elastic and trial elastic stress tensor trial σ is calculated based on the previously converged solution σ n : If the trial state violates the yield condition, the plastic correction step is performed assuming trial state as initial condition.In this step the final stress and evolution equations for the plastic strain tensor and hardening parameter are integrated to restore the consistency condition.By application of implicit Euler backward integration procedure and assuming non-associated flow rule, increment of plastic strain tensor is approximated as follows 1 ( ( , )/ ) Consistently, the increment of the hardening parameter defined by Eq. ( 5) is approximated as follows 1 ( ( , ) / ( , )) In above approximations, l D is an incremental consistency parameter that obeys discrete form of the complementary conditions By using Eqs.( 16)-( 19), incremental form of the constitutive model can be stated by the following system of four non-linear algebraic equations (tensorial and scalar) ) can be obtained using iterative Newton-Raphson procedure.At each iteration, denoted as k, above four equations are linearized around the current values of state variables to obtain : ( / ) 0 : ( ) : : ( ) : where Finally, updated state variables for non-associated formulation are defined as  If the associated flow rule is assumed, according to Eq. ( 19) p e l D = D holds and integration scheme proposed by Cvitanić et al. (2016) is reproduced.On the other hand, if anisotropy evolution is neglected orthotropic parameters of yield function/plastic potential are not altered by hardening parameter and integration scheme conform the one developed by Cvitanić et al. (2008).

NUMERICAL ANALYSIS OF ACCURACY. ISO-ERROR MAPS
In order to estimate the accuracy of the proposed algorithm based on non-associated flow rule and distortional evolution of the yield and the plastic potential function, iso-error maps (Simo and Hughes, 1988) are calculated.The iso-error maps are generally accepted as an effective and reliable tool for assessing the accuracy of the constitutive integration algorithms under realistic strain steps.For two-dimensional implementation of the developed algorithms, iso-error maps are generated following known procedure.For the selected stress point on the yield surface sequence of two-component strain increments is applied and for each strain increment numerical solution is obtained according to the tested algorithm.The solution obtained by the same algorithm by dividing the considered strain increment into large number of sub-increments of the same size is adopted as exact solution.Furthermore, for each considered strain increment, based on the obtained numerical and exact solution, the error is calculated as the percentage of the relative root mean square of errors between the computed stress tensor σ and exact stress tensor * σ .This error measure is defined as Finally, for the considered stress point iso-error map is drawn as the contour plot of the error field presented over the two-component strain space.By inspection of these plots obtained for the representative stress points, dependence of the algorithm accuracy on the strain increment size can be assessed.
For the analyzed formulations, iso-error maps are calculated at three representative stress points on the yield surface: A-uniaxial, B-balanced biaxial and C-pure shear as shown in Figure 8. applies.As the exact solution, the stresses obtained by 100 substeps of each strain increment are utilized for the Hill formulations.For the Karafillis-Boyce formulations solutions obtained by 50 sub-steps are adopted.
In the tested algorithms, the polynomial relations for the anisotropic parameters of the Hill stress functions and Karafillis-Boyce stress functions corresponding to DC06 steel sheet that are presented in Section 3 are utilized.Furthermore, combined Swift-Voce stress-strain relation for the rolling direction reported by Safaei et al. (2014)  Figures 9 and 10 show calculated iso-error maps for three representative stress points obtained by the Hill and Karafillis-Boyce formulations, respectively.Calculated iso-error maps are compared to the maps obtained by the algorithm based on associated flow rule and distortional hardening concept coupled with the yield function adjusted to r-values (Cvitanić et al., 2016).For the analyzed formulations, the magnitudes of errors are comparable to the results reported in literature obtained by the associated/non-associated formulations based on the utilized orthotropic stress functions and isotropic hardening concept (Cvitanić et al. 2008).Commonly, considering points B and C, for the isotropic formulations the exact solutions are obtained for the loading directed along the yield surface symmetry axis.As shown in Figures 9 and 10, considering these points, for the analyzed orthotropic formulations the axes of exact solution are shifted to the stress symmetry axes.From the comparison of the iso-error maps for different stress points, it can be seen that the errors are relatively smaller for the balanced biaxial stress state for all analyzed formulations.Considering different flow rule approach, for the Hill formulations there is no significant difference between isoerror maps obtained by the associated or non-associated flow rule.For the Karafillis-Boyce formulations, the errors are relatively smaller if associated flow rule is utilized.Furthermore, from the comparison of the iso-error maps for different orthotropic stress functions it can be seen that although there is a difference in maps configuration there is no prominent difference in error magnitudes.Finally, it can be stated that for the analyzed formulations, a reduction in error magnitude can be expected when the strain size increment is reduced.

CONCLUSIONS
In the present paper, constitutive formulations for sheet metal forming based on non-associated flow rule that enable distortion of the yield function/plastic potential with on-going deformation process are analyzed.The formulations are developed considering evolution of the yield stress ratios and rvalues with straining obtained in uniaxial tensile tests for DC06 steel sheet samples with different alignment to the rolling direction.As yield function/plastic potential, simple orthotropic quadratic Hill (1948) or non-quadratic Karafillis-Boyce (1993) stress functions are utilized.According to the non-associated flow rule approach, the anisotropy parameters of the yield function are adjusted to the data related to the evolution of the yield stress ratios while the anisotropy parameters of the plastic potential are adjusted to the data related to evolution of the r-values.Considering these yield stress ratio and r-value evolutions, based on the principle of the plastic work equivalence, the anisotropy parameters of the utilized stress functions are set as 4 th order polynomial functions of the equivalent plastic strain.For the considered DC06 sheet sample, the analyzed non-associated formulations result in acceptable predictions of the yield stress and r-value directional dependences and their evolutions with on-going deformation.Besides, considering yield function/plastic potential contours, evolution of the plastic potential is more pronounced especially for the quadratic plastic potential function.Therefore, for the considered material, the presented results clearly indicate advantage of the non-associated over the associated flow rule approach.Furthermore, in the constitutive formulation, for the equivalent plastic strain as the internal variable that governs the evolution of the yield function/plastic potential, evolution equation consistent with the same principle of plastic work equivalence is utilized.
In the present paper, the algorithmic formulations of the analyzed non-associated constitutive descriptions with yield function/plastic potential evolution are developed based on the implicit Euler backward integration procedure.Application of this procedure results in the system of four nonlinear algebraic equations (tensorial and scalar) that present incremental form of the constitutive equations.For the obtained system, the solution procedure based on the iterative Newton-Raphson procedure is developed.Considering data for DC06 sheet sample, accuracy of the derived stress integration procedures has been estimated by calculating iso-error maps for three representative stress points.For the analyzed non-associated formulations with yield function/plastic potential evolution, the error magnitudes are comparable to the error magnitudes obtained for the associated/non-associated formulations based on the same orthotropic stress functions and the isotropic hardening concept.It has been revealed that although there is a difference in maps configuration there is no considerable difference in error magnitudes for the formulations based on the different flow rule approach and the quadratic or the non-quadratic stress function.Furthermore, based on the obtained results it can be stated that for the analyzed formulations, a reduction of error magnitude can be expected when the strain size increment is reduced.
plastic flow are described by the same function and Eq.(3) becomes the so-called associated flow rule.In the present formulation, the parameter p e is considered as equivalent plastic strain that obeys the following form of the principle of plastic work equivalence process is elastic, the incremental changes of the internal variables (plastic strain tensor ε d p and parameter p e ) vanish and 0 dl = .
s and b s .The anisotropy parameters of the Hill or Karafillis-Boyce plastic potential are defined in terms of experimental data indicating plastic flow such as Lankford parameter that reads

Figure 1 :
Figure 1: Yield stress ratio directional dependences corresponding to several values of the equivalent plastic strain.Predictions of yield stress ratios obtained by Hill yield function.

Figure 2 :
Figure 2: Lankford parameter directional dependences corresponding to several values of the equivalent plastic strain.Predictions of Lankford parameters obtained by Hill plastic potential.

Latin
American Journal of Solids and Structures 14 (2017) 1853-1871yield stress at balanced biaxial stress state is the averaged value of the yield stresses corresponding to the rolling and transverse direction of the Hill and Karafillis-Boyce plastic potentials are calculated using r-values corresponding to the orientations 0, 45

Figure 3 :
Figure 3: Yield stress ratio directional dependences predicted by associated model and a) Hill ; b) Karafillis-Boyce function adjusted to r-values corresponding to several values of the equivalent plastic strain.
By solving above linearized system, explicit expressions for the increments of state variables are obtained

Figure 8 :
Figure 8: Plane stress yield surface and points A, B and C for iso-error maps.

Figure 9 :Figure 10 :
Figure 9: Iso-error maps based on Hill associated formulation for points: a) A; b) B; c) C and on Hill non-associated formulation for points: d) A; e) B; f) C.

Latin American Journal of Solids and Structures 14 (2017) 1853-1871
stress function is derived as an extension of the isotropic von Mises yield function.For plane stress condition Hill stress function can be written in the following form n and r are anisotropic material parameters that can be adjusted to experimental data and y s is the yield stress for the referent direction.For Hill stress function calculation proce- dure results in the explicit expressions for anisotropy parameters.The parameters of the Hill stress function adjusted to the yield stresses read Latin American Journal of Solids and Structures 14 (2017) 1853-1871

Table 1 :
Parameters of Poly4 fit for Hill and Karafillis-Boyce yield function/plastic potential.
for the considered material is used to govern the yield surface expansion