Study on Constitutive Relation of Nickel-Base Superalloy Inconel 718 Based on Long Short Term Memory Recurrent Neural Network

: The high temperature tensile test of Inconel 718 under the conditions of deformation temperature of 950 ◦ C–1100 ◦ C and strain rate of 0.0005 s − 1 –0.1 s − 1 was carried out, and its true stress–true strain curve was drawn. Through the analysis of the flow stress of Inconel 718 under different conditions, it can be seen that the high-temperature rheological behavior of Inconel 718 is affected by the coupling of strain hardening effect and dynamic softening effect, and has significant loading history correlation. By applying the stretched data, a long short term memory (LSTM) recurrent neural network was trained to characterize the constitutive relationship of Inconel 718. The experimental results show that the prediction results of the LSTM constitutive model are extremely consistent with the experimental data, which is significantly better than the modified Johnson–Cook (M-JC) model. Finally, high temperature tensile experiments under variable strain rates were carried out to verify the feasibility of the LSTM constitutive model in the complex loading and unloading stages.


Introduction
Inconel 718 is currently the most widely used nickel-based superalloy, it has good comprehensive mechanical properties and mechanical stability below 650 • C [1,2]. The main strengthening mechanism of Inconel 718 is precipitation strengthening. The main strengthening phase is γ (Ni3Nb), supplemented by the weak strengthening phase γ (Ni3AlTi) with a face-centered cubic structure [3]. The γ phase is metastable, and its equilibrium phase is the δ phase (Ni3Nb) with orthogonal structure [4]. Compared with other types of precipitation hardening alloys, Inconel 718 has a larger coherent distortion between the precipitation strengthening phase γ phase and the matrix, which makes the structure and properties of the alloy more stable [5,6]. Due to its good high-temperature strength, oxidation resistance, thermal corrosion resistance and fatigue performance, Inconel 718 is widely used in the manufacture of hot parts of aero engines and high temperature parts of aerospace rocket engines [7]. Meanwhile, Inconel 718 is also recognized as a difficult-to-form material. It has extremely high resistance to high temperature deformation, and its alloy structure is particularly sensitive to the forming process. In most hot forming and welding processes, Inconel 718 has the problems of high process temperature requirement and narrow selection of process parameters [8]. At present, in order to save material and experimental costs, finite element analysis of the forming process has become the main method to evaluate the formability of Inconel 718. The reliability of the finite element method mainly depends on the constitutive model that can characterize the dynamic rheological properties of the material [9].
In the deformation process of Inconel 718, there are two competing mechanisms: The work hardening and the softening induced by dynamic recovery and dynamic recovery stabilization [10]. These two mechanisms are related to strain, strain rate and temperature, and the constitutive model can characterize the correlation. In review of high temperature constitutive models for metals and alloys, Lin divided the current constitutive models into three categories: Phenomenological constitutive model, physical based constitutive model and artificial neural network (ANN) constitutive model [11]. The basis of the phenomenological constitutive model is the experimental observation of the flow stress under different loading conditions, it does not have strict physical meaning. This type of model is widely used in finite element analysis (FEA) of metal forming process because it contains fewer material parameters and is easy to measure [12]. In order to facilitate the solution of material constants, early studies on the phenomenological constitutive model of Inconel 718 assumed that the effects of strain, strain rate and temperature on rheological behavior are independent, which significantly reduces the prediction accuracy of the model [13][14][15]. In order to study the coupling effect of strain rate and temperature, Wang The weight of each function is solved by double multivariate nonlinear regression [17]. Lin et al. established a modified Arrhenius model of Inconel 718 taking the effects of multi-factor coupling into consideration. The material constants in this model are the function of the Zener-Hollomon parameter. These studies showed that the coupling effect between strain, strain rate, temperature and other parameters is nonlinear, and it is inaccurate to use traditional regression methods to obtain material constants. Therefore, many researchers used neural network methods to analyze the rheological properties of materials to overcome the limitations of traditional regression methods in solving highly nonlinear problems. Huang and Blackwell measured the flow stress data of Inconel 718 at a temperature range of 950 • C to 980 • C and strain rates of 10 −4 s −1 and 10 −3 s −1 , and used the data to train an ANN constitutive model. Although this model is limited by a small sample size, it has potential in characterizing material work hardening and dynamic softening phenomena [18]. Yao et al. conducted a single impact compression test of Inconel 718 at a temperature of 950-1050 • C and a strain rate of 0.004-10 s −1 . The prediction ability of the single hidden layer BP (Back propagation) neural network model based on the experimental results is significantly better than the modified Arrhenius constitutive model [19]. Bhekisipho Twala introduced the aging time of Inconel 718 into the input variable to characterize the microstructure of the material. The ANN constitutive model trained by this four-dimensional input data set can accurately predict the flow stress of Inconel 718 in multiple initial microstructure states [20].
With the promotion of the phenomenological constitutive model in FEA of various processes, some complex porous lattice structure parts in the aerospace field have also begun to be manufactured using Inconel 718, such as honeycomb structure, lattice truss structure and corrugated board structure [21,22]. In order to meet the lightweight standards in the aerospace field and some special heat dissipation requirements, these structures are usually designed with complex flow channels and micro-support structures, and these structures may become unstable if they undergo slight thermal deformation. In order to weigh the forming quality and the degree of deformation in the hot forming process, modern hot forming processes use complex loading paths to achieve multiple process target control. For example, the loading process of the diffusion bonding process of Inconel 718 includes multiple stages such as plastic deformation, creep and stress relaxation. As mentioned above, the current phenomenological constitutive model and ANN constitutive model of Inconel 718 are mostly based on the total strain theory of plasticity. They mainly focus on the flow stress σ under Metals 2020, 10, 1588 3 of 14 the coupling action of strain ε, temperature T and strain rate . ε. The corresponding relationship can be expressed as: Obviously, these models are only suitable for simple loading situations where the strain component changes proportionally, ignoring the deformation history effect. However, Inconel 718 is a typical viscoplastic material at high temperatures, and its rheological properties are obviously affected by the loading history. Under the complex loading path, the temperature T and strain rate . ε should be expressed as time-dependent dynamic functions T(t) and . ε(t), then the flow stress at any time t 0 can be expressed as follows: Ψ is an internal variable function which characterizes the microstructure change of material from the beginning of loading to the time t 0 , which is related to the integral of strain, strain rate and temperature along the loading curve. If the time increment on the loading path is assumed to be ∆t, then at a certain time t 0 the flow stress is not only related to the loading conditions at that time, but also to the loading conditions of (t 0 − ∆t). It is difficult for ANN model to obtain satisfactory results for the prediction of time-dependent sequence data. This is because there is no correlation between the inputs of ANN model, and its outputs are independent of each other. The recurrent neural network (RNN) model is a kind of neural network which takes sequence data as input and recurses in the evolution direction of sequence. Different from ANN model, the hidden layer nodes of the model are connected, and the input of the hidden layer is the superposition of the output of the current input layer and the output of the previous hidden layer [23]. In other words, when RNN model calculates the output of a certain time on the time series, it not only considers the input of the time point, but also considers the memory of the output of the previous time point. Obviously, many RNN models have advantages in learning nonlinear features of sequences. Among them, the long short term memory (LSTM) recurrent neural network, as a specially designed RNN model, effectively avoids the long-term dependency problem of RNN model, so it is widely used in the prediction of various time series [24].
In this study, several groups of Inconel 718 were created and their uniaxial tensile tests at different temperatures and strain rates were carried out, and a LSTM constitutive model was trained based on the experimental results. Compared with the prediction results of modified Johnson-Cook (M-JC) model obtained by experimental parameters fitting, the root mean square error (RMSE) of rheological properties prediction results of the two models under different loading conditions were obtained. In addition, based on the hot forming process curve of Inconel 718 in project, a high temperature tensile test under a complex loading path was designed. The obtained experimental data is compared with the prediction results of LSTM constitutive model under the same path, and the possibility of LSTM constitutive model applied to the flow stress prediction of Inconel 718 under complex loading path was studied.

