Next Article in Journal
Vertical and Lateral Dynamics of 4L Freight Bogie
Previous Article in Journal
Dynamics of Interacting Colloidal Particles Using the IIR Recursive Digital Filter Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Model of Structural Phase Transitions in Al-Cu Solid Solutions under Dynamic Loading Using Machine Learning

Department of General and Theoretical Physics, Chelyabinsk State University, Chelyabinsk 454001, Russia
*
Author to whom correspondence should be addressed.
Dynamics 2024, 4(3), 526-553; https://doi.org/10.3390/dynamics4030028
Submission received: 21 April 2024 / Revised: 5 July 2024 / Accepted: 9 July 2024 / Published: 12 July 2024

Abstract

:
The development of dynamic plasticity models with accounting of interplay between several plasticity mechanisms is an urgent problem for the theoretical description of the complex dynamic loading of materials. Here, we consider dynamic plastic relaxation by means of the combined action of dislocations and phase transitions using Al-Cu solid solutions as the model materials and uniaxial compression as the model loading. We propose a simple and robust theoretical model combining molecular dynamics (MD) data, theoretical framework and machine learning (ML) methods. MD simulations of uniaxial compression of Al, Cu and Al-Cu solid solutions reveal a relaxation of shear stresses due to a combination of dislocation plasticity and phase transformations with a complete suppression of the dislocation activity for Cu concentrations in the range of 30–80%. In particular, pure Al reveals an almost complete phase transition from the FCC (face-centered cubic) to the BCC (body-centered cubic) structure at a pressure of about 36 GPa, while pure copper does not reveal it at least till 110 GPa. A theoretical model of stress relaxation is developed, taking into account the dislocation activity and phase transformations, and is applied for the description of the MD results of an Al-Cu solid solution. Arrhenius-type equations are employed to describe the rates of phase transformation. The Bayesian method is applied to identify the model parameters with fitting to MD results as the reference data. Two forward-propagation artificial neural networks (ANNs) trained by MD data for uniaxial compression and tension are used to approximate the single-valued functions being parts of constitutive relation, such as the equation of state (EOS), elastic (shear and bulk) moduli and the nucleation strain distance function describing dislocation nucleation. The developed theoretical model with machine learning can be further used for the simulation of a shock-wave structure in metastable Al-Cu solid solutions, and the developed method can be applied to other metallic systems, including high-entropy alloys.

1. Introduction

The theoretical description of the mechanical behavior of materials is still an urgent problem, especially under dynamic loading conditions, despite the great research activity over several decades. The gradual investigation of the dominant mechanisms leading to the evolution of material structure during loading has led to the creation of numerous models predicting material properties. Under dynamic loading, the equivalent shear stress exceeds the quasi-static yield stress due to a finite rate of defect microstructure evolution, and the purpose of a dynamic plasticity model is to take this finite rate of relaxation into account. The simplest way to describe the dynamic plasticity is to introduce a dynamic yield stress as a material parameter [1,2], which is, unfortunately, valid only for specific test conditions. Semi-empirical constitutive models [3,4,5] were created to determine the strain rate and temperature dependence of material yield behavior from the consideration of defect microstructure evolution, typically the thermally activated dislocation motion with parameters fitting to the experiments [6]. In general, the dislocation-mediated plasticity is the main mechanism for most crystalline metals, which is supported by both molecular dynamics (MD) simulations [7] and experimental studies [8]. The incorporation of a crystal plasticity approach with an explicit description of the dislocation activity in different slip planes is a fruitful approach for quasi-static [9] and dynamic [10,11,12] applications. This approach is successfully used to describe the shock-wave processes in solid metals [13,14,15].
In addition to the dislocation plasticity, the plastic relaxation can be caused by several other mechanisms, which can act simultaneously with the dislocation activity and with each other. One can mention the interfacial slip and diffusion processes [16], mechanical twinning [17,18,19,20] in the case of coarse-grained solids, as well as additional grain boundary migration [21,22,23] and grain boundary sliding with grain rotation [24,25,26] in the case of fine-grained solids. In addition, a structural phase transition can have a significant effect on the dynamic response of metals. The phase transition of FCC (face-centered cubic) to BCC (body-centered cubic) leads to excessive stress relaxation according to MD simulations [27] at high compressive strain rates. Shock-induced phase transitions are revealed by means of MD simulations in single-crystal aluminum [28,29], copper [30], iron [31], cerium [32], magnesium [33] and high-entropy alloys [34,35]. The FCC-to-BCC phase transition in copper was experimentally revealed in the quasi-static measurements by Dewaele et al. [36] with the diamond anvil cell, as well as in dynamic measurements with shock-wave compression [37] and ramp compression [38]. The large number of experimentally detected phase transitions in materials has led to the development of mechanical models that take into account phase transitions. Levitas and Javanbakht [39,40] proposed a phase field model with in-deep thermodynamics consideration for the description of martensitic transformations. Yeddu and Lookman [41] developed a variant of the phase field model based on the time-dependent Ginzburg–Landau kinetics equation for the assessment of the β(BCC)-to-ω(hexagonal) phase transformation in Zr–Nb alloys. Consideration of multiple possible variants of austenite-to-martensite transformation in a quenching and partitioning steel QP1180 allowed Yang et al. [42] to describe correctly the strain path and orientation dependencies of the transformation process. In spite of these efforts, there is still a lack of a simple and robust model of plastic relaxation with the interplay of dislocation activity and phase transitions. Combining MD data, theoretical framework and machine learning (ML) methods is a prospective approach in the field.
The equation of state (EOS) with the additional consideration of the deviatoric part of the stress tensor becomes the core part of many constitutive models for shock-wave problems. At a very high impact level, the deviatoric part becomes negligibly low [43] and the EOS completely describes the behavior of matter. In general, the construction of the EOS consists of proposing a theoretically based functional form and fitting the relevant parameters to the available experimental data and first-principles calculations [44,45]. This appears to be costly for complex substances and multi-component alloys. The application of ML methods, in particular artificial neural networks (ANNs), is a perspective method of solving this problem [46].
ML is increasingly being used in material science and in the mechanics of materials [47,48]. ML methods are an efficient way to account for changes in mechanical properties of materials and are quite difficult to implement by classical methods and require an analysis of a large amount of experimental data and simulation results. The problem statement such as prediction, recognition and classification determines the choice of the most effective ML algorithm such as linear and logical regression, decision trees and clustering [49,50,51]. Among various ML methods, ANNs and Bayesian global optimization of model parameters are of particular interest for the present study. ANNs can be used to approximate the strain, strain rate and temperature effect on the yield surface more precisely in comparison with the known analytical dependencies [52,53], to establish structure–property relationships [48,54,55] and to construct a surrogate model of the material [56,57]. Bayesian calibration of model parameters was fruitfully used to calibrate the empirical Johnson–Cook model using the data on plate impact [58] and Taylor tests [59], as well as to fit a dislocation plasticity model to the Taylor test data [60,61,62].
In this work, we simulate the behavior of pure metals, aluminum and copper, as well as Al-Cu solid solutions with different contents of elements under uniform high-strain-rate compression. The mechanical behavior of the material is described by a coupled model of phase transition and plasticity with the Arrhenius-type equations for the rates of phase transitions and a modified Maxwell-type relaxation model [63] with a variable relaxation time that takes into account the evolution of the dislocation structure. The parameters of the stress relaxation and phase fraction evolution model are fitted by the Bayesian method to the results of the MD simulations. In addition to the coupled model of phase transition and plasticity, we have developed a deep forward-propagation ANN to describe the EOS of alloys. The ANN relates pressure as the hydrostatic part of stresses, bulk modulus and temperature with density and internal energy simultaneously for different concentrations of copper in the Al-Cu solid solution ranging from 0% Cu to 100% Cu. Another forward-propagation ANN allows the calculation of the shear modulus and dislocation nucleation threshold for the same continuous range of Cu concentration.

2. MD Study of Uniaxial Compression of Al-Cu Solid Solutions

2.1. MD Setup

We performed MD simulations of uniaxial dynamic compressive and tensile deformations of single crystals of Al, Cu and Al-Cu solid solutions with different component contents with an engineering strain rate of 109 s−1, which is close to that achieved in state-of-art experiments on irradiation of thin metal foils with powerful sub-picosecond laser pulses [64,65,66,67]; thus, this strain rate is experimentally attainable. The single crystal is oriented so that the [100], [010] and [001] crystallographic directions are along the x-, y- and z-axes (Figure 1). The MD system contains half a million atoms and is originally cubic with 50 × 50 × 50 lattice periods in size. Periodic boundary conditions are set in each of the three directions.
It is known that increasing the concentration of the alloying element affects the mechanism of dislocation motion and leads to higher stresses [68,69,70,71]. Usually, the proportion of the additive element in alloys is small—less than 5% Cu in Al matrix and up to 18% of Al in Cu matrix, but the alloying element can locally reach high concentrations in material. It is shown in [72] that the segregation of solid solution copper atoms reaches 15.3 wt % at grain boundaries during cryogenic surface mechanical grinding treatment. The formation of solid solution precipitates with high content of Cu in Al alloy during cyclic hardening is also possible [73]: the Cu content can reach almost 50 wt % according to the electronic microscopy images in this paper. Components of high concentrations having equal atomic fractions form high-entropy alloys with unique mechanical properties. These facts motivated us to consider copper concentrations of 10, 20, 30, 50, 70 and 80% in order to investigate theoretically the deformation behavior of the solid solution, although a part of the considered solid solutions can be thermodynamically unstable in traditional metallurgical casting. On the other hand, if they reveal unique dynamic-protecting properties, novel production techniques, such as additive manufacturing, can be applied to produce such metastable alloys; therefore, theoretical analysis of their properties is relevant. Aluminum atoms in the system were replaced by copper atoms in random order to form the solid solution. The compression and tension of the samples were modeled at constant temperature in the range of 100 to 900 K in steps of 100 K. Thus, 72 MD simulations were performed in total. The calculations were performed using the LAMMPS software package [74] with an interatomic potential of the angle-dependent potential (ADP) type [75], which is one of the most reliable classical force fields. The Nose–Hoover thermostat [76] and barostat were used to relax the stresses for 10 ps before the deformation stage. The crystal was deformed along the x-axis with fixed transverse dimensions. The deformation lasted for 300 ps, achieving an engineering strain value of 0.3. For comparison, auxiliary modeling was also carried out in the adiabatic approximation, that is, without a thermostat. Also, several MD systems were subjected to compression with other strain rates of 108 s−1 and 3·108 s−1 to reveal the rate effect.

2.2. Results of MD Simulations: Dislocation Activity and Phase Transitions

