Improved Approach to Determine the Material Parameters for a Combined Hardening Model

For the simulation of isothermal mechanically loaded components, it is indispensable to have a material model, which describes the material behavior very accurately. In this case, a combined hardening model was chosen in order to reflect the prevalent deformation behavior. The combined hardening model enables simulation independent of the number of load cycles and the chosen strain amplitude. The main point is the declaration of the parameters from the chosen material model. This work deals with the estimation of the parameters. For validation and as input data of the here defined approach low cycle fatigue (LCF) tests were performed on cast aluminum and at 250 ̊C. The comparison of the test results and the simulations indicated that σmax from the simulated hysteresis lies inside a range of ±5% referred to the test results.


Introduction
The result of the pursuit for more efficient and lighter components is that the loading conditions and the requirements on the materials are increasing.If the loading conditions exceed the elastic material behavior plastic strain occurs [1].This means that material behavior changes from an elastic into a non-proportional inelastic behavior.Further factors, which influence the material behavior, are the temperature and the rate how an external load is applied.If components are subjected to such loading conditions, it is necessary to understand the material's response.With the results from isothermal LCF tests, it is possible to examine the inelastic material behavior.For designing and simulation the material behavior has to be pictured out by mathematically models.One popular model, which covers such effects, is the Ramberg-Osgood-model [2].The model deals with the description of elastic and inelastic areas from stress strain curves.Two linear kinematic models, which take the plastic range of the material behavior into account, are the Prager [3] and Ziegler [4] hardening rules.Based on the linear kinematic model from Prager, Amstrong and Frederic [5] developed a nonlinear kinematic model.The time independent behavior from the yield stress is considered in the isotropic part of the Chaboche-model [6], which also includes the nonlinear kinematic model from Amstrong and Frederic.
In this paper, the material behavior of cast aluminum is declared by using a At the end, finite element simulation with Zebulon ® was carried out.The present work will focus on the estimation of the parameters based on the LCF test results.

Material Model
A combined hardening model has been chosen for the reproduction of the material behavior.Based on the Chaboche model the combined hardening law include a kinematic and isotropic hardening law.The main equation that expresses the used model is: where F < 0 represents the elastic part and if F = 0 plastic flow takes place.The evolution of the yield surface is taken into account in the back stress α and in the drag stress R [7] [8].Thereby the back stress α represents the translation and the drag stress R the cyclic hardening behavior.The yield surface obeys the Von Mises criterion: The back stress R represent the isotropic part, which deals with the increase or decrease of the yield surface, depending on the accumulated plastic strain.R 0 conforms to the initial yield stress by zero accumulated plastic strain p  .The back stress R is described by following equation: For the consideration of the translation of the yield stress the non-linear kintmatic hardening law of Frederick-Armstrong was used.Thereby the Prager hardening rule is replaced through the Ziegler hardening rule.In this case the Materials Sciences and Applications kinematic hardening is given through following equation: ( ) As shown in [9] [10] and under consideration of [11] the explicit integration of the kinematic hardening Equation ( 4) over an uniaxial strain cycle with the first data pair (0, α 1 ) gives following equation:

Experimental Procedure
As input data for the CHP routine LCF test results were used.The tests were carried out on a servo-hydraulic testing machine type 8802 produced by the company Instron ® GmbH.An Instron ® GmbH controller type Fast Track 8800 handled the data collection and the controlling.All tests were made strain rate controlled and under a strain-ratio of R ε = −1.The tested material was cast aluminum and Figure 1 shows the test specimen geometry.
Tests were carried out at 250˚C, different heat treatments and by different strain amplitudes, as Table 1 shows.All LCF investigations were performed in consideration of the recommended standards for LCF testing, according to DIN EN 3988 and ISO 12106.

Parameter Definition
If a material model was selected a major rule plays the determination of the parameters, in regard to be able to generate exact simulation.To deal with the estimation of the material parameters the CHP routine was developed.The routine can be divided into two parts.The first part handles the estimation of the parameters and the simulation from the cycle on which the kinematic parameters were declared.If this simulation match with the test results, the second part of the CHP routine can be follow.The second part looks whether the data set from the first part can also describe the hysteresis at any other cycle.In the case that the simulation don't match, the isotropic parameters were adjusted.Every CHP routine controlled simulation were realized with Zebulon ® .The creation of the input file and the editing of the simulation results were done by the CHP routine.It has to be pointed out that for the use of the CHP routine following boundary conditions has to be complied from the test data: • The shape from the hysteresis has to be symmetric.
• The stress curve in the nonlinear range has to be identical at each cycle.
The most aluminum and steel alloys show such a behavior [12] [13].For materials (e.g.: grey cast iron), which not conform to the first determined material behavior, other material models has to be chosen, as demonstrated in [14] [15].This cycle serves also as initial cycle for the simulation step.The fifth step is the simulation step.Therefore, the parameters from the third and fourth step had been combined to one material model and saved in the input file.After simulation, the output file is read by the CHP routine and compared with the test data from the simulation cycle.
The most critical point in the estimation of the material parameters are the definition of the R values, which is done in the second step of CHP routine.For the definition of R the young modulus of each cycle has to be declared.Zhou et al. [16] and Fournier et al. [8] use a linear regression on the points located between two offset parameters.The two offset parameters were declared at the beginning of the second step.In this case, the offset parameters are determined by σ offsetmin and σ offsetmax as shown in Figure 3(a).After this, the regression line will be displaced with the offset parameter ε offset and cut with the hysteresis arm.This is plotted in Figure 3(a).The point of intersection represents the crossover from elastic to the plastic material behavior.Furthermore also the R value will be specified by this point.