Experiments
Commercial hot-rolled Inconel 718 sheet material with a thickness of 3 mm was selected as the blank of the high temperature tensile test specimen, which was purchased from Sunrise Metal Technology (Beijing, China) Co. Ltd. The material element composition and ratio are shown in Table 1. Wire cut electrical discharge machining (EDM) equipment was used to prepare tensile test pieces, the dimensions of these specimens are shown in Figure 1, where the axial direction is the rolling direction of the sheet. To prevent Inconel 718 from being oxidized at high temperature, zirconia powder was applied to the sample surface. min after the three temperature sensors reached the set temperature, and then the specimen was loaded at a constant strain rate. After the experiment, it was taken out when it cooled to room temperature.  Figure 1. Schematic diagram of specimen size for tensile test (unit: mm).

Flow Stress Curve Analysis of Inconel 718
Through the high temperature tensile experiments of Inconel 718, engineering stressengineering strain curves under different temperature and strain rate conditions were obtained. For the convenience of research, equations (3) and (4) are usually used to convert the engineering stress and engineering strain of the material into true stress and true strain , respectively.
Many studies have pointed out that this method has low accuracy in predicting the rheological behavior of materials after necking. According to previous research results [25], an FE-based inverse modeling procedure was used in this paper to derive the true stress-strain curve of Inconel 718 as shown in Figure 2.
It can be seen from the Figure 2 that the stress-strain curves of Inconel 718 at various temperatures have a close changing trend. In the initial stage of deformation, the flow stress increases significantly with the increase of strain, and the work hardening effect caused by the increase of dislocation density is significant. As the deformation continues, the dynamic softening effect inside the material gradually increases, and the dislocation density begins to show a downward trend. The two effects reach equilibrium at a certain deformation moment, at which time the flow stress will tend to stabilize. Obviously, the level of flow stress is affected by the coupling of temperature and strain rate. Under the same temperature conditions, the deformation resistance of Inconel 718 decreases significantly as the strain rate decreases. This is because the deformation time is shortened and the nucleation number of dynamic recrystallization is reduced under high strain rate, which leads to insufficient dynamic softening process of the material. On the other hand, the increase of temperature promotes the slip and climb of dislocations, thereby it improves the nucleation rate of dynamic recrystallization and promotes the dynamic softening process of the material. This also explains the phenomenon that the flow stress increases significantly with the decrease of temperature under the same strain rate. The high temperature tensile testing machine (Strong Association Instrument Technology Co. Ltd, Shanghai, China) was used for the test, and the maximum strain obtained was 0.6. The experimental temperature was set at 900 • C, 950 • C, 1000 • C, 1050 • C and 1100 • C, and the strain rate of 0.1 s −1 , 0.01 s −1 , 0.001 s −1 and 0.0005 s −1 . The experimental equipment adopts three-stage temperature measurement technology, which can effectively ensure that the parts are heated evenly. The temperature was increased at a rate of 20 • C/min firstly, and the temperature was kept for 10 min after the three temperature sensors reached the set temperature, and then the specimen was loaded at a constant strain rate. After the experiment, it was taken out when it cooled to room temperature.