Data from MD simulations of the uniaxial compression of crystals for a wide temperature range were analyzed using the “Polyhedral template matching” algorithm [77] integrated into the OVITO software [78]. This method classifies the local structural environment of the particles. The phase fraction was calculated as the ratio of the number of atoms of a certain phase to the total number of atoms of the modeled crystal. Figure 2 shows the evolution of shear stress τ = σ x x σ y y / 2 and atomic structure for pure aluminum and copper at 300 K. Both metals are characterized by an elastic response at the beginning of deformation. The appearance of the metastable BCC phase at an engineering strain of about 0.11 leads to an elastically nonlinear behavior. The subsequent appearance of the HCP (hexagonal close-packed) phase reflects the nucleation and movement of partial dislocations (Figure 3a) along the close-packing planes and leads to a sharp drop in stresses. Thus, the HCP phase here is essentially stacking faults left as a result of the slip of partial dislocations. The shear stress level increases in copper during the following plastic flow due to strain hardening. In the case of aluminum, strain hardening also takes place, but the formation of the BCC phase is observed almost throughout the entire volume of the crystal, which leads to rapid stress relaxation; shear stress even drops down to a negative value. The pressure curve reflects the effect of polymorphic phase transition in the case of pure aluminum (Figure 4): The pressure growth slows down above 35 GPa, which is associated with a change in the structure of the loaded material. The results of the MD data give a slightly higher pressure value compared to the experimental data [36] on compression in diamond anvil cells and DFT calculations [79]. In general, all considered data are consistent.
Let us consider in detail the behavior of an Al-Cu solid solution by an example of the crystal with 10% Cu in aluminum at the temperature of 300 K. Figure 5 shows the atomic structure of the crystal and the evolution of phase fractions and stresses. At the beginning of the simulation, the crystal has a primarily FCC crystal structure and deforms elastically. The shear stress reaches higher values than in pure aluminum due to the strengthening by copper atoms. Then, dislocation loops are formed in the material (Figure 3b). At this stage, the fraction of the metastable BCC phase is significantly less than in pure aluminum and copper. Active multiplication of Shockley partial dislocations takes place. The shear stress relaxes down to 2 GPa due to the plastic flow at an engineering strain of 0.11. In the course of further deformation, there is a gradual increase in the fraction of the BCC phase, while the FCC and HCP phases completely disappear. The transformation of the crystal structure entails a decrease in shear stresses to negative values. After that, the shear stress increases again due to the following deformation of the stable BCC phase. For the cases of solid solutions with higher copper concentrations (Figure 6), the process of dislocation plasticity plays an insignificant role in the stress relaxation in the material. Stress relaxation occurs primarily due to the transition of the initial structure into the BCC phase and the regions of disordered atoms (“OTHER” or amorphous phase). Phase transitions completely suppress dislocation activity. For most of the modeled crystals, the process of BCC phase formation from disordered atoms is common for the considered range of engineering strain. When the Cu fraction increases up to 50%, the amplitude of the maximum shear stress decreases and the beginning of the phase transition occurs at lower strains. The distribution of phase fractions for a solid solution with 80% Cu is similar to the case with 20% Cu. Figure 3 shows an increasing trend in the dislocation density with deformation change by a decreasing trend for pure Al and all investigated Al-Cu solid solutions; this is because the phase transition decreases the volume of the FCC phase and, thus, the total length of the dislocations. The maximal dislocation density increases with the addition of copper atoms up to the Cu concentration of 20%, but decreases with a further increase in Cu concentration, see Figure 3b.
Let us consider the effect of temperature on the stress–strain state of the material (Figure 7a). The material softens, and the amplitude of the shear stresses decreases with increasing temperature. The MD data analysis shows an increase in the fraction of disordered atoms and a decrease in the phase transition threshold with increasing temperature. To determine the influence of the strain rate on the phase transition process, additional modeling of solid solution compression with strain rates of 108 and 3·108 s−1 at a temperature 300 K was carried out (Figure 7b). The change in strain rate by an order of magnitude does not significantly affect the evolution of the phase composition. On the contrary, the plasticity process is sensitive to changes in the strain rate, which is evident from the strong variation in the content of the HCP phase presenting the traces of dislocation gliding, Figure 7b. Consequently, the phase transition is strain-dependent instead of strain-rate sensitive in the studied dynamic loading range. It should be noted that there is a transition from the BCC phase to the region of disordered atoms (“OTHER” phase) when the strain is greater than 0.3. The crystal is almost completely disordered at the engineering strain of 0.45.
As we show further, it is possible to formulate a theoretical model of coupled phase transition and plasticity for alloys on the basis of the obtained MD data on phase transformation under dynamic compressive deformation. Such a model must take into account the formation of the BCC lattice from the FCC and “OTHER” structures, as well as the increase in the fraction of the “OTHER” structure due to the decrease in the fractions of the FCC and BCC lattices.

3. The Constitutive Equations in the Form of the Artificial Neural Network

Forward-propagation artificial neural networks (ANNs) are used to approximate the constitutive equations of the Al-Cu alloy with a solid solution of copper atoms with different concentrations of components. The stress–strain relation for solids in the elastic domain can be presented in the form of the generalized Hook’s law [83]:
σ i j = C i j p q   ε p q ,
in which the components of the tensor C i j p q of elastic moduli are not constants and are introduced as functions of the current strain tensor C i j p q = C i j p q   ( ε ) . We considered that the isotropic formulation and the tensor C i j p q   ( ε ) has only two independent components, i.e., Lamme coefficients. These coefficients can be resolved through other quantities that can be calculated from numerical calculations and experiments: the shear modulus G and the bulk modulus K, which will be considered in the remaining part of the paper. The constitutive equations in this work are based on uniaxial compression data. In this case, the deformation was determined by the ratio of the deformed and initial volumes ε = V / V 0 1 or equivalent ε = ρ 0 / ρ 1 . Therefore, we redefined the strain dependence on ε in Equation (1) to the density ρ dependence, which will greatly simplify the use of these constitutive equations in the macroscopic simulations in both mesh- and particle-based methods. The constitutive equations in the ANN form can be used in any other numerical schemes by implementing a forward propagation on the ANN parameters (weights and biases). The first ANN, whose scheme is presented in Figure 8, is essentially a scalar EOS and describes the following vector function:
y = f ( α C U , T , ρ ) ,
where the output layer of the ANN is y = P , E , K . P is the pressure; E is the specific internal energy of material, and K is the bulk modulus. The input layer of the ANN consists of the arguments of the vector function (2). T is the temperature; ρ is the density, and α C U 0 , 1 is the atomic fraction of Cu atoms in the solution. The training and test data sets were constructed from MD simulation data of uniaxial compression and tension of the representative volume of the Al-Cu alloy with the solid solution of copper atoms. In the case of tension, the deformation trajectories of MD data were restricted by the beginning of the fracture. The values of density, pressure, total energy and temperature were extracted directly from the MD simulation data, and the bulk modulus was estimated from the pressure curve as follows:
K i = j = i n j = i + n P j × ε j j = i n j = i + n ε j 2 1 ,
where i is the state point for the training and test data set, n is the averaging interval, Pj is the uniaxial tension/compression pressure and ε j is the uniaxial tension/compression strain value; the step of strain in MD output data is 10−4.
The second ANN (Figure 9) was used to describe the second part of the constitutive equation of the Al-Cu alloy with a solid solution of copper atoms, namely, the value of the shear modulus and dislocation nucleation threshold during uniaxial deformation. The function arguments were also determined by Equation (2), but the function values were set as y = ( G , Q ε ) , where G is the shear modulus and Q ε is the nucleation strain distance function introduced in [65]. The shear modulus at state point i was calculated as
G i = j = i n j = i + n S x x   j × ε j j = i n j = i + n ε j 2 1 ,
where S x x   j is the component of the stress tensor deviator and i, n and j are indices, which are defined in the same way as for Equation (3). The function Q ε characterizing the nucleation threshold was constructed as a continuous one, since an ANN is capable of approximating only continuous dependencies with high precision:
Q ε = ε j ε n u c ,
where ε j is the current strain points of the system under uniaxial tension/compression and ε n u c is the threshold value for the onset of dislocation nucleation under uniaxial deformation. The nucleation of dislocations under uniaxial compression at the fixed strain rate begins when the condition Q ε 0 is satisfied. The inclination of Q ε in the strain space allows one to take into account the strain-rate dependence of the nucleation threshold [84].
All results obtained by MD simulations were divided into training and test data sets. The test data were used only to check the accuracy and generalizability of the ANN by the cross-validation method and did not take part in the training. The data sets contained a total of 285,715 data vectors ( α C U , ρ , T, P, E and K) for the first ANN and 133,763 data vectors ( α C U , ρ , T, G and Q ε ) for the second ANN. The data distribution and the data ranges are given in Table 1. All the data underwent a normalization procedure using the standard definition of MinMax normalization:
x i j n = x i j x min j x max j x min j , i = 1 , , n r o w s , j = 1 , , n c o l s ,
where x i j n is the normalized value, index j defines the data column, e.g., P, E, K, α C U , T and ρ for the first ANN, index i defines the data row, x max j is the maximum value in the j-th column and x min j is the minimum value in the j-th column. The range of values was expanded by 20% for each column relative to its real minimum and maximum. That is, the values x max j and x min j were shifted by 10% relative to the real minimum and maximum values in the original data sets. This is because the derivative of sigmoid transfer function is small at the extreme values close to 0 and 1 and the learning process becomes slower. Another motivation is to prevent going beyond the limits.
The ANNs were constructed as classical fully connected neural networks [85] in which the values obtained through the activation function of neurons were transmitted to the next layers of the network (Figure 8 and Figure 9 show ANNs in the form of computational graphs, where the arrows indicate the direction and sequence of the calculations). At the current layer l of the ANN, the activation function a j l depends on the input z j l to the neuron:
z j l = k w j k l a k l 1 b j l , l = 1 , 2 , , n ,
where n is the number of layers in the ANN, j is the index by neurons in the l-th layer, k is the index by neurons in the l − 1 layer, w j k l are the weight coefficients for the j-th neuron, a k l 1 is the neuron function value from the previous layer and b j l are biases values for the j-th neuron. The ANN contains several types of layers; the first one is the input layer for which l = 0, a k 0 = x 1 , x 2 , , x m , where x 1 , x 2 , , x m are the arguments of the approximating function. The arguments are α C U , T and ρ for both ANNs. The second type are the hidden layers, where the main transformation of the input data takes place, for them l = 1 , , n 1 . The last type is the output layer, where l = n = L . The output layer gives the values of the function P, E and K for the first ANN and G and Q ε for the second ANN.
The first ANN (Figure 8) has five hidden layers of 20 neurons each and the output layer that contains 3 neurons. The artificial neurons in hidden layers are described by the PReLU activation function [86]:
a j l = z j l i f   z j l > 0 a p r   z j l i f   z j l 0 ,
where a p r is the parameter of the PReLU function. The sigmoid activation function is used to describe the output layer in this ANN:
a j L = 1 1 + e z j L ,
Such a combination of activation functions makes it possible to describe highly non-linear dependencies, despite the fact that the PReLU function in the positive and negative regions is separately linear [85]. The second ANN (Figure 9) contains six hidden layers of 20 neurons each and the output layer that contains 2 neurons. Swish activation function is used in hidden layers [87]:
a j l = z j l 1 + e β z j l ,
where β is the parameter of Swish function. The sigmoid function in Equation (9) is also used as the activation function on the output layer. These ANN architectures were chosen because they show the least deviation from the MD simulation results.
The next step was to train the ANNs after building the training and test data sets. The stochastic gradient descent (SGD) method was used to train the neural networks: this method minimizes the cost function that defines the accuracy of the ANN. The cross-entropy cost function was chosen for both ANNs, which together with the output sigmoid activation function (9), gave the following expression:
L = 1 N i = 1 N y i ln y i ^ + 1 y i ln 1 y i ^ ,
where N is the number of outputs of the ANN, y i are the real values of the function at the point and y i ^ are the ANN results at the point. The cross-entropy is typically used for classification problems, but our previous studies showed its efficiency for regression problems as well [62,88]. SGD is classically implemented in fully connected forward-propagation ANNs based on a backpropagation algorithm. In this work we also used the optimization of the SGD algorithm called Adam [89]. The variable learning rate is implemented in the Adam method (i.e., gradient descent steps change with training progress), which allows one to quickly “go down” to the minimum of the cost function (11). Also, the Adam method provides greater accuracy of the ANN results and, at the same time, avoids the instability of the ANN training process (exploding or vanishing gradient). Table 1 and Table 2 present the main characteristics and values of the hyperparameters for the first and second ANNs. The formulas for Adam’s algorithm are given in Appendix A.
The ANN allowed us to predict data for continuous values of input parameters. Figure 10 shows the pressure dependencies for the intermediate values of copper concentrations at the temperature of 300 K. These data were not included in the test and training data sets. The ANN qualitatively describes both an increase in pressure during the compression of a substance and a drop in pressure during the stretching process. There is an increase in the density of the substance and a deceleration in pressure growth with an increase in the concentration of copper.
Next, we compared the temperature dependencies of the bulk modulus B and shear modulus G (Figure 11 and Figure 12) predicted by the ANNs with data from other studies in the literature. The main experimental results are given for the single-crystal aluminum [90,91,92,93] and single-crystal copper [94,95,96]. A comparison of the temperature dependence of the shear modulus (Figure 12b) also contained the theoretical estimate [97] and experimental results for the polycrystal [98]. The bulk modulus of copper decreases with increasing temperature by about 20 GPa at temperatures close to 900 K, both according to our results and experimental data as can be seen in Figure 11a. There is some scatter in the experimental data of the bulk modulus of pure copper from different sources, but the ANN results generally agree with the considered data and lie between the experimental points (Figure 11a). The bulk modulus of pure aluminum (Figure 11b) decreases by 5 to 30 GPa depending on the considered studies with the increasing temperature up to 900. The results of our work show that the bulk modulus of pure aluminum drops by about 10 GPa, and these results also lie between the experimental points.
The comparison of the obtained values of the shear modulus from the trained ANN with the experimental data for pure copper (Figure 12a) shows the underestimation of the modulus value at the initial temperature of 100 K; the difference is about 1.5 GPa. In this case, the temperature dependence of the shear modulus predicted by the ANN shows the same trend as the experimental data: There is a decrease by approximately 8 GPa at a temperature of 800 K. In the case of pure aluminum (Figure 12b), the shear modulus predicted by the ANN describes the experimental data somewhat. The deviation manifests itself at high temperatures, where experimental data show a larger drop in the values of the shear modulus. Nevertheless, the trend toward a decrease in the shear modulus is observed both in our results and in the data from other works presented in Figure 12b. The deviation of the ANN prediction may be due to the used interatomic potential. In the recent work [99], the shear modulus is calculated for various temperatures by the MD method using EAM interatomic potential. The shear modulus values calculated for <110> crystallographic orientation and at the temperature 100 K sharply drop to approximately 15 GPa and decrease with a further increase in temperature, which does not correspond to any of the data presented in Figure 12b. This stimulates the further development of interatomic interaction potentials, such as NNP (neural network potential) [100] or PINN (physics-informed neural network potential) [101]. These types of interatomic potentials describe physical processes with sufficient accuracy, but at this stage of development require a lot of computational time in MD simulations. The accuracy of the constitutive relations in the form of ANNs trained on data obtained using such potentials will also increase.
The ANNs were tested by comparing the density ρ , bulk B and shear G modulus in the Al-Cu alloy with a copper solid solution, where the copper concentration of up to 6 at.% was considered (Figure 13 and Figure 14). The results of the ANN (this work), experimental data [102] and ab initio calculations [103] were compared. A comparison of the density of the alloy (Figure 13) shows that, both in the ANN predictions and in the experiment, a linear increase in density with an increase in the concentration of copper atoms is observed. Although the experiment shows the less sharp slope of the straight line, in general, all the considered data in Figure 13 are in agreement. The bulk modulus predicted by the ANN for the temperature of 300 K (Figure 14a) is somewhat overestimated in comparison with the experimental data and ab initio calculations. In this case, the experimental data [102] are presented for the polycrystal, which may give some discrepancy. The interatomic potential used can also lead to inaccuracies in the training and test data sets built on the basis of MD simulations. The shear modulus (Figure 14b) obtained from the ANN for the temperature of 300 K shows some average behavior between the experimental results and the ab initio calculation. Nevertheless, the shear modulus for the Al-Cu alloy grows faster with an increasing copper concentration according to the ANN predictions. This may be due to the fact that the density of this material predicted by the ANN increases somewhat faster compared to the experiment (Figure 13), which can explain the faster growth of both moduli.