Results and Discussion
To extract the combined hardening parameters, from the executed LCF tests, the CHP routine was used.All analysis were fitted at the lower hysteresis arm and the boundary conditions were chosen as following: • σ offsetmax = 0.5 * σ max • σ offsetmin = 0.05 * σ min • ε offset = 0.00001 In Figure 5 the tested hysteresis, of the 2000 th cycle with a heat treatment of T6, are printed out.The reference value is the mean value σ mean of the σ testmax of the currently investigated hysteresis.Figure 5 shows that σ testmax lies inside a scatter of ±5% regarding to σ mean .On basis of these test data the combined hardening parameter for the 2000 th cycle were generated with the CHP routine.The simulated hysteresis, based on the determined parameter, is shown in Figure 5.
The σ max of the simulated hysteresis lies within the range of ±5% regarding to σ mean .
In order to check out if the parameters from above can describe the material behavior also by another cycle, the second part of the OCP routine is used.In Table 2 the σ max of the selected test cycles are given.
The first simulation displayed that for some cycle the variance between σ max and σ mean is greater than 5%.Therefore, an optimization of the isotropic hardening parameters took place.The results of the simulation, with the optimized parameter, shows that the σ max for each cycle are inside the range of ±5%.The As input data for the isotropic fit all data up to the respective cycle N f were used to calculate the parameters.The kinematic parameters were determined by    the data set with the highest strain amplitude and cycle N f .The described simulated hysteresis in Figure 7 are based on the optimized material parameter set.
Table 3 shows that for each cycle and strain amplitude the σ max from the simulated hysteresis are inside a range of ±5% regarding to σ testmax .

Conclusions
The test results of the T6 heat-treated material showed that the variation of σ max is in a scatter of ±5%.Also the simulation result for the 2000 th cycle is inside this range.However, at the end of the CHP routine also the simulation results of the other chosen cycles are considering the default range of ±5%.
On the basis of the parameters, which were detected from the as cast specimens by strain amplitude of ε ampl = 0.4%, simulation by different strain amplitudes were carried out.The comparison of the different simulation results and the test data demonstrate that σ max of the simulated hysteresis always lies in a scatter band of ±5% regarding to σ mean of the tested hysteresis by the same cycles.
Overall, in relation to the CHP routine following conclusions can be made: • The simulations describe the material behavior in a range which is sufficient.
Thereby the critical point, during the estimation of the material parameters, is the determination of the boundary condition for fitting the young's modulus and the R values.
• The CHP routine determinate parameters that are able to simulate the stabilized cycle N f /2 or any cycle between one and N f .• The transportability of the estimated parameters is given, if the used finite-element-method includes the combined hardening law.
In the further research, the application of the highlighted routine should be widened so that visco-plastic material behavior could be considered.

B.
Seisenbacher et al.DOI: 10.4236/msa.2018.94024358 Materials Sciences and Applications combined hardening model.A modified Chaboche model represents the combined hardening model.The Chaboche model is changed by replacing the Prager hardening rule through the Ziegler hardening rule.The challenge is the estimation of the parameters of the combined hardening model.Therefore, a Matlab ® routine was established, which deals with the determination of the parameters from the combined hardening model.As input data for the combined hardening parameter (CHP) routine the hysteresis of the conducted LCF tests were taken.

Figure 1 .
Figure 1.Drawing of the test specimen, which was used for LCF tests.

Figure 2
Figure2present an overview of the steps which are executed in the first part of the CHP routine.At first, the input data must be loaded and the boundary condition for the second step has to be specified.The analysis of the input data takes place in the second step.Thereby the significant values, such as σ max , E, R, of each hysteresis were declared.Based on these values, in the third step the isotropic parameters were estimated.The fourth step deals with the determination of the kinematic hardening parameters.At the beginning of step four the cycle, at which the kinematic hardening parameters were declared, has to be defined.

Figure 2 .
Figure 2. Flow chat from the first part of the CHP routine.

Figure 3 .
Figure 3. Preparing data for the isotropic parameter fit.(a) Determination of Young's Modul and R values, (b) isotropic hardening.

Figure 4 .
Figure 4. Flow chat from the second part of the CHP routin.

Figure 5 .
Figure 5.The figure shows test results based on the heat treated material of the 2000 th cycle and the simulation result of this cycle.

Figure 6 .
Figure 6.Comparison of the Simulation and test results of heat treatment material and different cycles.

Figure 7 .
Figure 7.Comparison of the Simulation and test results by non heat treated material, different cycles and strain amplitudes.
the σ testmax values σ offsetmax , σ offsetmin parameters to the input data for fitting the young modulus σ testmax maximal stress of a tested hysteresis

Table 2 .
Test results from the heat treated specimens.

Table 3 .
Results of the non-heat treatment material.