Flow Stress Curve Analysis of Inconel 718
Through the high temperature tensile experiments of Inconel 718, engineering stress-engineering strain curves under different temperature and strain rate conditions were obtained. For the convenience of research, equations (3) and (4) are usually used to convert the engineering stress σ e and engineering strain ε e of the material into true stress σ and true strain ε, respectively.
Many studies have pointed out that this method has low accuracy in predicting the rheological behavior of materials after necking. According to previous research results [25], an FE-based inverse modeling procedure was used in this paper to derive the true stress-strain curve of Inconel 718 as shown in Figure 2.
It can be seen from the Figure 2 that the stress-strain curves of Inconel 718 at various temperatures have a close changing trend. In the initial stage of deformation, the flow stress increases significantly with the increase of strain, and the work hardening effect caused by the increase of dislocation density is significant. As the deformation continues, the dynamic softening effect inside the material gradually increases, and the dislocation density begins to show a downward trend. The two effects reach equilibrium at a certain deformation moment, at which time the flow stress will tend to stabilize. Obviously, the level of flow stress is affected by the coupling of temperature and strain rate. Under the same temperature conditions, the deformation resistance of Inconel 718 decreases significantly as the strain rate decreases. This is because the deformation time is shortened and the nucleation number of dynamic recrystallization is reduced under high strain rate, which leads to insufficient dynamic softening process of the material. On the other hand, the increase of temperature promotes the slip and climb of dislocations, thereby it improves the nucleation rate of dynamic Obviously, due to the complex superposition state of the material's internal hardening effect and dynamic softening effect in the time domain, Inconel 718 exhibits different deformation resistance levels under different deformation conditions, which can be expressed as a set of independent continuous curves. When the actual loading curve changes at a certain moment, the flow stress not only steps on different continuous curves, but also the influence of the superposition state of the material at that moment should be considered. Therefore, it is necessary to establish a constitutive model considering the loading history effect based on the tensile experimental data under different deformation conditions to predict the flow stress at any time in the complex loading curve.