4. Shear Stress Relaxations Due to Dislocation Activity and Phase Transitions

4.1. Theoretical Model of Dislocation Plasticity and Phase Transitions

The results of the MD simulations show the relaxation of shear stresses in the process of uniaxial compression of single crystals of the Al-Cu solid solution. The presence of a stacking fault and the formation of dislocations indicate the plastic response of the material mediated by dislocation activity, while the evolution of phase fractions accompanies phase transitions during deformation. The propagation of a plane shock wave caused, for example, by detonation or a high-velocity impact, increases pressure up to tens of gigapascals in the condition of uniaxial deformation. Such severe deformation should be analyzed in terms of finite strains. Following the idea of multiplicative decomposition [40,104,105,106], the macroscopic strain gradient F is represented as follows:
F = F e F t r F p ,
where F e is the elastic part of the deformation to be defined for the stress calculation, while F p is the contribution of plastic (irreversible) deformation due to the dislocation activity and F t r is the contribution of phase transitions. In the considered case of uniaxial deformation, the macroscopic strain gradient can be calculated as follows:
F = L x / L 0 0 0 0 1 0 0 0 1 ,
where F 11 is the ratio of the current length L x along the deformation axis to the initial length L 0 .
The atom rearrangement during phase transition is equivalent in a certain sense to an imposed deformation. Therefore, the contribution of phase transition can be represented as a linear combination of the atomic fractions of the emerging phases ( ω B C C for BCC and ω O T H E R for disordered atoms—“OTHER” phase):
F 11 t r = n 1 ω B C C + n 2 ω O T H E R + 1 , F 22 t r = F 33 t r = 1 F 11 t r .
Here, we assumed F 11 t r F 22 t r F 33 t r = 1 and did not consider the volume change at the phase transitions, because this volume change influencing the pressure was already taken into account in the ANN-based EOS. In Equation (14), n 1 and n 2 are the model parameters to be fitted to the MD data.
The next step is the equations of balance among the atomic fractions of the FCC, BCC and OTHER (disordered or amorphous) structures. The fraction of the FCC structure decreased due to the formation of the BCC and OTHER structures. An analysis of the MD data showed that the mutual transformation of the BCC and OTHER structures is observed during compression. Consequently, the rate of change in phase fractions in the material during deformation can be written as follows:
d ω F C C d ε = K 1 ω F C C K 2 ω F C C , d ω B C C d ε = K 1 ω F C C + K 4 ω O T H E R K 3 ω B C C , d ω O T H E R d ε = K 2 ω F C C + K 3 ω B C C K 4 ω O T H E R .
Here, ω F C C represents the sum of atomic fractions of both the FCC and HCP phases, because the HCP phase arises as a result of dislocation activity inside the FCC phase, but not due to a phase transition in the classical sense. The parameters in the expression for the coefficients in linear combinations depend on the formation of a certain structure in the phase transition. Based on the Arrhenius equation, it can be written as follows:
K β = K 0 β exp U β + V 0 ε 1 β σ 11 + 2 ε 2 β σ 22 k B T , β = 1 ,   at   the   transition   F C C B C C 2 ,   at   the   transition   F C C O T H E R 3 ,   at   the   transition   B C C O T H E R 4 ,   at   the   transition   O T H E R B C C ,
where K 0 β is the phase formation rate parameter, U β is the threshold energy, V 0 = 10 30 m3 is the activation volume, σ 11 and σ 22 are the components of the Cauchy stress tensor and ε 1 β and ε 2 β are the characteristic strains of phase transition along the loading direction and transverse direction, respectively, such that V 0 ε 1 β σ 11 + 2 ε 2 β σ 22 is a positive work of phase transition. The parameters of Equation (16) are also fitted to the MD data.
The components of the plastic strain gradient tensor are defined through the total plastic strain w :
F 11 p = exp w ,   F 22 p = F 33 p = 1 F 11 p ,
where F 11 p F 22 p F 33 p = 1 holds because of the volume conservation during the dislocation slip. Dislocation-mediated plastic deformation occurs, when shear stress exceeds the threshold yield stress Y b / 2 , but the inertness of the dislocation system holds the shear stress level above this threshold under dynamic loading. This inertness can be taken into account through the relaxation time approach [107]. Due to the evolution of the dislocation structure, the characteristic time of stress relaxation changes with the dislocation density [108], such that the following modified Maxwell model can be implemented [63,109]:
w ˙ = χ 1 1 2 S 1 3 Y b sign S H 1 2 S 1 3 Y b ω F C C , χ = 8 B 3 ρ D b 2 ,
where S is the stress deviator, H is the Heaviside step function, ρ D is the scalar density of dislocations and b is the modulus of Burgers vector. According to the results of the MD simulations [110], the dynamic drag coefficient of a dislocation B is directly proportional to temperature; however, following [111], a nonlinear growth of the friction coefficient is observed for copper and aluminum when approaching melting temperatures.
Now we considered the evolution equation for the dislocation density. The first process to be considered was the dislocation nucleation. The formation of a dislocation loop of radius R and length 2 π R in each elementary volume R 3 of a metal is a random process occurring in a single crystal with frequency c t / R due to thermal fluctuations, and its rate per unit volume can be represented as follows [79,84]:
Q n = 2 π R R 3 c t R exp W k B T H Q ε ,
where Q ε is the nucleation strain distance function [84], which shows the threshold of dislocation nucleation as discussed in Section 3, c t = G / ρ is the transverse speed of sound, G is the shear modulus, ρ is the material density and k B is the Boltzmann constant. Following [62], the nucleation energy barrier is presented as W = k n G b 3 , where k n is a parameter. Expressing the concentration of atoms through the density ρ of the solid solution and the average mass of one atom m 1 , we wrote down the dislocation nucleation rate as follows:
Q n = 2 π ρ c t m 1 exp k n G b 3 k B T H Q ε .
The nucleation threshold and shear modulus in our work were calculated using an ANN (see Section 3).
A certain part of plastically dissipated power was spent on the formation of new dislocations. This idea was used in [112] to formulate the dislocation multiplication rate, and we applied a similar approach:
Q m = b ε D 3 2 S w ˙ ,
where ε D is an efficient energy of the dislocation multiplication per length b of the new dislocation line. The annihilation rate of dislocations was written down in the standard form:
Q a = 2 K a ρ D w ˙ ,
where K a is the annihilation coefficient. As a result, the expression of dislocation kinetics took the following form:
ρ ˙ D = Q n + Q m Q a .
The static yield strength was calculated using Taylor’s law of hardening:
Y b = Y 0 + A G b ρ D ,
where A is the hardening parameter. The stresses can be represented through the spherical and deviatoric stress components
σ = P + S .
Usually, the pressure is calculated from the equation of state. In our work, we used an ANN (see Section 3). The input parameters of the ANN were the atomic concentration of copper in solid solution, temperature and current density of the substance:
ρ = ρ 0 2 E 11 e + 1 2 E 22 e + 1 2 .
From Equation (12), we expressed the elastic part of the strain gradient and found the Cauchy–Green strain tensor:
E e = 1 2 F e T F e I .
Thus, we obtained expressions for the stress deviator and shear stresses using Hook’s law [113]:
S = 4 3 G E 11 e E 22 e , τ = 3 4 S .
In the formulated model of shear stress relaxation, it was necessary to determine a set of parameters: K 0 β , U β , V 0 , ε 1 β and ε 2 β in the model of phase structure evolution (Equation (16)), as well as n 1 and n 2 in the model of phase transition-induced stress relaxation (Equation (14)) and B , Y 0 , A , k n , ε D and K a in the model of dislocation plasticity.

