1. Introduction
Nowadays, finite element modelling (FEM) has attained a great deal of importance when dealing with the design of new parts or even the design of the manufacturing process of these. To this end, it is absolutely essential to know the constitutive flow rules of the materials if reliable and precise results are required in the FEM simulations. Several mathematical expressions which describe the relation between the stress and the strain for the material (σ =
f(ε)) may be found in the existing prior literature, such as those by Ludwik, Hollomon, Swift or Voce, among others [
1]. Each of them is fitted to the specific type of material and to its specific conditions: metallic alloy, ductile, brittle, low strain values and so on. Nevertheless, when simulating a part under different conditions of temperature and strain rate is needed, the previously-mentioned equations do not accurately predict the material behaviour.
There are other mathematical models which improve those above-mentioned ones, where the material stress is then obtained as a function of strain, strain rate and temperature, that is,
[
2,
3]. In Khan et al. [
4], the authors develop a material model which has taken into account the variation of mechanical strength, ductility and elastic modulus as a function of grain size, in such a way that the Hall-Petch equation is considered. This leads to the fact that it is ever-more necessary to know a large number of intrinsic material parameters. Moreover, in Yanagida et al. [
5], the effect of dynamic recrystallization is taken into consideration. The method proposed by these authors makes it possible to estimate the kinetics of dynamic recrystallization directly from the flow curve calculated by one hot compression test, thus avoiding having to carry out a large number of hot compression tests. In this way, the kinetics of dynamic recrystallization for plain carbon steel is successfully obtained for a wide range of strain rate at elevated temperature.
A new constitutive equation is also proposed by El Mehtedi et al. [
6] in order to model the hot working behaviour of an AA6063 aluminium alloy for a temperature range of 450 °C to 575 °C. This new proposed model is a modification of the Hensel-Spittel constitutive equation, where this has been modified by replacing the applied flow stress σ with the sinh(ασ) term, which was first introduced in the Garofalo equation [
7].
In Abbasi-Bani et al. [
8], a comparative study is made on two constitutive equations such as Johnson-Cook and Arrhenius type ones in order to describe the flow behaviour of Mg-6Al-1Zn, which is a magnesium alloy, under hot deformation conditions. Johnson-Cook equation is applied for metals and the flow stress is considered to be a function of temperature, strain rate and plastic strain [
8]. In the case of the Arrhenius equation, a sine hyperbolic description is selected because it is considered that the latter gives a better approach for the stress flow in a wide range of hot work regimes. From the comparison of the two studied models, it is concluded that the Arrhenius type one predicts the flow behaviour of the material quite well. On the contrary, the Johnson-Cook equation is inappropriate in this case, specifically at the softening stage.
Chen et al. [
9] propose four different constitutive models in order to study the hot deformation behaviour of a homogenized AA6026 aluminium alloy. The four equations proposed by these authors are as follows: original Johnson-Cook model, modified Johnson-Cook model, Arrhenius model and strain compensated Arrhenius model. This modelling is carried out by hot compression tests at temperature values that range from 400 °C to 550 °C. The comparison between the models is made in terms of both their correlation coefficient and their average absolute relative error. In this sense, the original Johnson-Cook model turns out to be inadequate to describe the flow behaviour but both the modified one as well as the Arrhenius model greatly improve the accuracy. Finally, the compensated Arrhenius model is the most accurate in predicting the hot flow stress for the AA6026.
In line with this, Mirzadeh [
10] also carries out a comparative study on the modelling and the prediction of hot flow stress using three different equations, which are as follows: Hollomon, Johnson-Cook and Arrhenius with strain compensation. The Johnson-Cook equation fails to predict the hot flow stress adequately. Another of the most important conclusions of this present research work is that the modelling by strain compensation should not consider the deformation activation energy as a strain function.
Guan et al. [
11] also propose two constitutive equations, one phenomenological and the other empirical, in order to describe the superplastic behaviour of an Al-Zn-Mg-Zr alloy. Both equations are developed from tension tests at a temperature of 530 °C and they are compared with the experimental data. One of the main conclusions is that the empirical equation has a higher level of accuracy than the phenomenological equation when predicting the flow behaviour in a wide range of strain and strain rate values.
In line with this above-mentioned study, the research work by Shamsolhodaei et al. [
12] is also noteworthy. In this case, these researchers make a comparative analysis on a modified Zerilli-Armstrong model and an Arrhenius type constitutive model, where the aim is to examine the constitutive flow behaviour of a near equiatomic nickel titanium shape memory alloy at a temperature range of 700 °C to 1100 °C. The most important conclusion is that the Arrhenius equation may represent the flow behaviour at this temperature range more accurately than the modified Zerilli-Armstrong model. Nevertheless, this last model requires a shorter amount of time in order to evaluate the material constants for the model.
Mirzaie et al. [
13] evaluate the application of the Zerilli-Armstrong model to the modelling and prediction of hot deformation flow stress. The study is carried out in the case of face-centred cubic materials and one of their most important conclusions is that the original model needs to be modified by considering the peak strain into the flow stress formula and both hardening and softening phenomena.
Gan et al. [
14] develop a constitutive equation in order to model the flow stress of an AA6063 aluminium alloy at high temperature. To this end, isothermal compression tests are performed at temperature values which vary from 300 °C to 500 °C. In their research work, the authors use the Zener-Hollomon parameter in an exponent type equation to describe the combined effects of temperature and strain rate on the hot flow behaviour of the material. They come to the conclusion that the thus-developed constitutive equation can accurately predict the flow stress for AA6063 at high temperature and then, it may be used in the simulation of hot deformation processes, such as extrusion and forging.
One limitation when obtaining material constitutive equations is the range of the plastic strain values these flow rules are to be applied to. In most cases, the proposed material flow rule reaches up to strain values of ε ≈ 0.2, where this value is much lower than that achieved when the forming of a part with high deformation values is being studied [
15]. This occurs in most cases of forging processes. In the case of the ECAP (Equal Channel Angular Pressing) process [
16,
17], which is a severe plastic deformation (SPD) process, León et al. [
18] apply a methodology which allows the σ =
f(ε) relationship to be obtained for experimental strain values of ε ≈ 4 or even higher but in any case, not taking into account different compression velocities or work temperatures. Another way to obtain flow rules with high plastic strain values is through the use of compression tests, as outlined in Fereshteh-Saniee et al. [
19] and Wu et al. [
20]. In line with the modelling of SPD processes, Mohebbi et al. [
21] study the flow behaviour of AA1050 sheets processed by Accumulative Roll Bonding (ARB). In order to do this, plane strain compression tests along with stress relaxation are carried out, where the existence of dynamic recovery is validated with the help of transmission electron microscopy (TEM).
Another commonly-used methodology over the last decade is the employment of artificial neural networks (ANN) in order to predict the behaviour of a material from experimental compression tests. Thus, Ping et al. [
22] train a neural network so as to precisely predict the stress strain curve of a Ti-15-3 alloy, thus avoiding the problem of solving empirical equations with a large number of parameters. These authors use the predicted stress strain curves in order to optimise the hot forging processing conditions for the above-mentioned titanium alloy.
In the research work by Salcedo et al. [
23], the authors make a comparison on the mechanical properties of an AA6063 aluminium alloy after having been processed by isothermal forging starting from different initial strain values for the material. In order to model the material behaviour, the flow rule is determined by means of neural networks, where the equations obtained are of type σ = σ(ε,
T). With this methodology, Salcedo et al. [
23] are able to predict the material behaviour up to strain values of ε = 2.5.
In relation to the use of artificial neural network for modelling, the research study by Sabokpa et al. [
24] also merits a mention. Here the ANN methodology is employed to describe the high temperature flow behaviour of a cast AZ81 magnesium alloy. Isothermal compression tests are carried out for a temperature range of 250 °C to 400 °C up to a plastic strain value of ε = 0.6. Constitutive equations of Arrhenius type are compared with a feed forward back propagation artificial neural network model, this evaluation indicating that the trained ANN model is not only more accurate but also more efficient in terms of the calculation of the large number of material constants associated with constitutive equations.
In this present research work, a new constitutive law is proposed in order to predict the material behaviour. The procedure proposed is valid for predicting the behaviour of both materials with no previous strain and materials which have been previously subjected to severe plastic deformation (SPD) processes. In addition, the procedure is also valid to model both the hot and the cold behaviour of such materials, as will be shown below in
Section 5. In this study, experimental tests are carried out over a temperature range which varies from room temperature to 300 °C (specifically, 25 °C, 100 °C, 150 °C, 200 °C, 250 °C, and 300 °C) and by considering materials with or without any previous deformation accumulated. The results obtained are compared with the values predicted by other laws from the existing prior literature. Moreover, the study is complemented with finite element simulations and with a comparison between the results obtained from the simulations and those experimentally obtained. With the new flow rule proposed in this research work, it is possible to take the recrystallization phenomenon into account. Furthermore, this study presents a novelty in that the methodology here proposed is valid for both hot and cold working.
2. Proposal of a Novel Flow Rule
The new formulation proposed in this present research work is shown in Equation (1). It presents a novelty in that it is able to predict the behaviour of the materials analysed adequately over a wide temperature range, which varies from 25 °C to 300 °C. In order to demonstrate that the new flow rule proposed has an accurate fitting in relation to the true material behaviour, it is compared with two of the most widely-used flow rules: Hollomon (Equation (2)) and Voce (Equation (3)). As was previously mentioned, this equation was validated not only for materials in an annealed state (N0) but also for materials which present a certain amount of accumulated strain because they have been subjected to some kind of plastic deformation process. Specifically, it is the materials that have been previously ECAP-processed which were taken as a starting point.
Now, the procedure followed in order to obtain the five constitutive parameters of this new law (σ
2,
a,
b,
c,
d and
n) is explained. In order to select the parameters shown in Equation (1), firstly, it is necessary to observe the behaviour of the flow rule (see
Figure 1).
Figure 1a shows how the stress-strain curve is always growing and it is named as Case 1. This case takes place when the material is usually deformed at room temperature. Adversely, as can be observed in
Figure 1b, the stress-strain curve reaches a maximum value of stress after which it decreases down to a point in which it may remain constant or a new increase of the curve may even occur. This second case is named as Case 2. In general, this second case corresponds with materials which are deformed at high temperature. In function of the type of curve that takes place, Case 1 or Case 2, the parameters are calculated.
The first step is to locate the three characteristic points from the experimental curve (S
0, S
1 and S
2), which are shown in
Figure 1. Each of these three points is obtained depending on the study case (Case 1 or Case 2). When Case 1 is selected, the S
0 point corresponds with the yielding point. The S
1 point corresponds with the point where the change in the curve tendency occurs and it is calculated from the third derivative of the stress-strain curve. This derivative may be numerically calculated from the experimental data of the stress-strain curves.
Figure 2 shows the stress-strain curve as well as its three corresponding derivatives: first, second and third. Finally, the S
2 point is taken before barrelling happens, which may be determined since a variation in the slope of the stress-strain curve is observed.
For Case 2, the procedure followed in order to calculate the three characteristic points (S0, S1 and S2) is outlined as follows. The S0 point continues to be the yielding point. The S1 point corresponds with the stress highest value and the S2 with the lowest one before a change in the tendency for the stress-strain curve is observed, where this means that the curve starts to increase and therefore, barrelling begins.
Once the characteristic points are obtained from the stress-strain curve, the above-mentioned parameters (a, b, c, d and n) from Equation (1) are calculated. In order to calculate d and n, a least squares fitting is performed taking different values of these within a specific interval. The range of values for d and n that is considered in this study is as follows: 0.1 ≤ d ≤ 3.1 and 0.1 ≤ n ≤ 10.1. The variation in d and n is carried out by increments of 0.1. It may be pointed out that the ranges considered can be widened with no loss of generality. For each pair of values d and n, values for a, b and c are found by solving a system of equations which establishes that the flow rule considered has to pass through the previously selected three characteristic points (S0, S1 and S2). The pair of values for d and n which gives the lowest value for the mean squared error (as Equation (4) shows) will be that selected and, therefore, the parameters a, b and c from Equation (1) will also be determined.
The values of the parameters obtained for the new flow rule proposed in the case of the three aluminium alloys considered in this present study and along with their different work intervals will be shown in the sections below.
3. Experimental Set-Up
In order to validate the flow rule proposed in this present research work, an experimental study of three aluminium alloys with different starting states is to be carried out, as is shown in
Table 1. The aluminium alloys under consideration are as follows: AA5754, AA3103 and AA5083. The different starting states studied are: N0, N2, N2 + TT and N4. N0 represents the annealed initial state. In the N2 state, the aluminium alloy is subjected to a severe plastic deformation using the ECAP (Equal Channel Angular Pressing) process which consists of two passages with route C (where this means rotating the billet 180° after each ECAP passage). In the N2 + TT state, after having ECAP-processed the material twice, a recovery thermal treatment is carried out, where this consists in heating the material up to 300 °C with a heating slope of 12 °C/min. In the N4 state, the aluminium alloy undergoes four ECAP passages with route Bc, which consists in rotating the billet 180°, 90° and 180°, respectively, after each passage.
In order to perform the experimental tests, a universal testing machine, which belongs to the Public University of Navarre, Spain (see
Figure 3a), was used. To carry out the compression test of the billets, it was necessary to design and to manufacture the plane-shape dies shown in
Figure 3b. Finally, a set of electrical resistances was employed in order to be able to perform the compression tests isothermally. With the aim of increasing the heat flow towards the inner part of the dies and, specifically, at the contact zone with the test billet, a steel sheet was attached to the outer part of the resistances, as can be observed in
Figure 3c. In this way, it is possible to carry out tests at higher temperature values, as the convection heat evacuation outwards is lower. The resistance temperature is controlled by two thermocouples, which are in contact with each of the two dies. This temperature value may be modified by the two digital PDI controllers, which are shown in
Figure 3c. Once the test temperature for the resistances is programmed, the billet temperature is monitored by an external digital thermocouple until the moment when the billet reaches the desired stationary temperature value. Subsequently, the billet is placed in the proper position and it is maintained between the plane-shape dies, as is shown in
Figure 3b, until it reaches the same temperature value as the dies and thus, to perform an isothermal compression test. By following this procedure it means there is a higher degree of test temperature control during the deformation process.
In order to obtain the load-stroke curve for all the aluminium alloys under consideration, six different test temperature values (25 °C, 100 °C, 150 °C, 200 °C, 250 °C and 300 °C) were studied at a compression velocity of 60 mm/min. In order to carry out this experimental study, cylindrical billets with a diameter of 8 mm and a length of 16 mm were employed. These billets are placed on the lower die and both the billet and the die are lubricated with Teflon®. In the case of the tests at a temperature higher than 25 °C, the billet is maintained on the lower die during five minutes so that the billet may reach the same temperature as the dies and thus, the compression is isothermal.
Once the compression tests for the different aluminium alloys are carried out at the different temperature values selected, the load-stroke curve obtained for each of them is transformed into a stress-strain curve, as a previous step to the attainment of, on the one hand, Hollomon’s and Voce’s laws and, on the other hand, the new flow rule proposed in this research work.
Figure 4a,c,e shows the complete stress-strain curves for AA3103 taking into account the barrelling effect that the billet undergoes when it is compressed, whereas
Figure 4b,d,f shows the selected zone of the experimental curves. From this selected zone, the three flow rules under consideration are assessed. The cutting zone is established at a strain value of 0.7 although there are several cases for the ECAP-processed material where this cutting zone is lower than this former value. The reason for this is due to the change in the existing slope in the stress-strain curve for the N2 and N4 starting materials. This slope change is because of the recrystallization process, which takes place in the material at some specific temperature values, as well as the barrelling effect due to friction.
Figure 5a,c,e shows the stress-strain curves obtained for AA5754 from the experimental compression tests between plane-shape dies, whereas
Figure 5b,d,f shows the zones selected in order to assess the three above-mentioned flow rules in all the cases of temperatures and starting states under consideration. In addition, it may be observed that the curves at the three lowest temperature values (that is, 25 °C, 100 °C, and 150 °C) are closer to each other than the other three remaining ones (200 °C, 250 °C, and 300 °C). This is due to the temperature effect on the material strength, which is more significant above 150 °C.
Finally, the stress-strain curves for AA5083 at its three starting states can be observed in
Figure 6a,c,e. Just as occurs with AA5754, the curves at 25 °C and at 100 °C show a similar behaviour, whereas, from a temperature of 150 °C, a change in their behaviour is observed since the stress-strain curves are found to be more spaced out.
From the flow curves shown in
Figure 4,
Figure 5 and
Figure 6 and following the previously-mentioned procedure, it is possible to obtain the parameters to be considered in Equation (1). These parameters are shown in
Table 2,
Table 3 and
Table 4 for AA3103, AA5754 and AA5083, respectively, where these alloys may be without any previous deformation or after having been previously ECAP-processed.
5. Discussion of Results
Once the experimental procedures used have been shown, the results obtained in the present research work are now outlined in this section. First of all,
Figure 7 shows the approach which the three flow rules present both at low (25 °C) and at high (250 °C) temperature in relation to the experimental stress-strain curve. It may be observed that, at low temperature, the approach is acceptable for Hollomon’s and Voce’s laws and it is very good for the new law proposed, as it is practically superimposed on the experimental curve. Nevertheless, at a temperature of 250 °C, where dynamic recrystallization takes place for the previously ECAP-processed material, the curve corresponding with Hollomon’s law increases up to a stress value of 240 MPa (about 40 MPa above the experimental curve) and then undergoes a sudden decrease which crosses the experimental curve at about a strain value of 0.04. On the other hand, the curve corresponding with Voce’s flow rule reaches the same maximum stress value as that from the experimental curve but then remains constant as strain increases, where this means that the approach seems not to be good enough either, as will be numerically demonstrated below. However, although the new law proposed does not show such a very good approach as found at low temperature values, the approach along the decreasing part is good enough. Therefore, it may be concluded that the approach of this new law to reality is very good and, thus, much better than in the cases of both Hollomon’s and Voce’s laws.
With the aim of quantifying the precision of fit between the different flow rules studied here and the true behaviour of the material during the compression test between plane-shape dies, the root mean square error (RMSE) was assessed in each of the experiments with the three aluminium alloys under consideration, as is depicted in Equation (5).
For AA3103, as can be observed in
Figure 8, the flow rule proposed fits the behaviour of the alloy better in the temperature range considered. Moreover, it is observed that Voce’s flow rule is that with the highest values for RMSE. For the material with N0 state (see
Figure 8a,b), the error fluctuates between positive and negative values, whereas with the ECAP-processed material (see
Figure 8c–f), the error value is negative and it increases as the material is compressed, where this tendency turns out to be greater at higher temperature values. In the case of Hollomon’s flow rule, it is observed that there is a fitting which is better than in the case of Voce’s flow rule. For N0 starting state, the error value is negative at first and then turns positive and decreases as compression progresses. For N2 and N4 state starting materials (see
Figure 8c–f), the error value is positive at first until it turns to negative around a strain value of 0.35. In the case of the new flow rule and for N0 state material, the error with the highest value occurs at the beginning of the compression process, whereas, as the material is being compressed, the error value decreases. For N2 and N4 starting states (
Figure 8c–f), the error obtained with the new flow rule is clearly smaller than in the case of the two remaining ones. It may be stated that its approach to the experimental curve is very good, where this turns out to be better at lower temperature values, as can be observed in
Figure 8c,e.
Figure 9 shows SEM (Scanning Electron Microscopy) micrographs for billets forged at 250 °C from different deformation states. As can be observed in N0 material, the grain size is bigger than in the two other cases due to a lower level of strain accumulation.
If RMSE graphs for AA5754 are analysed at the three starting states under consideration (N0, N2 and N2 + flash), it may be observed that the curve with the lowest error value (that is, with the best fitting) is that corresponding to the new flow rule (see
Figure 10). In the case of Hollomon’s flow rule, it may be observed in
Figure 10a that for N0 state at low temperature, the error shows a behaviour which is the contrary to that of Voce’s flow rule. For N2 starting state at low and high temperature (see
Figure 10c,d) as well as for N2 + flash starting state at high temperature (see
Figure 10f), the error graph from Hollomon’s flow rule is similar. Nevertheless, for N2 + flash starting state at low temperature (see
Figure 10e), the curve is similar to that from N0 at low temperature (see
Figure 10a). In the case of Voce’s flow rule, it is observed in
Figure 10d that for N2 starting material at high temperature, the error graph shows the same behaviour as its equivalent from AA3103 (
Figure 8d), which, at the same time, is contrary to the rest of the starting states and temperatures. This is due to the fact that the temperature value coincides with the recrystallization temperature of AA5754 for a N2 starting state. In the case of the new flow rule proposed, it may be observed that the fit to the experimental curve is better.
Figure 11 shows SEM micrographs for AA5754 billets forged at 250 °C from different deformation states. It is observed that recrystallization occurs in the previously ECAP-processed samples, as a consequence of the higher level of previous strain accumulated. Finally, with regard to the error graphs obtained from AA5083 (see
Figure 12), a similar behaviour to AA5754 may be observed.
Figure 13 shows SEM micrographs for AA5083 billets forged at 250 °C from different deformation states. As in the previous case, it is observed that for previously ECAP-processed billets, recrystallization takes place because of the higher value of previous strain.
Table 11 shows the error values obtained for AA3103. It may be observed that for N0 starting state at low temperature, the best approach turns out to be for Hollomon’s flow rule followed by the new flow rule at a close distance, whereas at high temperature, the new flow rule turns out to be the best. With regard to N2 and N4 starting materials, the best approach is achieved with the new flow rule and with a significant difference in relation to Hollomon’s and Voce’s flow rules.
In the case of the error values for AA5083 (see
Table 12), it is observed that the new flow rule provides a better approach for all the starting states in comparison with Hollomon and Voce. The only exception occurs in the case of N2 starting material at 25 °C and at 100 °C, where Voce’s flow rule approaches better. Lastly, it should be pointed out that for this aluminium alloy, the difference in the approach between the new flow rule proposed and those from Hollomon or Voce is rather significant, as can be observed in
Table 12.
Similar to what happens with AA5083, the flow rule with the best approach to the experimental curve for AA5754 is the new flow rule in all cases (see
Table 13). There is only one exception and it occurs with N2 starting material at 25 °C, where the best approach corresponds with Voce’s flow rule. As for AA5083, in the rest of the cases, the improvement achieved by the new flow rule in the approach to the experimental curve is rather significant in relation to Hollomon’s and Voce’s flow rules.
6. Validation by FEM
Now that the procedure followed in order to obtain the different material flow rules and the results obtained in the comparison with the experimental curve have previously been outlined, the same study is carried out by finite element modelling (FEM) with the aim of comparing the fitting of the three flow rules under consideration.
In order to do this, a model, which consists of the two plane-shape dies and the billet, is generated. A fine meshing for the billet is performed with an element size in height of 0.4 mm and a hexahedral type element with six control nodes, as can be observed in
Figure 14. After this, the material flow rules with the elastic part also being included are introduced with a Young’s modulus for aluminium of 70 GPa and a Poisson’s coefficient of 0.3. Furthermore, it is necessary to indicate the experimental test temperature that is required for the simulation. Therefore, fifty-four FEM simulations are needed in order to cover all the experimental study, which analyses six temperatures, three flow rules and three starting states.
Having explained the procedure to study the process by finite element simulations, the load-stroke curves obtained both for each flow rule and for each aluminium alloy are shown below. In the case of AA3103, as is shown in
Figure 15, at low temperature values, the approach of the three flow rules to the curve experimentally obtained is rather good. Only the curve from Voce’s law differs from the experimental curve as the process progresses, whereas in the cases of Hollomon’s law and the new law proposed, they follow the experimental curve very satisfactorily. Nevertheless, in the case of the starting states where the material is ECAP-processed, the approach of Voce’s law is worse and the curve with the best approach is that from the proposed new flow rule. The curve from Hollomon’s law also has a good approach but nevertheless, at the initial part, it differs from the experimental curve more than the curve, which corresponds with the new law (see
Figure 15c–f).
Figure 16 shows the load-stroke curves obtained from the FEM simulations for AA5754 at the three starting states and at temperatures of 25 °C and 250 °C. It may be observed that the approach of the proposed new law is very good until the compression final part, which is the zone with the highest error value. This also happens in the cases of Hollomon’s and Voce’s laws and it is due to the barrelling effect which takes place during the compression process. For this aluminium alloy, the approach of the three laws is acceptable and only in the cases of Hollomon’s law at low temperature (see
Figure 16a) and Voce’s law at high temperature and N2 state (see
Figure 16d), they present appreciable differences in relation to the experimental curve. Similarly for the AA3103, the curve with the best approach to the experimental one is the new flow rule.
In the case of AA5083 (see
Figure 17), the approach of the load-stroke curves for each flow rule is very similar to which was obtained for AA5754.
Figure 17a shows a deviation of Hollomon’s law at low temperature with respect to the experimental curve, as occurs in
Figure 17e. In addition, in the case of the material ECAP-processed twice and at high temperature, it may be observed that Voce’s law differs from the rest practically from the beginning of the compression process.
In the same way as for the experimental part, the error of each of the curves obtained from FEM was quantified for all the flow rules studied in relation to the experimental curve. As can be observed in
Table 14, in most cases, the curve with the best approach for AA3103 is that which corresponds with the new flow rule proposed in this research work. Only in three cases at low temperature values, does the curve from Hollomon’s law have a better approach than the new law.
For AA5083 (
Table 15), the new flow rule proposed shows a better approach in the FEM results practically in all cases, except for the case at 100 °C with the N2 starting state. Moreover, it may be appreciated that the error value in the cases of Hollomon’s and Voce’s laws is about ten times more than that obtained for the new law proposed in some cases.
As occurs in the two previous aluminium alloys, the flow rule with the best approach to the experimental curve for AA5754 turns out to be the new flow rule in all cases, as can be observed in
Table 16. There are only three exceptions: the two cases with a N2 starting material at 25 °C and at 100 °C, where the best approach is for Voce’s curve, and the case with a N4 starting material at 200 °C, where the curve which has the best approach to the experimental one is that from Hollomon. In the rest of the cases, the approach of the new law to the experimental curve is superior to Hollomon’s and Voce’s laws.