The Structure of LSTM Constitutive Model
The cyclic neural network is a fully connected neural network with a feedback mechanism. There are node connections not only between layers, but also across the timeline. Its structure is shown in Figure 3. The value of the hidden layer neuron of the RNN at a certain moment will be stored inside the neuron, which is called memory ℎ . When calculating the output of the neuron at the next moment, the memory will work together with the input at that moment. It can be expressed as: Obviously, due to the complex superposition state of the material's internal hardening effect and dynamic softening effect in the time domain, Inconel 718 exhibits different deformation resistance levels under different deformation conditions, which can be expressed as a set of independent continuous curves. When the actual loading curve changes at a certain moment, the flow stress not only steps on different continuous curves, but also the influence of the superposition state of the material at that moment should be considered. Therefore, it is necessary to establish a constitutive model considering the loading history effect based on the tensile experimental data under different deformation conditions to predict the flow stress at any time in the complex loading curve.

The Structure of LSTM Constitutive Model
The cyclic neural network is a fully connected neural network with a feedback mechanism. There are node connections not only between layers, but also across the timeline. Its structure is shown in Figure 3. Obviously, due to the complex superposition state of the material's internal hardening effect and dynamic softening effect in the time domain, Inconel 718 exhibits different deformation resistance levels under different deformation conditions, which can be expressed as a set of independent continuous curves. When the actual loading curve changes at a certain moment, the flow stress not only steps on different continuous curves, but also the influence of the superposition state of the material at that moment should be considered. Therefore, it is necessary to establish a constitutive model considering the loading history effect based on the tensile experimental data under different deformation conditions to predict the flow stress at any time in the complex loading curve.

The Structure of LSTM Constitutive Model
The cyclic neural network is a fully connected neural network with a feedback mechanism. There are node connections not only between layers, but also across the timeline. Its structure is shown in Figure 3. The value of the hidden layer neuron of the RNN at a certain moment will be stored inside the neuron, which is called memory ℎ . When calculating the output of the neuron at the next moment, the memory will work together with the input at that moment. It can be expressed as: The value of the hidden layer neuron of the RNN at a certain moment will be stored inside the neuron, which is called memory h t . When calculating the output y t+1 of the neuron at the next moment, the memory will work together with the input x t+1 at that moment. It can be expressed as: y t+1 = Uh t+1 (6) in the equation, W is the input transfer function, V is the memory transfer function, U is the output transfer function and ϕ is the activation function. The memory is updated with iteration, but it is assigned different weights. Since the output of the neuron is obtained through the activation of function, whenever the memory is iterated with other inputs, part of the information it contains will be lost. In fact, after four to five iterations, all the memory information will be lost. This is the long-term dependencies problem of the RNN model [24]. For the complex loading curve of Inconel 718, a sudden change and recovery of conditions such as temperature and strain rate occur at a certain moment. When the RNN model is used to forecast the flow stress and if the prediction time is far from the sudden change time, the memory information of the sudden change is lost. The historical effect of this mutation is ignored.
The LSTM network is a cyclic neural network with a special structure design, which is specialized in the ability to remember long-term information, so as to avoid the loss of information in the propagation process. The LSTM network is also composed of multiple groups of repeated neuron modules, and its internal structure is shown in Figure 4.
in the equation, is the input transfer function, is the memory transfer function, is the output transfer function and is the activation function. The memory is updated with iteration, but it is assigned different weights. Since the output of the neuron is obtained through the activation of function, whenever the memory is iterated with other inputs, part of the information it contains will be lost. In fact, after four to five iterations, all the memory information will be lost. This is the long-term dependencies problem of the RNN model [24].
For the complex loading curve of Inconel 718, a sudden change and recovery of conditions such as temperature and strain rate occur at a certain moment. When the RNN model is used to forecast the flow stress and if the prediction time is far from the sudden change time, the memory information of the sudden change is lost. The historical effect of this mutation is ignored.
The LSTM network is a cyclic neural network with a special structure design, which is specialized in the ability to remember long-term information, so as to avoid the loss of information in the propagation process. The LSTM network is also composed of multiple groups of repeated neuron modules, and its internal structure is shown in Figure 4.
The most notable feature of the LSTM network is that there is a data route that runs through all neuron modules, which is used to transfer the memory state. The memory route is transmitted in parallel with the input-output route of the neuron, and they interact linearly through a specially designed "gate" structure to ensure the stability of the memory information state. Through the "gate" structure, the memory information can be partially deleted or updated, and stored as long-term memory (cell state) and short-term memory (hidden state) by LSTM network. Taking the LSTM network model structure in Figure 4 as an example, when the input value at a certain moment and the memory ℎ at the previous moment enter the neuron structure, the forgetting vector that determines the degree of memory retention and is generated through the memory gate is: Selecting and uploading those memory states that need to be updated in input gate: in the equation, is the saved vector that determines the update state and is the value required to be updated. Through simple linear calculation of memory gate and input gate, the updated memory state can be obtained: The most notable feature of the LSTM network is that there is a data route that runs through all neuron modules, which is used to transfer the memory state. The memory route is transmitted in parallel with the input-output route of the neuron, and they interact linearly through a specially designed "gate" structure to ensure the stability of the memory information state. Through the "gate" structure, the memory information can be partially deleted or updated, and stored as long-term memory (cell state) and short-term memory (hidden state) by LSTM network. Taking the LSTM network model structure in Figure 4 as an example, when the input value x t at a certain moment and the memory h t−1 at the previous moment enter the neuron structure, the forgetting vector f t that determines the degree of memory retention and is generated through the memory gate is: Metals 2020, 10, 1588 7 of 14 Selecting and uploading those memory states that need to be updated in input gate: in the equation, i t is the saved vector that determines the update state and C t is the value required to be updated. Through simple linear calculation of memory gate and input gate, the updated memory state C t can be obtained: Finally, the output of the neuron structure comes from the linear transformation of C t and [x t , h t−1 ] through the output gate: According to the forward propagation principle of the aforementioned LSTM network, this article built a LSTM network model for predicting the constitutive relationship of Inconel 718. The model consists of an input layer containing three nodes, a LSTM hidden layer containing 100 nodes, two fully connected layers containing 50 nodes, a dropout layer and a regression layer.