4.2. Optimization and Results of Phase Transition Model

The parameter identification method based on the Bayesian approach is a powerful tool of machine learning [85,114]. It involves calculating the probability for many sets of randomly chosen parameter values. The probability characterizes the degree of agreement of the calculated values with the reference data. The probability is taken to be equal to one at the beginning of the modeling. In the process of modeling, the difference between the compared values is accumulated and the probability decreases correspondingly. The parameters are selected in such a way that the resulting probability is maximal for the best parameters of the model. In this case, the modeled curve of the compared quantity will be as close as possible to the reference curve. For the phase composition evolution model, the probability characterizes the degree of correspondence of the BCC and OTHER (disordered or amorphous) phase fractions calculated from the system of Equations (15) and (16) with MD data. The probability was calculated as the product of the probabilities at each modeling step i :
Probability = i exp k j = 1 2 ω i j model ω i j MD Δ ω j MD 2 .
where ω i 1 and ω i 2 are the atomic fractions of the BCC and OTHER phases, Δ ω j MD are the ranges of values of fractions in MD simulations and k is the normalization factor. We considered the HCP phase as a defective structure in the FCC lattice. Distribution maps were constructed to determine the accuracy, range and step of parameter variation. Figure 15 shows the parameter and probability distributions in a solid solution with 70% Cu at 300 K. It can be seen from the graphs that the set of parameters with maximum likelihood is not unambiguously defined. That is, insignificant fluctuations in the value of parameters will not significantly affect the shape of the curves of the distribution of the phase fractions. The set of parameters corresponding to the highest probability was further used to solve Equations (15) and (16). The calculated values of the phase fractions during uniaxial compression for the alloy with 70% Cu are in agreement with the MD modeling data at 300 K (Figure 16). The model describes the transition from the FCC to the BCC structure and regions with disordered atoms, as well as the further disordering of the BCC phase.

4.3. Optimization and Results of Stress Relaxation Model

The degree of agreement of shear stresses and pressures with the results of the MD calculations was estimated in order to determine the contribution of phase transitions to the stress relaxation process (parameters n 1 and n 2 ) and the parameters of the dislocation plasticity model. The shear stresses were calculated from the stress relaxation model and the pressure from the equation of state in the form of the ANN. The expression for calculating the probability has the following form:
Probability = i exp k τ i model τ i MD Δ τ MD 2 + P i model P i MD Δ P MD 2 .
where Δ τ MD and Δ P MD are the ranges of shear stress and pressure values, respectively, in MD simulations. Figure 17 shows the probability maps for different parameter ranges for an alloy with 10% Cu at 300 K. The variation in values of yield strength Y 0 and dislocation friction coefficient B almost does not affect the probability value. A similar low sensitivity of the results of the dislocation-based model of pore growth to the dislocation friction coefficient was previously reported in [115] in the case of iron. On the contrary, favorable values of the parameter n 1 are localized in a narrow range. For the hardening coefficient A, nucleation parameter k n and dislocation loop formation energy ε D , there are regions for which the probability takes zero value. In this case, the shear stress plot is far from the reference curve.
Let us consider the results of the shear stress relaxation model with the fitted parameters by the Bayesian method by an example with a copper concentration of 10% presented in Figure 18. The model qualitatively well describes the stresses in the material under uniaxial compression (Figure 18a). The difference is in the amplitude of the maximum shear stresses. The model also indicates hardening in the material at strains 0.1–0.2 and a sharper drop in stresses at the initiation of dislocation activity (at a strain of about 0.09) and phase transitions (at a strain of about 0.17). The ANN describes the pressure in the material with good accuracy (Figure 18b). Another ANN calculates the shear modulus and nucleation threshold of dislocations. After the nucleation of dislocations, there is an active multiplication of dislocations at a rate of the order of 1027 s−1m−2 (Figure 19). Due to the dynamic loading with an ultra-high strain rate applied to the material, the dislocation density exceeds 1016 m−2, which is in the range of the experimental data for Cu [116] and MD modeling [117,118]. Our model value of the dislocation density in Figure 19a is an order of magnitude lower than that determined by the MD in Figure 3, but this discrepancy is explained by the fact that we took into account only mobile dislocations in Equation (18) of the relaxation plasticity model, while most of the dislocation segments in the MD are immobilized due to intersections with other dislocations. More complex dislocation kinetics with an explicit separate accounting of mobile and immobilized dislocations similar to [63] can be applied in future work. The plastic strain rate increases as the dislocation density increases dramatically, but after an initial spike at a strain of about 0.1, it drops by an order of magnitude.

5. Conclusions

The MD study of uniaxial compression of Al, Cu and Al-Cu solid solutions was performed. The stress relaxation process is provided by (a) dislocation plasticity and (b) phase transition. There is a sharp drop in stresses in pure metals due to the initiation of plastic deformation induced by dislocation nucleation and motion. The formation of a metastable BCC phase is also observed, which precedes the onset of the plasticity process. In contrast to copper, aluminum reveals a phase transition from the FCC to the BCC structure almost throughout the entire volume of the simulated crystal at pressures of about 36 GPa. The phase transition is accompanied by a decrease in shear stresses to negative values. A similar behavior is observed for Al-Cu solid solutions, but there is a difference in that the phase transition becomes the dominant stress relaxation mechanism. The phase transition completely suppresses the dislocation activity for alloys with Cu concentrations in the range 30–80%, while the dislocation plasticity prevails again for pure Cu. An increase in temperature leads to a softening of the material and a decrease in the amplitude of the shear stresses.
A theoretical model of stress relaxation was developed, taking into account the main mechanisms of plasticity, the dislocation activity and phase transitions, by an example of an Al-Cu solid solution. The Bayesian method was applied to identify the model parameters. Two forward-propagation ANNs trained by MD data for uniaxial compression and tension were used to approximate the constitutive relations. The ANNs were used to describe the equation of state and determine the beginning of plasticity and the elastic modulus.
Further development of the present work is connected with using the developed model for the simulation of a shock-wave structure in an Al-Cu solid solution, as well as with the application of the model to other metallic systems, including high-entropy alloys. Most of the solid solutions considered here are thermodynamically unstable in traditional metallurgical casting, and experimental data are still lacking, but if they reveal unique dynamic-protecting properties, novel production techniques, such as additive manufacturing, can be applied to produce such metastable alloys. Therefore, a theoretical analysis of their properties is relevant.

Author Contributions

N.G.: methodology, software, validation, formal analysis, investigation, writing—original draft preparation and visualization; E.F.: methodology, software, validation, formal analysis, investigation, writing—original draft preparation and visualization; A.M.: conceptualization, methodology, resources, data curation, writing—review and editing, supervision, project administration and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC were funded by the RUSSIAN SCIENCE FOUNDATION, grant number 20-11-20153-P, https://rscf.ru/en/project/23-11-45024/ (accessed on 10 July 2024).

Data Availability Statement

Data are available upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The Adam algorithm can be described in terms of several calculations:
g = 1 m W , b L f x ( i ) , W , b , y ( i ) , s ( n + 1 ) = ρ 1 s ( n ) + 1 ρ 1 g , r ( n + 1 ) = ρ 2 r ( n ) + 1 ρ 2 g g .
First, the gradient of lost function and gradient momentum of the first s ( n + 1 ) and r ( n + 1 ) second orders were calculated using Formula (A1). Second, the gradient values were corrected using the formulas:
s ^ ( n + 1 ) = s ( n + 1 ) 1 ρ 1 t , r ^ ( n + 1 ) = r ( n + 1 ) 1 ρ 2 t .
Third, after correcting the gradient values using Formula (A2), the neural network parameters were updated:
Δ W = ε s ^ W ( n + 1 ) r ^ W ( n + 1 ) + δ , Δ b = ε s ^ b ( n + 1 ) r ^ b ( n + 1 ) + δ , W ( n + 1 ) = W ( n ) + Δ W .
The following parameters were used in the Adam algorithm (Equations (A1)–(A3)): ρ 1 = 0 , 1 , ρ 2 = 0 , 1 and t is the current epoch number and δ is the small parameter for numerical stabilization. The parameters W and b in Equations (A1)–(A3) denote all weights and biases in the neural network.

References

  1. Steinberg, D.J.; Cochran, S.G.; Guinan, M.W. A constitutive model for metals applicable at high strain rate. J. Appl. Phys. 1980, 51, 1498–1504. [Google Scholar] [CrossRef]
  2. Gust, W.H. High impact deformation of metal cylinders at elevated temperatures. J. Appl. Phys. 1982, 53, 3566–3575. [Google Scholar] [CrossRef]
  3. Zhao, H. A constitutive model for metals over a large range of strain rates Identification for mild-steel and aluminium sheets. Mater. Sci. Eng. A 1997, 230, 95–99. [Google Scholar] [CrossRef]
  4. Follansbee, P.S.; Kocks, U.F. A constitutive description of the deformation of copper based on the use of the mechanical threshold stress as an internal state variable. Acta Metall. 1988, 36, 81–93. [Google Scholar] [CrossRef]
  5. Gourdin, W.H.; Lassila, D.H. Flow stress of OFE copper at strain rates from 10−3 to 104 s−1: Grain-size effects and comparison to the mechanical threshold stress model. Acta Metall. Mater. 1991, 39, 2337–2348. [Google Scholar] [CrossRef]
  6. Zerilli, F.J.; Armstrong, R.W. Dislocation mechanics-based constitutive relations for material dynamics calculations. J. Appl. Phys. 1987, 61, 1816–1825. [Google Scholar] [CrossRef]
  7. Bryukhanov, I.A.; Emelyanov, V.A. Shear stress relaxation through the motion of edge dislocations in Cu and Cu–Ni solid solution: A molecular dynamics and discrete dislocation study. Comput. Mater. Sci. 2022, 201, 110885. [Google Scholar] [CrossRef]
  8. Zhang, G.; Zhu, Z.; Ning, J.; Feng, C. Dynamic impact constitutive relation of 6008-T6 aluminum alloy based on dislocation density and second-phase particle strengthening effects. J. Alloys Compd. 2023, 932, 167718. [Google Scholar] [CrossRef]
  9. Khan, A.S.; Liu, J. A deformation mechanism-based crystal plasticity model of ultrafine-grained/nanocrystalline FCC polycrystals. Int. J. Plast. 2016, 86, 56–69. [Google Scholar] [CrossRef]
  10. Austin, R.A.; McDowell, D.L. A dislocation-based constitutive model for viscoplastic deformation of fcc metals at very high strain rates. Int. J. Plast. 2011, 27, 1–24. [Google Scholar] [CrossRef]
  11. Luscher, D.J.; Mayeur, J.R.; Mourad, H.M.; Hunter, A.; Kenamond, M.A. Coupling continuum dislocation transport with crystal plasticity for application to shock loading conditions. Int. J. Plast. 2016, 76, 111–129. [Google Scholar] [CrossRef]
  12. Yao, S.-L.; Pei, X.-Y.; Yu, J.-D.; Bai, J.-S.; Wu, Q. A dislocation-based explanation of quasi-elastic release in shock-loaded aluminum. J. Appl. Phys. 2017, 121, 035101. [Google Scholar] [CrossRef]
  13. Austin, R.A. Elastic precursor wave decay in shock-compressed aluminum over a wide range of temperature. J. Appl. Phys. 2018, 123, 035103. [Google Scholar] [CrossRef]
  14. Djordjevic, N.; Vignjevic, R.; Kiely, L.; Case, S.; De Vuyst, T.; Campbell, J.; Hughes, K. Modelling of shock waves in fcc and bcc metals using a combined continuum and dislocation kinetic approach. Int. J. Plast. 2018, 105, 211–224. [Google Scholar] [CrossRef]
  15. Yao, S.; Pei, X.; Liu, Z.; Yu, J.; Yu, Y.; Wu, Q. Numerical investigation of the temperature dependence of dynamic yield stress of typical BCC metals under shock loading with a dislocation-based constitutive model. Mech. Mater. 2020, 140, 103211. [Google Scholar] [CrossRef]
  16. Mori, T.; Huang, J.H.; Taya, M. Stress relaxation by plastic flow, interfacial sliding and diffusion in an inclusion bearing material. Acta Mater. 1997, 45, 429–438. [Google Scholar] [CrossRef]
  17. Borodin, E.N.; Mayer, A.E. Structural model of mechanical twinning and its application for modeling of the severe plastic deformation of copper rods in Taylor impact tests. Int. J. Plast. 2015, 74, 141–157. [Google Scholar] [CrossRef]
  18. Li, P.; Huang, Y.; Wang, K.; Xiao, S.; Wang, L.; Yao, S.; Zhu, W.; Hu, W. Crystallographic-orientation-dependence plasticity of niobium under shock compressions. Int. J. Plast. 2022, 150, 103195. [Google Scholar] [CrossRef]
  19. Feng, X.; Liu, X.; Bai, S.; Tang, Y.; Ye, Y. Investigation of dynamic tensile mechanical responses and deformation mechanism at high strain rates in a TWIP steel. J. Mater. Res. Technol. 2023, 26, 639–653. [Google Scholar] [CrossRef]
  20. Zhang, M.; Ma, A.; Lorca, J. Simulation of dislocation slip and twin propagation in Mg through coupling crystal plasticity and phase field models. Mech. Mater. 2023, 180, 104608. [Google Scholar] [CrossRef]
  21. Fomin, E.V.; Mayer, A.E. Slip of low-angle tilt grain boundary (110) in FCC metals at perpendicular shear. Int. J. Plast. 2020, 134, 102843. [Google Scholar] [CrossRef]
  22. Huang, Q.; Zhao, Q.; Zhou, H.; Yang, W. Misorientation-dependent transition between grain boundary migration and sliding in FCC metals. Int. J. Plast. 2022, 159, 103466. [Google Scholar] [CrossRef]
  23. Verma, A.; Johnson, O.K.; Thompson, G.B.; Chesser, I.; Ogata, S.; Homer, E.R. Insights into factors that affect non-Arrhenius migration of a simulated incoherent Σ3 grain boundary. Acta Mater. 2023, 258, 119210. [Google Scholar] [CrossRef]
  24. Admal, N.C.; Po, G.; Marian, J. A unified framework for polycrystal plasticity with grain boundary evolution. Int. J. Plast. 2018, 106, 1–30. [Google Scholar] [CrossRef]
  25. Masuda, H.; Tobe, H.; Hara, T.; Sato, E. Three-dimensional characterization of superplastic grain boundary sliding inside Al–Zn–Mg–Cu alloy sheet. Scr. Mater. 2019, 164, 82–85. [Google Scholar] [CrossRef]
  26. Borodin, E.N.; Mayer, A.E.; Gutkin, M.Y. Coupled model for grain rotation, dislocation plasticity and grain boundary sliding in fine-grained solids. Int. J. Plast. 2020, 134, 102776. [Google Scholar] [CrossRef]
  27. Li, G.; Wang, Y.; Wang, K.; Xiang, M.; Chen, J. Shock induced plasticity and phase transition in single crystal lead by molecular dynamics simulations. J. Appl. Phys. 2019, 126, 075902. [Google Scholar] [CrossRef]
  28. Jiang, D.-D.; Shao, J.-L.; Wu, B.; Wang, P.; He, A.-M. Sudden change of spall strength induced by shock defects based on atomistic simulation of single crystal aluminum. Scr. Mater. 2022, 210, 114474. [Google Scholar] [CrossRef]
  29. Zhu, Y.; Hu, J.; Wei, Q.; Zhang, J.; Sun, Y.; Luo, G.; Shen, Q. Enhanced spall strength of single crystal aluminum by temperature rise mitigation and structural phase transition under shock pulse. Mech. Mater. 2023, 186, 104809. [Google Scholar] [CrossRef]
  30. Bolesta, A.V.; Fomin, V.M. Molecular dynamics simulation of shock-wave loading of copper and titanium. AIP Conf. Proc. 2017, 1893, 020008. [Google Scholar] [CrossRef]
  31. Shao, J.-L.; Guo, X.-X.; Lu, G.; He, W.; Xin, J. Influence of shear wave on the HCP nucleation in BCC iron under oblique shock conditions. Mech. Mater. 2021, 158, 103878. [Google Scholar] [CrossRef]
  32. Zhao, K.; Gu, Y.; Zhao, F.; Chen, Q. Revisiting the effect of shear stress on the γ→α phase transition of cerium under shock loading. Mech. Mater. 2023, 184, 104743. [Google Scholar] [CrossRef]
  33. Jian, Z.; Chen, Y.; Xiao, S.; Wang, L.; Li, X.; Wang, K.; Deng, H.; Hu, W. Shock-induced plasticity and phase transformation in single crystal magnesium: An interatomic potential and non-equilibrium molecular dynamics simulations. J. Phys. Condens. Matter 2022, 34, 115401. [Google Scholar] [CrossRef] [PubMed]
  34. Wu, Y.-C.; Shao, J.-L. FCC-BCC phase transformation induced simultaneous enhancement of tensile strength and ductility at high strain rate in high-entropy alloy. Int. J. Plast. 2023, 169, 103730. [Google Scholar] [CrossRef]
  35. Xie, H.; Ma, Z.; Zhang, W.; Zhao, H.; Ren, L. Phase transition in shock compressed high-entropy alloy FeNiCrCoCu. Int. J. Mech. Sci. 2023, 238, 107855. [Google Scholar] [CrossRef]
  36. Dewaele, A.; Loubeyre, P.; Mezouar, M. Equations of state of six metals above 94 GPa. Phys. Rev. B 2004, 70, 094112. [Google Scholar] [CrossRef]
  37. Sharma, S.M.; Turneaure, S.J.; Winey, J.M.; Gupta, Y.M. Transformation of shock-compressed copper to the body-centered-cubic structure at 180 GPa. Phys. Rev. B 2020, 102, 020103. [Google Scholar] [CrossRef]
  38. Fratanduono, D.E.; Smith, R.F.; Ali, S.J.; Braun, D.G.; Fernandez-Pañella, A.; Zhang, S.; Kraus, R.G.; Coppari, F.; McNaney, J.M.; Marshall, M.C.; et al. Probing the solid phase of noble metal copper at terapascal conditions. Phys. Rev. Lett. 2020, 124, 015701. [Google Scholar] [CrossRef]
  39. Javanbakht, M.; Levitas, V.I. Interaction between phase transformations and dislocations at the nanoscale. Part 2: Phase field simulation examples. J. Mech. Phys. Solids 2015, 82, 164–185. [Google Scholar] [CrossRef]
  40. Levitas, V.I.; Javanbakht, M. Interaction between phase transformations and dislocations at the nanoscale. Part 1. General phase field approach. J. Mech. Phys. Solids. 2015, 82, 287–319. [Google Scholar] [CrossRef]
  41. Yeddu, H.K.; Lookman, T. Phase-field modeling of the beta to omega phase transformation in Zr-Nb alloys. Mater. Sci. Eng. A 2015, 634, 46–54. [Google Scholar] [CrossRef]
  42. Yang, H.; Zhou, G.; Wang, H.; Wu, P.; Lee, M.-G.; Li, D. Experimental and numerical study on strain path and crystal orientation dependences of phase transformation behaviors of QP1180 steel. Int. J. Plast. 2023, 160, 103499. [Google Scholar] [CrossRef]
  43. Binqiang, L.; Guiji, W.; Jianjun, M.; Hongpin, Z.; Fuli, T.; Jianheng, Z.; Cangli, L.; Chengwei, S. Verification of conventional equations of state for tantalum under quasi-isentropic compression. J. Appl. Phys. 2014, 116, 193506. [Google Scholar] [CrossRef]
  44. Shang, S.-L.; Wang, Y.; Kim, D.E.; Liu, Z.-K. First-principles thermodynamics from phonon and Debye model: Application to Ni and Ni3Al. Comput. Mater. Sci. 2010, 47, 1040–1048. [Google Scholar] [CrossRef]
  45. Matar, S.F.; Al Alam, A.F.; Ouaini, N. First principles account for large changes in electronic structure and bonding from LaCu to LaCuMg and LaCuMg4. Comput. Mater. Sci. 2015, 97, 231–236. [Google Scholar] [CrossRef]
  46. Gracheva, N.A.; Lekanov, M.V.; Mayer, A.E.; Fomin, E.V. Application of neural networks for modeling shock-wave processes in aluminum. Mech. Solids 2021, 56, 326–342. [Google Scholar] [CrossRef]
  47. Bonatti, C.; Berisha, B.; Mohr, D. From CP-FFT to CP-RNN: Recurrent neural network surrogate model of crystal plasticity. Int. J. Plast. 2022, 58, 103430. [Google Scholar] [CrossRef]
  48. Heidenreich, J.N.; Gorji, M.B.; Mohr, D. Modeling structure-property relationships with convolutional neural networks: Yield surface prediction based on microstructure images. Int. J. Plast. 2023, 163, 103506. [Google Scholar] [CrossRef]
  49. Casalino, G.; Minutolo, F.M.C. A model for evaluation of laser welding efficiency and quality using an artificial neural network and fuzzy logic. Proc. Inst. Mech. Eng. B J. Eng. Manuf. 2004, 218, 641–646. [Google Scholar] [CrossRef]
  50. Du, Y.; Mukherjee, T.; DebRoy, T. Conditions for void formation in friction stir welding from machine learning. NPJ Comput. Mater. 2019, 5, 68. [Google Scholar] [CrossRef]
  51. Gaja, H.; Liou, F. Defect classification of laser metal deposition using logistic regression and artificial neural networks for pattern recognition. Int. J. Adv. Manuf. Technol. 2017, 94, 315–326. [Google Scholar] [CrossRef]
  52. Li, X.; Roth, C.C.; Mohr, D. Machine-learning based temperature- and rate-dependent plasticity model: Application to analysis of fracture experiments on DP steel. Int. J. Plast. 2019, 118, 320–344. [Google Scholar] [CrossRef]
  53. Li, X.; Roth, C.C.; Mohr, D. Neural network based rate- and temperature-dependent Hosford–Coulomb fracture initiation model. Int. J. Mech. Sci. 2023, 260, 108643. [Google Scholar] [CrossRef]
  54. Nakka, R.; Harursampath, D.; Ponnusami, S.A. A generalised deep learning-based surrogate model for homogenisation utilising material property encoding and physics-based bounds. Sci. Rep. 2023, 13, 9079. [Google Scholar] [CrossRef]
  55. Zhu, L.; Wang, X.; Zhong, W.; Li, M. Learning microstructure–property mapping via label-free 3D convolutional neural network. Vis. Comput. 2024. [Google Scholar] [CrossRef]
  56. Latypov, F.T.; Fomin, E.V.; Krasnikov, V.S.; Mayer, A.E. Dynamic compaction of aluminum with nanopores of varied shape: MD simulations and machine-learning-based approximation of deformation behavior. Int. J. Plast. 2022, 156, 103363. [Google Scholar] [CrossRef]
  57. Latypov, F.T.; Fomin, E.V.; Krasnikov, V.S.; Mayer, A.E. Modeling of shock wave propagation in porous magnesium based on artificial neural network. Mech. Mater. 2024, 191, 104953. [Google Scholar] [CrossRef]
  58. Walters, D.J.; Biswas, A.; Lawrence, E.C.; Francom, D.C.; Luscher, D.J.; Fredenburg, D.A.; Moran, K.R.; Sweeney, C.M.; Sandberg, R.L.; Ahrens, J.P.; et al. Bayesian calibration of strength parameters using hydrocode simulations of symmetric impact shock experiments of Al-5083. J. Appl. Phys. 2018, 124, 205105. [Google Scholar] [CrossRef]
  59. Sjue, S.; Ahrens, J.; Biswas, A.; Francom, D.; Lawrence, E.; Luscher, D.; Walters, D. Fast strength model characterization using Bayesian statistics. AIP Conf. Proc. 2020, 2272, 070043. [Google Scholar] [CrossRef]
  60. Nguyen, T.; Francom, D.C.; Luscher, D.J.; Wilkerson, J.W. Bayesian calibration of a physics-based crystal plasticity and damage model. J. Mech. Phys. Solids 2021, 149, 104284. [Google Scholar] [CrossRef]
  61. Nguyen, T.; Fensin, S.J.; Luscher, D.J. Dynamic crystal plasticity modeling of single crystal tantalum and validation using Taylor cylinder impact tests. Int. J. Plast. 2021, 139, 102940. [Google Scholar] [CrossRef]
  62. Rodionov, E.S.; Pogorelko, V.V.; Lupanov, V.G.; Mayer, P.N.; Mayer, A.E. Modified Taylor impact tests with profiled copper cylinders: Experiment and optimization of dislocation plasticity model. Materials 2023, 16, 5602. [Google Scholar] [CrossRef] [PubMed]
  63. Popova, T.V.; Mayer, A.E.; Khishchenko, K.V. Evolution of shock compression pulses in polymethacrylate and aluminum. J. Appl. Phys. 2018, 123, 235902. [Google Scholar] [CrossRef]
  64. Ashitkov, S.I.; Komarov, P.S.; Struleva, E.V.; Agranat, M.B.; Kanel, G.I. Mechanical and optical properties of vanadium under shock picosecond loads. JETP Lett. 2015, 101, 276–281. [Google Scholar] [CrossRef]
  65. Kanel, G.I.; Zaretsky, E.B.; Razorenov, S.V.; Ashitkov, S.I.; Fortov, V.E. Unusual plasticity and strength of metals at ultra-short load durations. Phys. Usp. 2017, 60, 490–508. [Google Scholar] [CrossRef]
  66. Zuanetti, B.; McGrane, S.D.; Bolme, C.A.; Prakash, V. Measurement of elastic precursor decay in pre-heated aluminum films under ultra-fast laser generated shocks. J. Appl. Phys. 2018, 123, 195104. [Google Scholar] [CrossRef]
  67. Khokhlov, V.A.; Zhakhovsky, V.V.; Inogamov, N.A.; Ashitkov, S.I.; Sitnikov, D.S.; Khishchenko, K.V.; Petrov, Y.V.; Manokhin, S.S.; Nelasov, I.V.; Shepelev, V.V.; et al. Melting of titanium by a shock wave generated by an intense femtosecond laser pulse. JETP Lett. 2022, 115, 523–530. [Google Scholar] [CrossRef]
  68. Jax, P.; Kratochvil, P.; Haasen, P. Solid solution hardening of gold and other f.c.c. single crystals. Acta Metall. 1970, 18, 237–245. [Google Scholar] [CrossRef]
  69. Monchoux, F.; Neuhäuser, H. An in situ transmission electron microscope study of dislocation glide in Cu-Ge alloys. J. Mater. Sci. 1987, 22, 1443–1449. [Google Scholar] [CrossRef]
  70. Roth, H.A.; Davis, C.L.; Thomson, R.C. Modeling solid solution strengthening in nickel alloys. Metall. Mater. Trans. A 1997, 28, 1329–1335. [Google Scholar] [CrossRef]
  71. Caillard, D.; Martin, J.L. Thermally Activated Mechanisms in Crystal Plasticity; Pergamon Press: Oxford, UK, 2003; Volume 8, pp. 3–433. [Google Scholar]
  72. Xu, W.; Liu, X.C.; Li, X.Y.; Lu, K. Deformation induced grain boundary segregation in nanolaminated Al–Cu alloy. Acta Mater. 2020, 182, 207–214. [Google Scholar] [CrossRef]
  73. Sun, W.; Zhu, Y.; Marceau, R.; Wang, L.; Zhang, Q.; Gao, X.; Hutchinson, C. Precipitation strengthening of aluminum alloys by room-temperature cyclic plasticity. Science 2019, 363, 972–975. [Google Scholar] [CrossRef] [PubMed]
  74. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef]
  75. Apostol, F.; Mishin, Y. Interatomic potential for the Al-Cu system. Phys. Rev. B 2011, 83, 054116. [Google Scholar] [CrossRef]
  76. Hoover, W.G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695–1697. [Google Scholar] [CrossRef] [PubMed]
  77. Larsen, P.M.; Schmidt, S.; Schiøtz, J. Robust structural identification via polyhedral template matching. Model. Simul. Mat. Sci. Eng. 2016, 24, 055007. [Google Scholar] [CrossRef]
  78. Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO—The Open Visualization Tool. Model. Simul. Mater. Sci. Eng. 2010, 18, 015012. [Google Scholar] [CrossRef]
  79. Mayer, A.E.; Mayer, P.N.; Lekanov, M.V.; Panchenko, B.A. Incipience of plastic flow in aluminum with nanopores: Molecular dynamics and machine-learning-based description. Metals 2022, 12, 2158. [Google Scholar] [CrossRef]
  80. Grachyova, N.A.; Fomin, E.V.; Mayer, A.E. Application of Machine Learning in the Study of Phase Transformations of Al-Cu Alloys under Deformation. Available online: https://www.researchgate.net/publication/366432358_APPLICATION_OF_MACHINE_LEARNING_IN_THE_STUDY_OF_DEFORMATION_BEHAVIOR_OF_AL-CU_ALLOYS (accessed on 10 July 2024).
  81. Grachyova, N.A.; Fomin, E.V.; Mayer, A.E. Theoretical description of phase transitions in Al-Cu alloys using machine learning. In XIII All-Russian Congress on Theoretical and Applied Mechanics. Collection of Abstracts: In 4 Vol.; Ministry of Science and Higher Education of the Russian Federation; Russian Academy of Sciences; Russian National Committee on Theoretical and Applied Mechanics; Peter the Great St. Petersburg Polytechnic University: St. Petersburg, Russia, 2023; pp. 358–359. [Google Scholar]
  82. Grachyova, N.A.; Fomin, E.V.; Mayer, A.E. Application of machine learning in the study of deformation behaviour of Al-Cu alloys, Physical Mesomechanics of Materials. In Physical Principles of Multilevel Structure Formation and Mechanisms of Nonlinear Behaviour: Proc. of Intern. Conf., Tomsk, Russia, 5–8 September 2022; Institute of Strength Physics and Materials Science SB RAS; IPC NSU: Novosibirsk, Russia, 2022; p. 556. ISBN 978-5-4437-1353-3. [Google Scholar]
  83. Hashiguchi, K. Elastoplasticity Theory; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar] [CrossRef]
  84. Mayer, A.E.; Krasnikov, V.S.; Pogorelko, V.V. Homogeneous nucleation of dislocations in copper: Theory and approximate description based on molecular dynamics and artificial neural networks. Comput. Mater. Sci. 2022, 206, 111266. [Google Scholar] [CrossRef]
  85. Goodfellow, I.; Bengio, Y.; Courville, A. Deep Learning; MIT Press: Cambridge, MA, USA, 2016. [Google Scholar]
  86. He, K.; Zhang, X.; Ren, S.; Sun, J. Delving deep into rectifiers: Surpassing human-level performance on imagenet classification. In Proceedings of the 2015 IEEE International Conference on Computer Vision (ICCV), Santiago, Chile, 7–13 December 2015; pp. 1026–1034. [Google Scholar] [CrossRef]
  87. Ramachandran, P.; Zoph, B.; Le, Q.V. Searching for Activation Functions. arXiv 2017, arXiv:1710.05941. [Google Scholar] [CrossRef]
  88. Pogorelko, V.V.; Mayer, A.E.; Fomin, E.V.; Fedorov, E.V. Examination of machine learning method for identification of material model parameters. Int. J. Mech. Sci. 2024, 265, 108912. [Google Scholar] [CrossRef]
  89. Kingma, D.P.; Ba, J. Adam: A method for stochastic optimization. arXiv 2014, arXiv:1412.6980. [Google Scholar] [CrossRef]
  90. Sutton, P.M. The variation of the elastic constants of crystalline aluminum with temperature between 63 °K and 773 °K. Phys. Rev. 1953, 91, 816–821. [Google Scholar] [CrossRef]
  91. Kamm, G.N.; Alers, G.A. Low temperature elastic moduli of aluminum. J. Appl. Phys. 1964, 35, 327–330. [Google Scholar] [CrossRef]
  92. Gerlich, D.; Fisher, E.S. The high temperature elastic moduli of aluminum. J. Phys. Chem. Solids 1969, 30, 1197–1205. [Google Scholar] [CrossRef]
  93. Tallon, J.L.; Wolfenden, A. Temperature dependence of the elastic constants of aluminum. J. Phys. Chem. Solids 1979, 40, 831–837. [Google Scholar] [CrossRef]
  94. Overton, W.C., Jr.; Gaffney, J. Temperature variation of the elastic constants of cubic elements. I. Copper. Phys. Rev. 1955, 98, 969–977. [Google Scholar] [CrossRef]
  95. Chang, Y.A.; Himmel, L. Temperature dependence of the elastic constants of Cu, Ag, and Au above room temperature. J. Appl. Phys. 1966, 37, 3567–3572. [Google Scholar] [CrossRef]
  96. Ledbetter, H.M.; Naimon, E.R. Elastic properties of metals and alloys. II. Copper. J. Phys. Chem. Ref. Data. 1974, 3, 897–935. [Google Scholar] [CrossRef]
  97. Raju, S.; Sivasubramanian, K.; Mohandas, E. The high temperature bulk modulus of aluminium: An assessment using experimental enthalpy and thermal expansion data. Solid State Commun. 2002, 122, 671–676. [Google Scholar] [CrossRef]
  98. Wawra, H. Die Kröner-Grenzen der Elastizitätsmoduln technisch wichtiger Werkstoffe: Teil II: Zahlenwerte zur Temperaturabhängigkeit der Moduln der Elemente. Int. J. Mater. Res. 1978, 69, 518–523. [Google Scholar] [CrossRef]
  99. Wen, Y.; Yunxin, W. Effects of temperature and pressure on elastic properties of single crystal aluminum in different crystal orientations. Phys. Status Solidi B 2020, 257, 2000434. [Google Scholar] [CrossRef]
  100. Marchand, D.; Curtin, W.A. Machine learning for metallurgy IV: A neural network potential for Al-Cu-Mg and Al-Cu-Mg-Zn. Phys. Rev. Mater. 2022, 6, 053803. [Google Scholar] [CrossRef]
  101. Purja Pun, G.P.; Batra, R.; Ramprasad, R.; Mishin, Y. Physically informed artificial neural networks for atomistic modeling of materials. Nat. Commun. 2019, 10, 2339. [Google Scholar] [CrossRef]
  102. Senoo, M.; Hayashi, T. Elastic constants of Al-Cu solid-solution alloys and their variations by aging treatments. JSME Int. J. Ser. 1 Solid Mech. Strength Mater. 1988, 31, 664–670. [Google Scholar] [CrossRef]
  103. Ma, D. First-Principles Investigations of Solid Solution Strengthening in Al Alloys; Shaker Verlag GmbH: Aachen, Germany, 2012. [Google Scholar]
  104. Lee, E.H.; Liu, D.T. Finite–strain elastic–plastic theory with application to plane wave analysis. J. Appl. Phys. 1967, 38, 19–27. [Google Scholar] [CrossRef]
  105. Lee, E.H. Elastic-plastic deformation at finite strains. J. Appl. Mech. 1969, 36, 1–6. [Google Scholar] [CrossRef]
  106. Nemat-Nasser, S. Decomposition of strain measures and their rates in finite deformation elastoplasticity. Int. J. Solids Struct. 1979, 15, 155–166. [Google Scholar] [CrossRef]
  107. Zhao, S.; Petrov, Y.V.; Volkov, G.A. The modified relaxation plasticity model and the non-monotonic stress–strain diagram. Int. J. Mech. Sci. 2023, 240, 107919. [Google Scholar] [CrossRef]
  108. Selyutina, N.; Borodin, E.N.; Petrov, Y.; Mayer, A.E. The definition of characteristic times of plastic relaxation by dislocation slip and grain boundary sliding in copper and nickel. Int. J. Plast. 2016, 82, 97–111. [Google Scholar] [CrossRef]
  109. Mayer, A.E.; Lekanov, M.V.; Grachyova, N.A.; Fomin, E.V. Machine-learning-based model of elastic—Plastic deformation of copper for application to shock wave problem. Metals 2022, 12, 402. [Google Scholar] [CrossRef]
  110. Krasnikov, V.S.; Mayer, A.E. Influence of local stresses on motion of edge dislocation in aluminum. Int. J. Plast. 2018, 101, 170–187. [Google Scholar] [CrossRef]
  111. Kuksin, A.Y.; Yanilkin, A.V. Atomistic simulation of the motion of dislocations in metals under phonon drag conditions. Phys. Solid State 2013, 55, 1010–1019. [Google Scholar] [CrossRef]
  112. Krasnikov, V.S.; Mayer, A.E.; Yalovets, A.P. Dislocation based high-rate plasticity model and its application to plate-impact and ultra short electron irradiation simulations. Int. J. Plast. 2011, 27, 1294–1308. [Google Scholar] [CrossRef]
  113. Landau, L.D.; Lifshitz, E.M. Theory of Elasticity. In Course of Theoretical Physics; Pergamon Press: New York, NY, USA, 1986; Volume 7. [Google Scholar]
  114. Tohme, T.; Vanslette, K.; Youcef-Toumi, K. A generalized Bayesian approach to model calibration. Reliab. Eng. Syst. Saf. 2020, 204, 107141. [Google Scholar] [CrossRef]
  115. Pogorelko, V.V.; Mayer, A.E. Dynamic tensile fracture of iron: Molecular dynamics simulations and micromechanical model based on dislocation plasticity. Int. J. Plast. 2023, 167, 103678. [Google Scholar] [CrossRef]
  116. Meyers, M.A.; Gregori, F.; Kad, B.K.; Schneider, M.S.; Kalantar, D.H.; Remington, B.A.; Ravichandran, G.; Boehly, T.; Wark, J.S. Laser-induced shock compression of monocrystalline copper: Characterization and analysis. Acta Mater. 2003, 51, 1211–1228. [Google Scholar] [CrossRef]
  117. Bringa, E.M.; Rosolankova, K.; Rudd, R.E.; Remington, B.A.; Wark, J.S.; Duchaineau, M.; Kalantar, D.H.; Hawreliak, J.; Belak, J. Shock deformation of face-centred-cubic metals on subnanosecond timescales. Nat. Mater. 2006, 5, 805–809. [Google Scholar] [CrossRef]
  118. Hu, J.; Liu, Z.; Chen, K.; Zhuang, Z. Investigations of shock-induced deformation and dislocation mechanism by a multiscale discrete dislocation plasticity model. Comput. Mater. Sci. 2017, 131, 78–85. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of crystal loading in MD modeling. The compression direction [100] is shown by arrows.