Training of the LSTM Constitutive Model
Then, 20 sets of true stress-true strain data obtained by Inconel 718 uniaxial tensile experiment were used for training of LSTM network. Two sets of data under the conditions of 1000 • C, 0.1 s −1 and 1100 • C, 0.0005 s −1 were selected as the validation set, and the remaining data were used as the training set. In the tensile experiments, in order to ensure the same time interval between each data point, the experimental equipment used a fixed acquisition frequency, which lead to a large difference in the length of the experimental data under different strain rates. The data length under the strain rate of 0.001 s −1 at the same temperature is 100 times that of 0.1 s −1 . To facilitate data processing, the experimental data was truncated with 700 as the length threshold, and 18 groups of experimental data were divided into 72 training set samples.
In order to alleviate the over-fitting problem of the LSTM recurrent neural network, the dropout method was used to randomly invalidate 50% of the hidden nodes in each iteration. The invalidated neuron nodes no longer participated in the forward calculation in this iteration and the weight and bias terms were not updated.
The Adam optimization algorithm, that can update the weight of LSTM iteratively based on the training set, and design independent adaptive learning rate for different parameters by calculating the first and second moment estimation of gradient, was adopted as the network training algorithm. The training initial learning rate was set as 0.01, and the maximum number of iterations was 1000 times. Through 603 iterations, the RMSE of LSTM network reached 2.429.

Prediction Results of LSTM Constitutive Model
The rheological properties of Inconel 718 at different temperatures and strain rates were predicted using the LSTM constitutive model after training. The results are shown in Figure 5.
It can be seen that the fitting results of LSTM constitutive model are in good agreement with the measured results. The correlation coefficient (R) of the fitting results and measurement results under different deformation conditions was calculated by formula (13) and listed in Table 2.

Prediction Results of LSTM Constitutive Model
The rheological properties of Inconel 718 at different temperatures and strain rates were predicted using the LSTM constitutive model after training. The results are shown in Figure 5. It can be seen that the fitting results of LSTM constitutive model are in good agreement with the measured results. The correlation coefficient (R) of the fitting results and measurement results under different deformation conditions was calculated by formula (13) and listed in Table 2.    In order to further verify the prediction accuracy of LSTM model, the prediction results at 1000 • C, 0.1 s −1 and 1100 • C, 0.0005 s −1 were compared with the results of M-JC model. The Johnson-Cook model is a typical phenomenological constitutive model, which can effectively characterize the independent effects of temperature, strain rate and strain on flow stress. Taking the coupling effect between deformation conditions into consideration, the modified Johnson-Cook model [26,27] was proposed. It can be expressed as: where σ is the flow stress of the material, ε is the strain, and . ε and . ε re f is strain rate and reference strain rate respectively, T and T re f is temperature and reference temperature respectively, A 1 , B 1 , B 2 , C 1 , λ 1 and λ 2 are material constants. Based on the measured data of Inconel 718 high temperature tensile tests, the unknown material constants of M-JC model can be obtained by linear fitting method. Taking 900 • C as the reference temperature and 0.001 s −1 as the reference strain rate, the M-JC model using the procedure modified by Lin et al. [26] of Inconel 718 can be deduced: Figure 6 shows the prediction curves of LSTM constitutive model and M-JC model under different temperature and strain rates. It can be seen that the LSTM constitutive model has high prediction accuracy in the elastic stage of high temperature deformation. Its RMSE is 0.034 at 950 • C, 0.0005 s −1 and 0.028 at 1050 • C, 0.001 s −1 . In the plastic deformation stage, especially before the necking of specimen, LSTM constitutive model has significant advantages in predicting the flow stress. Table 3 shows the root mean square error and the correlation coefficient of the plasticity stage predicting results of two constitutive models at 950 • C, 0.0005 s −1 and 1050 • C, 0.001 s −1 .