Figure 1. Schematic representation of crystal loading in MD modeling. The compression direction [100] is shown by arrows.
Dynamics 04 00028 g001
Figure 2. Shear stresses and evolution of phase fractions for pure aluminum (a) and pure copper (b) during compression.
Figure 2. Shear stresses and evolution of phase fractions for pure aluminum (a) and pure copper (b) during compression.
Dynamics 04 00028 g002
Figure 3. Dislocation density evolution for (a) pure aluminum and copper and (b) Al-Cu solid solutions during compression at the temperature of 300 K obtained from the MD. The dislocation length in FCC phase was divided by the total volume of the system; therefore, it tends to zero after the phase transition from FCC to BCC and “Other” structures.
Figure 3. Dislocation density evolution for (a) pure aluminum and copper and (b) Al-Cu solid solutions during compression at the temperature of 300 K obtained from the MD. The dislocation length in FCC phase was divided by the total volume of the system; therefore, it tends to zero after the phase transition from FCC to BCC and “Other” structures.
Dynamics 04 00028 g003
Figure 4. Comparison of MD data on pressure–density curves for (a) aluminum and (b) copper with experimental data [36] (Dewaele (2004)) and DFT calculations [79] (DFT Panchenko (2022)).
Figure 4. Comparison of MD data on pressure–density curves for (a) aluminum and (b) copper with experimental data [36] (Dewaele (2004)) and DFT calculations [79] (DFT Panchenko (2022)).
Dynamics 04 00028 g004
Figure 5. Evolution of phase fractions and shear stresses in a solid solution with 10% Cu in aluminum. The atoms are colored according to the legend in the graph; the atoms forming FCC structure are not shown. Adapted from ref. [80].
Figure 5. Evolution of phase fractions and shear stresses in a solid solution with 10% Cu in aluminum. The atoms are colored according to the legend in the graph; the atoms forming FCC structure are not shown. Adapted from ref. [80].
Dynamics 04 00028 g005
Figure 6. Shear stresses and evolution of phase fractions of solid solutions with Cu concentration of (a) 20%, (b) 30%, (c) 50% and (d) 80%. Adapted from refs. [81,82].
Figure 6. Shear stresses and evolution of phase fractions of solid solutions with Cu concentration of (a) 20%, (b) 30%, (c) 50% and (d) 80%. Adapted from refs. [81,82].
Dynamics 04 00028 g006
Figure 7. (a) Distribution of shear stresses at different temperatures in solid solution with 20% Cu and (b) evolution of crystal structure of solid solution with 10% Cu at different strain rates.
Figure 7. (a) Distribution of shear stresses at different temperatures in solid solution with 20% Cu and (b) evolution of crystal structure of solid solution with 10% Cu at different strain rates.
Dynamics 04 00028 g007
Figure 8. The first part of constitutive equation of Al-Cu solid solution in the form of fully connected forward-propagation ANN that calculates the pressure, total energy and bulk modulus.
Figure 8. The first part of constitutive equation of Al-Cu solid solution in the form of fully connected forward-propagation ANN that calculates the pressure, total energy and bulk modulus.
Dynamics 04 00028 g008
Figure 9. The second part of constitutive equation of Al-Cu solid solution in the form of fully connected forward-propagation ANN that calculates the shear modulus and nucleation strain distance function.
Figure 9. The second part of constitutive equation of Al-Cu solid solution in the form of fully connected forward-propagation ANN that calculates the shear modulus and nucleation strain distance function.
Dynamics 04 00028 g009
Figure 10. The dependence of pressure on density for intermediate values of copper concentrations predicted by the ANN.
Figure 10. The dependence of pressure on density for intermediate values of copper concentrations predicted by the ANN.
Dynamics 04 00028 g010
Figure 11. Temperature dependencies of the bulk modulus for (a) pure copper and (b) pure aluminum predicted by the trained ANN and from experimental data and theoretical calculations: The data of Bridgman (1923), Frederick (1947) and Gruneisen (1910) are taken from [96]; Chang–Himmel (1966) [95]; Gerlich–Fisher (1969) [92]; Kamm–Alers (1964) and Tallon–Wolfenden (1979) [91,93]; Sutton (1953) [90]; Wawra (1978) [98]; Raju (2002) [97].
Figure 11. Temperature dependencies of the bulk modulus for (a) pure copper and (b) pure aluminum predicted by the trained ANN and from experimental data and theoretical calculations: The data of Bridgman (1923), Frederick (1947) and Gruneisen (1910) are taken from [96]; Chang–Himmel (1966) [95]; Gerlich–Fisher (1969) [92]; Kamm–Alers (1964) and Tallon–Wolfenden (1979) [91,93]; Sutton (1953) [90]; Wawra (1978) [98]; Raju (2002) [97].
Dynamics 04 00028 g011
Figure 12. Temperature dependencies of the shear modulus for (a) pure copper and (b) pure aluminum predicted by the trained ANN and from experimental data: Chang–Himmel (1966) [95]; Overton–Gaffney (1955) [94]; Tallon–Wolfenden (1979) [93]; Gerlich–Fisher (1969) [92]; Kamm–Alers (1964) [91].
Figure 12. Temperature dependencies of the shear modulus for (a) pure copper and (b) pure aluminum predicted by the trained ANN and from experimental data: Chang–Himmel (1966) [95]; Overton–Gaffney (1955) [94]; Tallon–Wolfenden (1979) [93]; Gerlich–Fisher (1969) [92]; Kamm–Alers (1964) [91].
Dynamics 04 00028 g012
Figure 13. The comparison of density of the Al-Cu alloy with solid solution of copper atoms obtained in this work and the experiment of Senoo–Hayashi (1988) [102].
Figure 13. The comparison of density of the Al-Cu alloy with solid solution of copper atoms obtained in this work and the experiment of Senoo–Hayashi (1988) [102].
Dynamics 04 00028 g013
Figure 14. Bulk modulus (a) and shear modulus (b) of Al-Cu solid solution for low concentrations of copper atoms in the aluminum matrix predicted by the ANN and compared with the experimental data at room temperature from Senoo–Hayashi (1988) [102] and ab initio calculations from Ma (2012) [103].
Figure 14. Bulk modulus (a) and shear modulus (b) of Al-Cu solid solution for low concentrations of copper atoms in the aluminum matrix predicted by the ANN and compared with the experimental data at room temperature from Senoo–Hayashi (1988) [102] and ab initio calculations from Ma (2012) [103].
Dynamics 04 00028 g014
Figure 15. Probability distribution in the parameter space in the case of 70% Cu for the following couples of parameters: (a) ε 23 ε 12 and (b) ε 11 U 1 .
Figure 15. Probability distribution in the parameter space in the case of 70% Cu for the following couples of parameters: (a) ε 23 ε 12 and (b) ε 11 U 1 .
Dynamics 04 00028 g015
Figure 16. Comparison of the evolution of phase fractions with parameters fitted by Bayesian method with the results of MD calculations for an alloy with 70% Cu, where the sum of FCC and HCP phases is in green, BCC phase is in blue and disordered (OTHER) structure is in grey.
Figure 16. Comparison of the evolution of phase fractions with parameters fitted by Bayesian method with the results of MD calculations for an alloy with 70% Cu, where the sum of FCC and HCP phases is in green, BCC phase is in blue and disordered (OTHER) structure is in grey.
Dynamics 04 00028 g016
Figure 17. Probability distribution in the parameter space in the case of 10% copper for the following couples of parameters: (a) dislocation formation energy–hardening coefficient, (b) dislocation formation energy–initial yield strength, (c) dislocation friction coefficient–dislocation nucleation parameter and (d) dislocation friction coefficient–parameter of shear stress relaxation due to FCC-to-BCC phase transition.
Figure 17. Probability distribution in the parameter space in the case of 10% copper for the following couples of parameters: (a) dislocation formation energy–hardening coefficient, (b) dislocation formation energy–initial yield strength, (c) dislocation friction coefficient–dislocation nucleation parameter and (d) dislocation friction coefficient–parameter of shear stress relaxation due to FCC-to-BCC phase transition.
Dynamics 04 00028 g017
Figure 18. Comparison of (a) shear stresses and (b) pressures from MD simulations with the results from a stress relaxation model with parameters fitted by Bayesian method for a copper concentration of 10% at 300 K.
Figure 18. Comparison of (a) shear stresses and (b) pressures from MD simulations with the results from a stress relaxation model with parameters fitted by Bayesian method for a copper concentration of 10% at 300 K.
Dynamics 04 00028 g018
Figure 19. (a) Plastic strain rate and dislocation density and (b) nucleation, dislocation multiplication and annihilation rates during uniaxial compression of Al-Cu solid solution with 10% Cu at 300 K. The dislocation density ρ D in FCC phase is multiplied by the phase fraction ω F C C in (a) in order to obtain an averaged density through the whole material volume.
Figure 19. (a) Plastic strain rate and dislocation density and (b) nucleation, dislocation multiplication and annihilation rates during uniaxial compression of Al-Cu solid solution with 10% Cu at 300 K. The dislocation density ρ D in FCC phase is multiplied by the phase fraction ω F C C in (a) in order to obtain an averaged density through the whole material volume.
Dynamics 04 00028 g019
Table 1. Main characteristics and hyperparameters of the ANNs and the distribution of training and test data sets.
Table 1. Main characteristics and hyperparameters of the ANNs and the distribution of training and test data sets.
ParametersANN 1ANN 2
a p r , β 0.10.6
Number of hidden layers56
Number of neurons in each hidden layer2020
μ (initial step of Adam)0.0010.001
Number of epochs420780
Mini batch size3010
n (averaging parameter for calculating shear and bulk moduli)150150
Amount of training data230,200 (80.5%)109,207 (81.6%)
Amount of test data55,515 (19.5%)24,556 (18.4%)
α C U   min / α C U   max (%)0/100
ρ min / ρ max (kg/m3)2160/12,700
Tmin/Tmax (K)100/900
Pmin/P max (GPa)−13.7/116
Emin/Emax (eV/atom)−3.60/−2.68
Kmin/Kmax (GPa)10.4/320.6
Gmin/Gmax (GPa)5.3/65.7
Q ε   min / Q ε   max 0.12/0.3
Table 2. The mean absolute percentage error (MAPE) in the final configuration of both ANNs.
Table 2. The mean absolute percentage error (MAPE) in the final configuration of both ANNs.
ANN OutputType of Data SetMAPE (%)
Pressuretraining data0.21
test data0.23
Total energytraining data0.26
test data0.35
Bulk modulustraining data0.32
test data0.44
Shear modulustraining data0.54
test data1.0
Nucleation distancetraining data0.52
test data0.9
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Grachyova, N.; Fomin, E.; Mayer, A. Theoretical Model of Structural Phase Transitions in Al-Cu Solid Solutions under Dynamic Loading Using Machine Learning. Dynamics 2024, 4, 526-553. https://doi.org/10.3390/dynamics4030028

AMA Style

Grachyova N, Fomin E, Mayer A. Theoretical Model of Structural Phase Transitions in Al-Cu Solid Solutions under Dynamic Loading Using Machine Learning. Dynamics. 2024; 4(3):526-553. https://doi.org/10.3390/dynamics4030028

Chicago/Turabian Style

Grachyova, Natalya, Eugenii Fomin, and Alexander Mayer. 2024. "Theoretical Model of Structural Phase Transitions in Al-Cu Solid Solutions under Dynamic Loading Using Machine Learning" Dynamics 4, no. 3: 526-553. https://doi.org/10.3390/dynamics4030028

APA Style

Grachyova, N., Fomin, E., & Mayer, A. (2024). Theoretical Model of Structural Phase Transitions in Al-Cu Solid Solutions under Dynamic Loading Using Machine Learning. Dynamics, 4(3), 526-553. https://doi.org/10.3390/dynamics4030028

Article Metrics

Back to TopTop