Validation of Prediction Performance of LSTM Constitutive Model under Complex Loading Paths
In order to verify the ability of the trained LSTM constitutive model to predict the flow stress changes under complex loading paths, the high temperature tensile experiment including loading and unloading stages under different strain rates were designed. The process parameter curve is shown in Figure 7a. The drawing experiment fully refers to the hot forming process curve of Inconel 718 in the actual project, and the comparison between the measured result and the predicted curve of LSTM rigid model is shown in Figure 8.
It can be seen that compared to the M-JC model and other phenomenological constitutive models, the LSTM constitutive model still has high prediction accuracy when the load, unloading and strain rates change suddenly. The RMSE of the prediction results are 2.763. In addition, the LSTM constitutive model has significant advantages in predicting the elastic deformation stage of Inconel 718, the RMSE is 0.075 and 0.064. It can also be found from the results that the prediction model has large errors in the prediction of Inconel 718's creep phase and stress relaxation phase. This is because the training data was obtained under the condition of constant rate stretching, and there is less data in the viscoelastic phase. By improving the cell structure of the LSTM model and incorporating the standard creep test data into the training set, the accuracy of the model's prediction of the Inconel 718 viscoelastic phase rheological properties can be significantly improved.
prediction accuracy in the elastic stage of high temperature deformation. Its RMSE is 0.034 at 950 °C, 0.0005 s −1 and 0.028 at 1050 °C, 0.001 s −1 . In the plastic deformation stage, especially before the necking of specimen, LSTM constitutive model has significant advantages in predicting the flow stress. Table 3 shows the root mean square error and the correlation coefficient of the plasticity stage predicting results of two constitutive models at 950 °C, 0.0005 s −1 and 1050 °C, 0.001 s −1 .

Validation of Prediction Performance of LSTM Constitutive Model under Complex Loading Paths
In order to verify the ability of the trained LSTM constitutive model to predict the flow stress changes under complex loading paths, the high temperature tensile experiment including loading and unloading stages under different strain rates were designed. The process parameter curve is shown in Figure 7a. The drawing experiment fully refers to the hot forming process curve of Inconel 718 in the actual project, and the comparison between the measured result and the predicted curve of LSTM rigid model is shown in Figure 8.    It can be seen that compared to the M-JC model and other phenomenological constitutive models, the LSTM constitutive model still has high prediction accuracy when the load, unloading and strain rates change suddenly. The RMSE of the prediction results are 2.763. In addition, the LSTM constitutive model has significant advantages in predicting the elastic deformation stage of Inconel 718, the RMSE is 0.075 and 0.064. It can also be found from the results that the prediction model has large errors in the prediction of Inconel 718's creep phase and stress relaxation phase. This is because the training data was obtained under the condition of constant rate stretching, and there is less data in the viscoelastic phase. By improving the cell structure of the LSTM model and incorporating the standard creep test data into the training set, the accuracy of the model's

Conclusions
The following conclusions can be derived from this study: (1) Uniaxial tensile experiments were conducted on Inconel 718 at different temperatures and strain rates. Under different deformation conditions, the material exhibits obvious work hardening and cyclic softening coupling characteristics. Through the analysis of the true stress-true strain curve, it is found that the rheological behavior of Inconel 718 at high temperature has significant historical correlation with loading. (2) A long short term memory recurrent neural network was trained using experimental data to characterize the constitutive relation considering the loading history effect. The model has a multi-hidden layer structure composed of LSTM layer and full connection layer. The Adam optimization method was adopted as the training algorithm, and the RMSE of training set after 603 iterations was 2.429. Funding: This research received no external funding.

Conflicts of Interest:
The authors declare no conflict of interest.