Download PDF
Research Article  |  Open Access  |  15 Nov 2024

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Views: 244 |  Downloads: 54 |  Cited:  0
J Mater Inf 2024;4:22.
10.20517/jmi.2024.31 |  © The Author(s) 2024.
Author Information
Article Notes
Cite This Article

Abstract

Crystal plasticity (CP) simulation is a powerful tool for studying and understanding the mechanical behavior of materials. A critical aspect of this method is the accurate determination of CP parameters, which ensures that the constitutive model accurately represents the real deformation behavior of a material, especially in high-fidelity simulations. However, identifying these parameters poses a significant challenge due to the high computational cost and the difficulty of finding optimal solutions within a vast and complex parameter space. To address these challenges, we propose a fast search strategy that leverages active learning (AL) and experimental data to accelerate the optimization of CP parameters. Using the Al-Cu eutectic materials as a case study, we introduced a quantitative index, Cecp, to measure the consistency between simulated and experimental stress-strain curves. We demonstrated that Gaussian process regression (GPR) serves as the most appropriate surrogate model for relating CP parameters and Cecp based on our dataset. After only six iterations guided by AL, the optimal CP parameters were successfully identified, resulting in a high-fidelity CP model for analyzing the mechanical behavior of Al-Cu eutectic materials. The machine learning-enhanced strategy is far superior to traditional methods in terms of both efficiency and accuracy. It advances our understanding of the macro-micro relationships of materials and accelerates the material design process.

Keywords

Parameter optimization, crystal plasticity simulation, active learning, stress-strain curve

INTRODUCTION

Materials design fundamentally relies on the comprehensive understanding and prediction of mechanical behaviors. Crystal plasticity (CP) simulation has become a powerful tool for this purpose, offering advantages such as low cost and high fidelity in accurately describing complex deformation behaviors across various scales[1-4].

The deformation behavior of a material can be elucidated by replicating the stress-strain curve based on high-fidelity CP simulation[2,5-8], the accuracy of which hinges on the precise calibration of material parameters. For instance, flow rules and hardening behaviors can be represented using physically or phenomenologically based expressions, including power-law or hyperbolic sine-type expressions[9,10]. These kinematic descriptions involve numerous adjustable material parameters that cannot be directly obtained from experiments. Utilizing imprecise parameters may lead to misleading predictions of mechanical properties from the stress-strain curve, including yield strength, elastic modulus, flow rules, and hardening behaviors, posing significant risks in practical applications. Moreover, the strong nonlinear relationship between CP parameters and the material’s deformation behaviors complicates the rapid optimization of these parameters. Therefore, efficiently identifying the optimal parameters in CP models is an urgent task.

Traditionally, a common approach for CP parameter identification is trial-and-error. Unlike exhaustive strategy, this method relies on expert intervention, and it is simplest to implement[2,6,8,11-13]. However, this random process often requires an unexpected amount of time and computational resources, making it challenging to find a global optimum quickly, especially with a large number of CP parameters. Inverse analysis is another commonly used approach. It utilizes experimental mechanical data to establish a mathematical function linking deformation behavior and CP parameters, which is then used to determine the optimal parameters. Typically, the solution of this function is non-unique, and various optimization methods are applied to find the best one. Guery et al. proposed a cost function based on two least squares criteria to correlate microstructural displacement fields with macroscopic stress, achieving a unique solution using a weighted finite element model[14]. However, the process can become ill-conditioned with even minor input changes. Additionally, it is generally suitable for optimizing 2-3 parameters. To identify more parameters in the CP model, Chakraborty et al. determined four adjustable parameters in the CP constitutive law for the face-centered cubic lattice structure by minimizing the deviation between the measured and simulated indentation response of individual grains[15]. They used a Nelder-Mead simplex algorithm to overcome the influence of crystal orientation on parameter selection. However, the applicability of this method remains unclear for more complex material structures, such as anisotropic crystals. Subsequently, Cauvin et al. identified the CP parameters of a Zn-Cu-Ti alloy, which exhibits highly anisotropic behavior, by an inverse analysis that compared the simulated and experimental results in terms of obtained tensile curves along three different directions[16]. They used the covariance matrix adaptation-evolution strategy genetic algorithm to ensure the solution uniqueness of the inverse problem. Compared to trial-and-error, inverse analysis is indeed more efficient. However, it requires a specific mathematical function for the inverse problem, which is often inadequate for multidimensional parameter selection in simulations of complex structures.

To avoid the randomness in selection and to establish the relationship between multidimensional parameters and the reliability of CP simulation for complex structures, it is necessary to explore a more intelligent approach capable of high-dimensional nonlinear modeling. Active learning (AL), an emerging machine learning paradigm, has great potential for efficiently addressing these issues[17-19]. It has been widely applied in the composition design of complex multi-element alloys, including high-entropy alloys[20,21], high-temperature alloys[22], steels[23,24], and perovskites[25,26]. Its effectiveness in achieving global optimization solutions from vast search spaces has been demonstrated. It is noteworthy that the parameter search in CP models is essentially analogous to the composition design, and the relationship between CP input parameters and the simulation output results of complex structures exhibits high-dimensional, nonlinear characteristics, which is an area where machine learning excels[27-31]. Recently, Goswami et al. and Samaniego et al. employed a neural network to quantitatively analyze the influence of all uncertain input parameters on all uncertain output parameters[32,33]. They also build a unified framework for optimization and inverse analysis and design. Although neural networks are expected to effectively optimize the CP parameters, their performance is highly dependent on the availability of extensive data, and their effectiveness in cases with sparse data sets is somewhat limited. Therefore, identifying a suitable approach to quickly optimize CP parameters remains a challenge in CP modeling.

In this study, we propose a fast parameter optimization strategy for CP simulation leveraging AL and in-situ experiments. Using Al-Cu eutectics as a case study, we establish a quantitative index that reflects the discrepancies between experimental and simulated stress-strain data. Then, we establish the relationship between the index and CP parameters using a high-precision machine learning model. Through several iterations, we identify the optimal CP parameters and examine the consistency between the simulations with these parameters and experimental data. Finally, we analyzed the mechanical behavior of Al-Cu eutectic materials. We believe this strategy proposed in this study can efficiently elucidate the macro-micro relationships of materials and accelerate materials design.

METHODS

Material preparation and characterization

A bulk directional solidification sample of Al-Cu alloys with 17.39 at.% Cu was initially prepared in a vacuum induction melting furnace, employing metal raw materials with a purity of 99.99%. A tensile specimen was then extracted from the bulk sample, and its schematic is depicted in Figure 1A, where the long axis is along the withdrawal direction of directional solidification. Macroscopic stretching of the specimen was conducted along the withdrawal direction of directional solidification (x-axis) with an initial strain rate of 10-4 s-1 using a TSMT EM6.504 tensile testing machine to obtain its macroscopic stress-strain curve, as shown in Figure 1B. Note that the starting point of the engineering stress-strain curve in Figure 1B does not originate from (0, 0). This is due to the fact that the initial force setting of the force sensor used in the testing procedure typically starts recording data when the force exceeds 54 N. Furthermore, Wei et al. and Ebrahimi et al. have also conducted tensile tests on materials using the same experimental equipment, and their stress-strain curves similarly do not originate from (0, 0)[34,35]. Their results proved evidence for the reliability of the experimental results obtained in this study. Additionally, it should be noted that all stress-strain curves presented in this study are of the engineering stress-strain curves. Furthermore, an in-situ tensile specimen was fabricated via electro-discharge machining with a gauge length of 8.0 mm, width of 3.2 mm and thickness of 1.0 mm, and its microstructure was observed using a TESCAN MIRA3 field emission scanning electron microscope (SEM), as shown in Figure 1C. The specimen was subsequently ground and finely polished using colloidal silica. A mono-layer of SiO2 nanoparticles was then deposited on the specimen surface using the drop-casting technique[36] to ensure a dense and uniform speckle morphology, and Figure 1D displays an example of an emission scanning electron microscope and microscopic-digital image correlation (SEM-DIC) speckle image. The microscopic deformation behavior was observed through in-situ tensile testing, with the specimen being deformed at an initial strain rate of 1 × 104 s-1 in a Kammrath & Weiss tensile stage. The parameters of the Kammrath & Weiss system can be found at https://www.kammrath-weiss.com/en/tensile-compression-modules/. After testing, the local strain distribution was analyzed using the GOM Correlate software, and details of the local strain distribution will be analyzed in Section “Sensitivity analysis between CP parameters and simulation output”.

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 1. (A) Schematics of the uniaxial tensile specimen, where the long axis is the directional solidification direction; (B) The engineering stress-strain curve of the tensile specimen; (C) Microstructure in the directional growth direction of the specimen; (D) An example of SEM-DIC speckle image (C). SEM-DIC: Scanning electron microscope and microscopic-digital image correlation.

CP model

To predict the deformation behavior of Al-Cu eutectic materials, the polycrystalline elasto-viscoplastic CP model proposed by Lebensohn et al. was employed. In this study, we introduced the model briefly, and the details were seen in the literature[37-40].

For an elastoplastic issue, the total strain for a system $$ \varepsilon _{ij}^\mathrm{t}(x) $$ can be decomposed into two parts: elastic strain $$ \varepsilon _{ij}^\mathrm{e}(x) $$ and plastic strain $$ \varepsilon _{ij}^\mathrm{p}(x) $$. Using Euler implicit time discretization and Hooke’s law, the total strain at point x with t + Δt can be expressed as

$$ \varepsilon_{i j}^{t, t+\Delta t}(x)=C_{i j k l}^{-1}: \sigma_{k l}(x)+\varepsilon_{i j}^{\mathrm{p}, t}(x)+\dot{\varepsilon}_{i j}^{\mathrm{p}, t+\Delta t}(x, \sigma) \Delta t $$

where $$ C_{ijkl}^{-1} $$ represents the inverse of the elastic tensor Cijkl. The constitutive relation between the plastic strain rate $$ \dot{\varepsilon} _{ij}^\mathrm{p}(x) $$ and the stress σkl is given by the classic incompressible rate-dependent CP equation:

$$ \dot{\varepsilon}_{i j}^{\mathrm{p}}(x)=\sum_{s} m^{s} \dot{\gamma}^{s}(x)=\dot{\gamma}_{0} \sum_{s=1}^{N} m_{i j}^{s}(x)\left(\frac{\left|m_{k l}^{s}(x): \sigma_{k l}(x)\right|}{\tau_{0}^{s}(x)}\right)^{n} \operatorname{sgn}\left(m_{k l}^{s}(x): \sigma_{k l}(x)\right) $$

where $$ \dot{\gamma}^s(x) $$ and $$ \dot{\gamma}_0 $$ are the shear rate and reference shear rate, respectively. $$ m_{ij}^s(x) $$ is symmetric Schmid tensor with slip system s at point x, n is the stress exponent, $$ \tau_0^s(x) $$ is the critical shear stress (CRSS) of the slip system s, and the symbol “:” represents double contraction of indices.

To solve Equation (2) efficiently, the algorithm of fast Fourier transform (FFT) is employed. In the FFT algorithm, the local deformation response of heterogeneous materials can be expressed by convoluted integrals between the Green function of a linear reference medium and a fluctuating polarization field. Thus, using the stiffness of a linear reference medium and a displacement gradient uk,l(x), the stress can be expressed as:

$$ \sigma_{i j}(x)=\sigma_{i j}(x)+C_{i j k l}^{0} u_{k, l}(x)-C_{i j k l}^{0} u_{k, l}(x) $$

$$ \varphi(x)=\sigma_{i j}(x)-C_{i j k l}^{0} u_{k, l}(x)=\sigma_{i j}(x)-C_{i j k l}^{0} \varepsilon_{k l}(x) $$

here, φ(x) is the polarization field, and the local strain field εkl(x) = (uk,l(x) + uk,l(x))/2. When the stress equilibrium is reached, σij,j(x) = 0, one can obtain:

$$ C_{i j k l}^{0} u_{k, l}(x)+\varphi_{i j, j}(x)=0 $$

As the polarization field described by Equation (5) is a function of an unknown strain field, some iterative procedures are required to solve them. To ensure the effectiveness and stability of solving during iteration, the augmented Lagrangian scheme[41] proposed by Michel was applied. In this scheme, two auxiliary fields λ and e were introduced, in which $$ \lambda_{ij}^{(i)} $$ and $$ e_{ij}^{(i)} $$ represent auxiliary guess stress field and strain fields at the i-th iteration, respectively. The polarization field can be rewritten as:

$$ \varphi_{i j}^{(i)}(x)=\lambda_{i j}^{(i)}(x)-C_{i j k l}^{0} e_{k l}^{(i)}(x) $$

Thus we obtain:

$$ e_{i j}^{(i+1)}(x)=E_{i j}+\mathrm{FT}^{-1}\left(\hat{\tilde{e}}_{i j}^{(i)}+\operatorname{sym}\left(\hat{\Gamma}_{i j k l}^{0}(k) \hat{\lambda}_{k l}^{(i+1)}(k)\right)\right. $$

The augmented Lagrangian scheme is suitable for solving Equation (7), which requires the elimination of a residual (R) for each material point x. The residual R is a function of the stress field σi+1 constitutively related to the strain field εi+1:

$$ R_{k}\left(\sigma_{i j}^{(i+1)}(x)\right)=\sigma_{k}^{(i+1)}+C_{k l}^{0} \varepsilon_{l}^{(i+1)}\left(\sigma_{i j}^{(i+1)}(x)\right)-\lambda_{k}^{(i)}-C_{k l}^{0} e_{l}^{(i+1)}=0 $$

Note that, for convenience, in Equation (8), the residual is represented by the contraction symbol for symmetric tensors. For example, σijσk, k = 1, 6; CijmnCkl, k, l = 1, 6, and so on. The nonlinear Equation (8) is solved using a Newton-Raphson (N-R) method, that is:

$$ \sigma_{k}^{(i+1, j+1)}=\sigma_{k}^{(i+1, j)}-\left(\left.\frac{\partial R_{k}}{\partial \sigma_{l}}\right|_{\sigma_{k}^{(i+1, j)}}\right)^{-1} \sigma_{k}^{(i+1, j)} $$

Combing Equations (8) and (1), the Jacobian in Equation (9) can be provided as:

$$ \left.\frac{\partial R_{k}}{\partial \sigma_{l}}\right|_{\sigma_{k}^{(i+1, j)}}=\delta_{k l}+C_{k q}^{0} C_{q l}^{-1}+\left.\Delta t C_{k q}^{0} \frac{\partial \dot{\varepsilon}_{q}^{p}}{\partial \sigma_{l}}\right|_{\sigma_{k}^{(i+1, j)}} $$

The derivative on the right-hand side of Equation (10) is the tangent compliance of Equation (2). Due to the dependence of plastic strain on the stress, the CRSS is a function of stress τ0s(εp(σij(x)) = τ0s(σij(x)). The approximate expression for this tangent compliance is:

$$ \left.\frac{\partial \dot{\varepsilon}_{q}^{\mathrm{p}}}{\partial \sigma_{l}}\right|_{\sigma_{k}^{(i+1, j)}} \cong n \dot{\gamma}_{0} \sum_{s}^{N} \frac{m_{q}^{s} m_{l}^{s}}{\tau_{0}^{s}\left(\sigma_{i j}(x)\right)}\left(\frac{\left|m_{k l}^{s}(x): \sigma_{k l}(x)\right|}{\tau_{0}^{s}\left(\sigma_{i j}(x)\right)}\right)^{n-1} $$

Combining Equations (10) and (11), an approximate Jacobian can be given as[37]:

$$ \left.\frac{\partial R_{k}}{\partial \sigma_{l}}\right|_{\sigma_{k}^{(i+1, j)}} \cong \delta_{k l}+C_{k q}^{0} C_{q l}^{-1}+\left(\Delta t n \dot{\gamma}_{0}\right) \sum_{s}^{N} \frac{m_{q}^{s} m_{l}^{s}}{\tau_{0}^{s}\left(\sigma_{i j}(x)\right)}\left(\frac{\left|m_{k l}^{s}(x): \sigma_{k l}(x)\right|}{\tau_{0}^{s}\left(\sigma_{i j}(x)\right)}\right)^{n-1} $$

Despite this approximation, the convergence rate is almost unaffected even in the case of strain hardening, compared to the case of without hardening.

According to the expression of the Jacobian in Equation (11), different hardening laws are obtained. For the strain hardening behavior of Al-Cu eutectics, the CRSS evolution of each slip system is a function of the total accumulated shear strain[37]:

$$ \tau_{0}^{s}(x)=\tau_{0}^{s *}+\Delta \tau_{0}^{s}=\tau_{0}^{s *}+H \Delta \Gamma=\tau_{0}^{s *}+H \Delta \sum_{s=1}^{N}\left|\dot{\gamma}^{s} \Delta t\right| $$

Where H is a positive scalar related to the strength of the work-hardening behavior. Note that the hardening law employed in this study is phenomenological, and it does not account for dislocation evolution. To encompass more intricate dislocation movements, it is essential to couple the dislocation evolution within the CP model, including dislocation multiplication and annihilation.

Optimization search methodology

Machine learning

Five commonly used machine learning models - linear regression[42], ridge regression[43,44], support vector regression (SVR)[45,46], random forest (RF)[47,48] and Gaussian Process Regression (GPR)[49-51] - were used to establish the relationship between CP parameters and mechanical behavior (stress-strain curves). It is important to note that both SVR and GPR models use an isotropic radial-base kernel function (RBF-kernel). The hyper-parameters of the five models were determined by a method of Bootstrap Sampling Cross-Validation (BSCV), and the prediction uncertainty of the models was estimated by 1,000 samples.

The performance of these models was quantified by the root mean square error (RMSE), defined as

$$ \text { RMSE }=\sqrt{\frac{1}{n} \sum_{i=1}^{n}\left(\mu_{i}-\hat{\mu}_{i}\right)^{2}} $$

where μi and $$ \hat{\mu}_i $$ are the experimental and predicted values at the i-th strain position, respectively. n is the number of discrete points. The smaller the RMSE, the closer the predictions are to experiments.

To accurately elucidate the relationship between CP parameters and a quantitative indicator Cecp, five candidate machine learning models - Linear, Ridge, SVR, RF and GPR - were trained, and the best model was subsequently selected. Here, Cecp is used to evaluate the quality of the CP parameter, which lies in characterizing the discrepancy between experimental and simulated stress-strain curves. Further details will be provided in Section “Quantitative evaluation of CP parameter quality”. The hyper-parameter of SVR and RF was determined through exhaustive search and BSCV. The regularization coefficient (C) and the RBF kernel’s length related parameter (γ) were identified as the principal parameters for the SVR model; the number of trees (N) and the maximum tree depth (H) were designated as key parameters for the RF model; the scale parameter(Length_scale) was primarily emphasized within the GPR model; the determination of these hyperparameters is shown in Supplementary Figure 1. The optimal hyperparameters for the SVR, RF and GPR models are listed in Table 1. Figure 2 presents a comparison between the actual and predicted Cecp by the five models using the initial training set. A significantly larger number of data points deviate from the diagonal line in Figure 2A and B compared to Figure 2C-E. This observation suggests that the Linear and Ridge models fail to capture complex relationships between CP parameters and Cecp, unlike SVR, RF, and GPR models. Furthermore, nearly all sample points in the GPR model are distributed near the diagonal, underscoring the superior performance of the GPR model relative to other models. Additionally, the RMSE comparison among the five models, depicted in Figure 2F, reveals that the GPR model achieves the lowest RMSE, indicating its highest accuracy in predicting Cecp. Consequently, the GPR model was chosen to construct the AL framework.

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 2. Accuracy evaluation of the trained ML models. (A) Linear regression; (B) Ridge regression; (C) SVR; (D) RF; (E) GPR; (F) RMSE comparison among the models. ML: Machine learning; SVR: support vector regression; RF: random forest; GPR: Gaussian Process Regression; RMSE: root mean square error.

Table 1

Optimal hyperparameters for SVR, RF and GPR models

ModelsHyperparameters
SVRC = 45, γ = 1.45
RFN = 95,H = 7
GPRLength_scale = 1.85

Search strategy

According to Equations (2) and (13), most CP parameters can be confirmed from Refs[2,52-55]. Only four parameters remain to be undetermined: $$ \dot{\gamma}_{\mathrm{0\_{Al}}} $$, $$ \dot{\gamma}_{\mathrm{0\_{Al_2Cu}}} $$, HAl and $$ H_{\mathrm{Al_2Cu}} $$. Note that these determined CP parameters come from precise experimental measurements, and can be directly used in CP simulations without the need for calibration. Optimization becomes necessary only for the parameters that are either unavailable from experiments or cannot be experimentally measured, that is $$ \dot{\gamma}_{\mathrm{0\_{Al}}} $$, $$ \dot{\gamma}_{\mathrm{0\_{Al_2Cu}}} $$, HAl and $$ H_{\mathrm{Al_2Cu}} $$. Among them, $$ \dot{\gamma}_{\mathrm{0\_{Al}}} $$ and $$ \dot{\gamma}_{\mathrm{0\_{Al_2Cu}}} $$ are related to the dislocation movements of α-Al and θ-Al2Cu, respectively[37,56,57], while HAl and $$ H_{\mathrm{Al_2Cu}} $$ are related to the work-hardening behavior of α-Al and θ-Al2Cu, respectively[37,57]. Generally, the range of $$ \dot{\gamma}_0 $$ is from 10-2 to 10-5, with larger values indicating easier dislocation movement. The range of H is from 0 to 1,000 MPa, with larger values indicating stronger work-hardening behavior[37,57]. If 10-5 and 10 MPa are the steps of $$ \dot{\gamma}_0 $$ and H within their respective ranges, then they have 103 and 102 possible values, respectively. Thus, the total number of parameter combinations is 103 × 103 × 102 × 102 = 1010. Undoubtedly, searching for the optimal combination in such a large space (denoted as Ω in this study) is a very challenging task. Therefore, we constructed an adaptive framework for CP parameter optimization by borrowing the concept of AL from the field of alloy composition design. The complete process of finding the optimal CP parameters in this study is shown in Figure 3.

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 3. Strategy of CP parameter optimization based on AL. CP: Crystal plasticity; AL: active learning.

Part B in Figure 3 illustrates an intelligent framework called AL to guide the iterative refinement of CP parameters. The main process is as follows:

(1) A total of 26 groups of initial parameters were selected in the search space Ω. Using these parameters, 26 simulated stress-strain curves were generated using the CP model, and the corresponding Cecp values were calculated.

(2) If the maximum Cecp (denoted as $$ C_{\mathrm{ecp}}^{\mathrm{max}} $$) exceeds 0.99, the corresponding inputs are the optimal CP parameters. If not, proceed to the next step.

(3) An appropriate machine learning model f was trained using the initial parameters as inputs and the corresponding 26 Cecp data as outputs. The model f was then used to predict Cecp for all parameter combinations in Ω, with the mean and variance of each prediction denoted as μ and σ, respectively.

(4) The efficient global optimization (EGO) utility function was applied to select three parameter combinations by maximizing expected improvement (EI) across the entire Ω. EI is defined as EI = σ[ϕ(z) + (z)], where z = (μ - μ*)/σ, μ* is $$ C_{\mathrm{ecp}}^{\mathrm{max}} $$, and ϕ(z) and ψ(z) denote the standard normal density and distribution functions, respectively.

(5) The three selected parameter groups were then used as input to the CP model, and their corresponding Cecp values were calculated.

(6) The initial set was updated with the new recommendations. Steps 2-5 were then repeated until $$ C_{\mathrm{ecp}}^{\mathrm{max}} $$ > 0.99, resulting in the identification of the optimal parameter combinations.

It should be noted that the AL strategy above is an adaptive loop that iteratively aligns the CP simulation with experimental results by maximizing Cecp. Once the optimal parameters are identified, the resulting CP model will have high fidelity. The use of machine learning in this process makes the parameter search intelligent and efficient.

Note that the mechanical properties and microstructure of the Al-Cu eutectic alloy at room temperature were displayed in Figure 1. In the CP modeling, the experimental microstructure serves as the input of the CP model, with comprehensive details provided in Ref[52]. Additionally, the initial configurations and material parameters of the CP simulation for the Al-Cu eutectic materials are from Ref[52].

Quantitative evaluation of CP parameter quality

The key to evaluating the quality of the CP parameter lies in characterizing the discrepancy between experimental and simulated stress-strain curves. To quantify this discrepancy, we discretized both the experimental and simulation curves into 17 data points, ensuring that the horizontal coordinates of each set of points were consistent. Then, the determination coefficient R2 was used to assess the consistency between discrete experimental and CP simulation data. Notably, to distinguish it from the accuracy evaluation of machine learning illustrated in Figure 2, we redefine it as Cecp, that is:

$$ C_{\mathrm{ecp}}=R^{2}=1-\frac{\sum_{i=1}^{n}\left(\mu_{i}-\hat{\mu}_{i}\right)^{2}}{\sum_{i=1}^{n}\left(\mu_{i}-\bar{\mu}\right)^{2}} $$

where $$ \bar{\mu}_i $$ represents the average of n discrete points. The closer Cecp is to 1, the closer the simulations are to the experiments.

To validate the effectiveness of Cecp in assessing CP parameter quality, Figure 4 shows the degree of alignment between two examples of simulated stress-strain curves (simulation 1 and simulation 2) and the experimental curve. The CP parameters of the two examples are listed in Table 2. Observations from Figure 4 revealed that, in the nearly linear portion of the stress-strain curve, the data of Simulations 1 and 2 are in good agreement with the experimental data, demonstrating that the elastic parameters of the simulations align with those of experiments. This concurrence is attributed to the fact that the elastic parameters utilized in this study come from experimental measurements. Upon departing from the initial linearity, the data points from Simulation 2 (Cecp = 0.65) significantly deviate more from the experimental data than those from Simulation 1 (Cecp = 0.86). This discrepancy implies that the dislocation movements (reflected by parameters $$ \dot{\gamma}_{0\_\mathrm{Al}} $$ and $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$) and the work-hardening behavior (reflected by parameters $$ \dot{\gamma}_{0\_\mathrm{Al}} $$ and $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$) in Simulation 2 are more deviated from the experiment than that in Simulation 1. Additionally, the curve of simulation 1 is closer to the experimental curve than that of simulation 2. The Cecp value between simulation 1 and experiment is 0.86, which is higher than that for simulation 2 (0.65), as shown in Table 2. This indicates that a higher Cecp value corresponds to a smaller disparity between simulation and experiment. This further demonstrates that Cecp is indeed a valid metric for quantifying the quality of CP parameters. In this study, 0.99 is selected as a threshold of Cecp.

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 4. Example of quantitative evaluation of CP parameters. CP: Crystal plasticity.

Table 2

Qualitative and quantitative assessment of the degree of agreement between simulation and experiment under two sets of parameters

$$ \dot{\gamma}_{0\_\mathrm{Al}} $$$$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$HAl (MPa)$$ H_{\mathrm{Al_2Cu}} $$ (MPa)The degree of agreementCecp
Simulation 10.010.0000110100Proximity0.86
Simulation 20.0070.0035400900Deviation0.65

RESULTS AND DISCUSSION

Iterative optimization of CP model parameter

At the early stage, a trial-and-error approach was performed to establish an initial training set, as detailed in the Supplementary Materials. The maximum Cecp in the initial training set is 0.95, indicating that the simulated stress-strain curve deviates far from the experimental curve. Following the AL framework shown in Figure 2, we developed a high-precision GPR model, which was utilized to iteratively improve the CP parameters. It is important to note that when using a GPR-based search strategy to recommend new parameter combinations, the constraint must be met that the parameter $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$ must be less than the parameter $$ \dot{\gamma}_{0\_\mathrm{Al}} $$. This is because, during the deformation of Al-Cu eutectic alloys at room temperatures, the activation of dislocations in the θ-Al2Cu phase is more difficult than that in the α-Al phase.

As illustrated in Figure 5, Cecp surpassed 0.99 in just four iterations, and a total of six iterations were executed. Detailed recommendations are provided in Table 3. In the first iteration, there is a significant difference among the three points, and the maximum Cecp is only 0.86, which is below the threshold. In the second iteration, the Cecp values noticeably increase compared to those from the first iteration, and the data range significantly reduces. This indicates that the GPR samples are centrally located at a specific position of data sparsity in Ω. Notedly, Cecp continues to rise until it reaches 0.99 in the fourth iteration. Interestingly, in this iteration, the three sets of CP parameters are very close to each other, suggesting that the global optimal may be successfully identified. As there is no more suitable data in Ω to further improve Cecp, its values start to decrease from the fifth iteration onward. By the sixth iteration, the data points began to disperse again, indicating that the exploration in Ω guided by GPR has resumed to improve the data distribution of the training set. Furthermore, the Cecp values of the recommended CP parameters in the 5th and 6th iterations are consistently lower than the Cecp value obtained after the 4th iteration, indicating that the optimal CP parameters have been identified. Notedly, Cecp values in the 5th and 6th iterations have been over the preset threshold of 0.99 but the iteration continues, which is just to make it more intuitive that the CP parameters in the 4th iteration are optimal. Additionally, as shown in Table 3, among the set of optimal CP parameters that correspond to the maximum Cecp value, the parameter HAl is approximately 100 MPa and $$ H_{\mathrm{Al_2Cu}} $$ is around 300 MPa. This consistency with the sensitivity analysis presented in Section “Comparative analysis between the AL strategy and the typical traditional methods” further confirms the discovery of the optimal CP parameters. Since the optimal has been discovered, the iterative search is terminated.

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 5. Cecp for the recommendations in a total of six iterations. Iteration 0 corresponds to the initial training data.

Table 3

Recommended CP parameters and Cecp by AL

Iteration$$ \gamma_{0\_\mathrm{Al}} $$$$ \gamma_{0\_\mathrm{Al_2Cu}} $$HAl (MPa)$$ H_{\mathrm{Al_2Cu}} $$ (MPa)Cecp
10.00010.0007640200.51
0.00030.0003640200.81
0.00040.0002640200.86
20.010.0003910200.91
0.010.0002910200.91
0.010.0001910200.91
30.00050.00012004000.96
0.0050.0005100010000.98
0.0030.00054005000.96
40.0010.00053005000.99
0.0010.00051002000.99
0.00150.00053004000.99
50.010.00485502900.84
0.010.00495502900.84
0.010.0055502900.84
60.00050.00012507100.95
0.00010.00092507100.59
0.00030.00022507100.88

For an AL process that maximizes EI, the initial sparser training data leads to greater prediction uncertainty of GPR on the new sample points and increases the randomness of data distribution[58], as shown in the first iteration in Figure 5. The prediction accuracy on each recommendation is depicted in Figure 6A. Evidently, GPR overestimates the Cecp for the three new sets of CP parameters, resulting in significant prediction errors and uncertainties. From the mathematical expression of the utility function EI that is used to recommend the next experimental points, it is easy to understand that given a large uncertainty σ, EI → σ. That is, in the limit of high uncertainty, EI is maximized by choosing the next best material with the largest uncertainty (“exploration”). It represents a region in our search space, where the regressor is most uncertain. Thus, for regions not covered by the training data, GPR will also centrally sample in the early stages. Consequently, the values of Cecp at the second iteration are almost consistent, but they deviate significantly from the predictions, as indicated by the orange points in Figure 6A. Once the training set’s sparse area is filled, prediction uncertainty decreases, model accuracy increases, and the search moves closer to the optimal solution[58]. On this occasion, EI → μ - μ*, where μ is the predicted average by the machine learning (ML) model, μ* is the current maximum value in the training set. This means that in the limit of no uncertainty or if uncertainties are almost the same for all points, EI is maximized by choosing the largest μ (i.e., “exploitation”). Therefore, the next best material for measurement is chosen by faithfully trusting the predictions from the regressor. This conclusion is supported by the trends of the green and red dots (at the fourth and fifth iterations) in Figure 6A. After identifying the global optimization, GPR explores other sparse regions in Ω, leading to a more dispersed distribution of new sampling points, yet with accurate predictions, as shown by the purple and brown spots in Figure 6A. Additionally, Figure 6B shows that EI rapidly decreases from the first to the fifth iteration. Once the maximum Cecp has been identified, the probability of new recommendations improving upon the best training data is reduced to 0, causing EI to drop to 0.

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 6. Prediction accuracy and EI variation during the optimization process of CP parameters. (A) Accuracy between actual and predicted Cecp; (B) Variation of EI with the number of iterations. EI: Expected improvement; CP: crystal plasticity.

Mechanical behavior analysis of Al-Cu eutectic materials

To validate the optimal CP parameters, we next input the parameter into the CP model, obtained a simulated stress-strain curve, and compared it with the experimental curve. It is crucial to emphasize that, in this study, the input microstructure for the CP model is firmly based on direct experimental observations[52]. For quantitative comparisons between the experimental data and the CP simulations, the horizontal axes of the experimental and simulated stress-strain curves have been aligned. For comparison, other curves for the best parameter combinations at each AL iteration are also included, as shown in Figure 7. Note that the CP simulation of a complex stress status with the xy shear deformation has been provided in Supplementary Figure 2 to ensure the validation of the CP simulation. From Figure 7, it is evident that the simulated stress-strain curve for the second iteration (orange curve) is closer to the experiment (black curve) than that for the first iteration (blue curve). This is also reflected in the comparison of the Cecp list in Table 3. Similar patterns can be observed in other curves, demonstrating that Cecp is a suitable quantitative index for evaluating the reliability of CP simulation. Notably, the red curve in Figure 7 nearly coincides with the black dots, indicating that the recommendations from the fourth iteration may be optimal. In addition, the deviation of the purple (fifth iteration) and brown (fifth and sixth iterations) curves in Figure 7 from the black dots indicates that the recommended parameters begin to deviate from the optimal solution after the fourth iteration. This further demonstrates the recommendations from the fourth iteration are indeed optimal, resulting in a high-fidelity simulation that closely matches reality. Note that the elastic part of CP simulated stress-strain curves with different iterations is consistent due to the fact that the elastic parameters in this study are sourced from literature[55,59,60].

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 7. Comparison of simulated and experimental stress-strain curves. The simulations are performed with the best parameter combinations at each AL iteration. AL: Active learning.

Figure 8 shows a comparison between the simulated local strain distributions and the results from in-situ DIC experiments. The CP simulation was performed using the optimal parameters. As shown in Figure 8A(i)-(iii), with rising total strain εt, experimental strain concentration at microstructural defects increases, marked by the black dashed lines. Similarly, in the CP simulations, as shown in Figure 8B(i)-(iii), an analogous trend in strain concentration at microstructural defects with increasing total strain is observed, also marked by the black dashed lines, demonstrating the consistency of the experiments and simulations in microscopic deformation. This further illustrates the reliability of the optimal parameters calibrated using our proposed search strategy. In addition, when the total strain is 0.001 or 0.0025, the minimum strain is not equal to 0 in CP simulation, as shown in Figure 8B(ii) and (iii), which is different from the experimental observation shown in Figure 8A(i) and (iii). The discrepancy arises from the fact that the CP model does not account for dislocation reactions between different slip systems, such as dislocation annihilation and multiplication. In the actual room-temperature tensile process of the Al-Cu eutectic material, there exist dislocation reactions (dislocation multiplication or annihilation) in the α-Al phase, and more significantly, these mobile dislocations can cross the two-phase interface and enter the Al2Cu phase, further affecting the local strain within the Al2Cu phase[52,61-63]. However, since our CP model ignored the dislocation multiplication or annihilation, there is a discrepancy between the simulation and experimental results as depicted in Figure 8. It is worth noting that although our CP model primarily investigates the situation where the Al2Cu phase is without plastic deformation, it remains applicable to cases where considerable plastic deformation is present in the Al2Cu phase. For further details, refer to Supplementary Figure 3.

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 8. Comparison between the simulated and experimental local strain distributions. The CP simulation was completed with the optimal parameters. CP: Crystal plasticity.

To further explain microscopic deformation behavior, we examined the stress and plastic strain distributions at different total strains, as shown in Figure 9. In Figure 9A(i)-(iii), the red phase and blue phase represent α-Al phase and θ-Al2Cu phase, respectively. Figure 9B(i)-(iii) shows that the minimum stress ranges from 4 to 8 MPa, while the maximum stress increases from 8.5 to 46 MPa as the total strain εt increases from 0.0005 to 0.0025. In Figure 9C(i)-(iii), the maximum plastic strain gradually increases, whereas the minimum plastic phase remains at 0. The minimum stress and the maximum plastic strain are observed in the α-Al phase, while the maximum stress and minimum plastic are found in the θ-Al2Cu phase. This indicates that plastic deformation occurs in the α-Al phase even under small stresses, whereas the θ-Al2Cu phase is subjected to large stresses without plastic deformation.

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 9. Phase field, stress, and plastic strain distributions with different total strains in CP simulations. (A) (i-iii) Phase field; (B) (i-iii) Stress distributions; (C) (i-iii) Plastic strain distributions. CP: Crystal plasticity.

To further support this conclusion, we analyzed linear stress and plastic strain distributions along the x-axis [labeled by the black dotted line in Figure 9B(i)], and the results are shown in Figure 10. It was found that as the total strain εt increases, the stresses in both the α-Al and θ-Al2Cu phases increase, with the stress in the θ-Al2Cu phase being significantly higher than that in the α-Al phase. As the total strain increases, the increased rate of the stress difference between σmax and σmin in the α-Al phase is significantly slower than that in the θ-Al2Cu phase, as shown in Table 4. This indicates that the stress variation in the θ-Al2Cu phase is much greater than that in the α-Al phase. Moreover, only the plastic strain in the α-Al phase shows a gradual increase, while plastic deformation in the θ-Al2Cu phase remains at 0. For dual-phase materials, the deformation behavior is not only related to the stress and strain distributions within each phase but also significantly influenced by the inhomogeneous deformation near the interfaces between the two phases. In this study, we quantified the non-uniform deformation near the α-Al and θ-Al2Cu phase interface under different total strains. To better investigate the inhomogeneous deformation near the phase interface, the maximum plastic strain difference of the α-Al and θ-Al2Cu phases near the interface with different total strains was counted. Since the plastic deformation of the θ-Al2Cu phase is zero, the inhomogeneous deformation of the α-Al and θ-Al2Cu phase dual-phase material is primarily affected by the α-Al phase. The statistical results show that under total strains of 0.0005, 0.001, and 0.0025, the maximum plastic strain in the α-Al phase is 0.000018, 0.000375, and 0.00317, respectively. This indicates that as the total strain grows, the inhomogeneous deformation of the Al-Cu eutectic alloy significantly increases. These findings are not only consistent with the DIC experiment of Figure 8 but also align with the previous literature reports[52,62-64], further confirming the reliability of our search strategy in optimizing CP parameters.

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 10. Line stress and plastic strain distributions with different total strains. (A) Line stress; (B) Plastic strain distributions. The location of the lines is labeled by the black dotted line in Figure 9B(i).

Table 4

The stress and increased rate of the stress difference in the α-Al and θ-Al2Cu phases as the total strain increases

εtσmin (MPa)σmax (MPa)Δσ (MPa)Increase rate of σ
αAl$$ \mathrm{\theta_{Al_2Cu}} $$αAl$$ \mathrm{\theta_{Al_2Cu}} $$αAl$$ \mathrm{\theta_{Al_2Cu}} $$αAl$$ \mathrm{\theta_{Al_2Cu}} $$
0.0005405455651511--
0.001749095125213540%220%
0.00259614315831862160200%360%

In summary, the AL-based search strategy successfully optimized the CP parameters for the Al-Cu eutectic alloy at room temperature. Although the search strategy in this study was only applied to the parameter optimization of a quite simple CP model with only four parameters, it can be anticipated that the strategy proposed in this study is also suitable for more complex systems, such as torsional deformation[65,66]. In this case, similar treatment for the parameters and the quantification of the simulation reliability should be performed, and more iterations may be required due to more undetermined parameters than those in our study, but even so, it is faster than the trial-and-error selection. However, we also acknowledge that for the parameter optimization of non-uniform computations[67,68], our strategy may need to be further improved.

Sensitivity analysis between CP parameters and simulation output

To quantify the impact of the CP input parameters ($$ \dot{\gamma}_{0\_\mathrm{Al}} $$, $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$, HAl and $$ H_{\mathrm{Al_2Cu}} $$) on the CP simulation output Cecp, we assumed four cases by a control variable method. The predicted results by GPR are given in Figure 11. Note that four cases in this study were determined randomly yet all meet the restriction that $$ \dot{\gamma}_{0\_\mathrm{Al}} $$ is greater than $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$ due to the slip resistance of α-Al phase being significantly lower than that of θ-Al2Cu phase at room temperature[2,53].

Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning

Figure 11. Sensitivity analysis between CP parameters and Cecp. (A) $$ \dot{\gamma}_{0\_\mathrm{Al}} $$vs. Cecp with a constant of $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$, HAl and $$ H_{\mathrm{Al_2Cu}} $$; (B) $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$vs. Cecp with a constant of $$ \dot{\gamma}_{0\_\mathrm{Al}} $$, HAl and $$ H_{\mathrm{Al_2Cu}} $$; (C) HAlvs. Cecp with a constant of $$ \dot{\gamma}_{0\_\mathrm{Al}} $$, $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$ and $$ H_{\mathrm{Al_2Cu}} $$; (D) $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$vs. Cecp with a constant of $$ \dot{\gamma}_{0\_\mathrm{Al}} $$, $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$ and HAl. CP: Crystal plasticity.

Figure 11A displays the variation of Cecp with respect to $$ \dot{\gamma}_{0\_\mathrm{Al}} $$ while $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$ (0.00001), HAl (500 MPa) and $$ H_{\mathrm{Al_2Cu}} $$ (500 MPa) are held constant. It is evident that as $$ \dot{\gamma}_{0\_\mathrm{Al}} $$ increases, Cecp decreases linearly. This is due to the fact that the dislocation movement in the α-Al phase is increasingly arduous with the increasing parameter $$ \dot{\gamma}_{0\_\mathrm{Al}} $$, resulting in elevating the yield strength of the α-Al phase. Considering that the yield strength of the α-Al phase is inferior to that of the Al-Cu eutectics, a reduction in $$ \dot{\gamma}_{0\_\mathrm{Al}} $$ brings the CP simulated outcomes closer to the experimental data, that is a larger Cecp. Similarly, as the parameter $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$ diminishes, the dislocation movement in the θ-Al2Cu phase is increasingly hindered, rendering CP simulation outcomes more congruent with the experimental findings (a greater Cecp), as shown in Figure 11B. Herein, $$ \dot{\gamma}_{0\_\mathrm{Al}} $$ (0.01), HAl (500 MPa) and $$ H_{\mathrm{Al_2Cu}} $$ (500 MPa) are maintained as fixed values.

Figure 11C displays the variation of Cecp with respect to HAl while $$ \dot{\gamma}_{0\_\mathrm{Al}} $$ (0.01), $$ \dot{\gamma}_{0\_\mathrm{Al_2Cu}} $$ (0.001) and $$ H_{\mathrm{Al_2Cu}} $$ (500 MPa) are held constant. As HAl increases, an initial rise is followed by a decline in Cecp, and the peak value appears when $$ \dot{\gamma}_{0\_\mathrm{Al}} $$ is near 100 MPa. HAl is related to the hardening behavior of phase α-Al, and the hardening behavior of the material remains essentially unchanged under a given deformation condition (room temperature). Therefore, with the variation of the HAl parameter, an optimal point of HAl exists at which Cecp attains its maximum value. Similarly, regarding the $$ H_{\mathrm{Al_2Cu}} $$ parameter, its correlation with the hardening behavior of phase θ-Al2Cu is noted, and a comparable trend between $$ H_{\mathrm{Al_2Cu}} $$ and Cecp is observed, as shown in Figure 11D. Notably, the peak Cecp value is observed when $$ H_{\mathrm{Al_2Cu}} $$ is around 300 MPa.

From the analysis above, it can be concluded that the AL-based CP parameter optimization strategy has effectively quantified the influence of CP input parameters on simulation reliability for Al-Cu eutectic alloys. This strategy can be easily extended to the optimization of other CP parameters, in which the parameters required to be optimized must lie within a certain range, and a quantitative indicator should be employed to evaluate the deviation between CP simulation results and the experiments. For example, as reported in Ref[2], the parameters $$ g_0^s $$, h, $$ \dot{\gamma}_{0} $$ and m can also be determined by our proposed strategy. Employing these parameters as inputs of the strategy and the discrepancies between simulations and experiments as the output, one can rapidly determine the most optimal group of these parameters in several iterations.

Comparative analysis between the AL strategy and the typical traditional methods

The traditional methods used for optimizing CP parameters mainly include trial-and-error and inverse analysis. The former is the most common and the simplest method. It involves randomly selecting an initial set of parameters to construct a CP model and evaluating the discrepancy between CP simulations and experimental results. This process is repeated until either an acceptable discrepancy is achieved or computing resources are excessively consumed[2,6,8,11-13]. Notably, the discrepancy is often qualitatively judged by visual observation of stress-strain curves or other patterns. Additionally, the arbitrary nature of determining when to stop CP computation can lead to excessive cost consumption[69,70]. Moreover, the results of the trial-and-error method depend on the operator’s experience and intuition. When a large number of parameters need to be determined, there is a significant risk of not quickly finding the optimal solution. Reviewing the AL strategy proposed in Figure 2, we observe a distinct advantage: it quantifies the discrepancy between simulated and experimental stress-strain curves by a simple index, Cecp, whose reliability has been demonstrated in Figures 3 and 7. Notedly, Cecp is mathematically identical to R2. The renaming to Cecp is to distinguish it from R2 for the performance evaluation of the machine learning algorithms used in this study. The usage of Cecp enhances the comprehensibility of our methodology and simultaneously ensures the ease of data handling. Moreover, Cecp provides a quantitative bridge that guides the simulation process from experimental data. Furthermore, by using our strategy, as shown in Figures 5 and 6, and Table 3, intelligent searching of high-dimensional parameters has been realized, which significantly reduces the computational overhead in terms of time and cost, requiring only six iterations.

The inverse analysis method begins by proposing a set of initial parameters, and then evaluates the deviation between CP simulation and experimental results within the trust region by a trust-region policy optimization algorithm[71,72]. If the deviation is large, the algorithm adjusts the trust region, repeating the process until the deviation is minimized. The inverse analysis allows for a direct connection between parameters and simulated results (e.g., mechanical behavior of materials), but it can lead to unnecessary resource wastage[73]. Additionally, the mathematical modeling of the algorithm in the trust region is highly complex, and the results significantly depend on the rationality of the initial guessed parameters[73]. Generally, optimization algorithms are also needed to ensure the uniqueness of the solution[73]. Unlike the inverse analysis, our proposed strategy achieves optimization by maximizing the EI of candidate parameter combinations. This approach balances “exploration” (improving data distribution) and “exploitation” (approaching the optimal) within only a few iterations[20,23,25,26]. The search results are unaffected by the initial dataset as the model input is updated from each iteration. Although we use a non-parametric model, GPR, to ensure accuracy and efficiency, the choice of the surrogate models in our strategy is flexible, including other ML algorithms and genetic algorithms.

We acknowledge that this search strategy has limitations and may only be suitable for some specific scenarios, such as when the material’s deformation behavior is adequately represented by stress-strain curves or other one-dimensional discrete data. Note that the CP parameters in this study were optimized by a single stress-strain curve from an experimental measurement. However, for higher accuracy, stress-strain curves should be based on multiple replicate experiments. So, to better optimize the CP parameters, it is necessary to optimize them using data from multiple replicate experiments in future work.

CONCLUSIONS

A search strategy for CP parameter optimization has been proposed with the assistance of AL and in-situ experiments. The strategy has been applied to an Al-Cu eutectic material, with its mechanical behavior simulated using the CP model with optimal parameters. The main conclusions are summarized as follows:

(1) A quantitative indicator, Cecp (renaming of R2), is used to describe the differences between CP simulations and experiments. It has a simple mathematical form and can quantify simulation reliability.

(2) An adaptive CP parameter optimization strategy based on AL and in-situ experiment was proposed. This strategy finds the optimal CP parameters through CP simulations of only 44 parameter combinations and six iterations of learning for Al-Cu eutectics, which greatly reduces the computational overhead in time and economy, and the selection of the surrogate models used in the strategy is flexible enough.

(3) GPR is the most appropriate model for correlating CP parameters and Cecp. The optimal CP parameters are successfully identified through only six AL iterations. Validation with experiments shows that Cecp for the optimal parameters exceeds 0.99. Using these parameters, a high-fidelity CP model is established.

(4) The deformation behavior of Al-Cu eutectics at room temperature was revealed. It is shown that the α-Al phase is subjected to small stress with large plastic deformation, while the θ-Al2Cu phase is subjected to high stress without plastic deformation.

The developed AL-based search strategy for CP parameter optimization in this study was solely employed to determine four parameters for the Al-Cu eutectics at room temperature. Obviously, this strategy can be easily extended to optimize CP parameters for different materials and deformation conditions. Using this strategy, we have efficiently calibrated the CP parameters of the Al-Cu eutectic material at various deformation temperatures, as shown in Supplementary Figure 4, and the associated research findings are pending publication.

DECLARATIONS

Authors’ contributions

Conceptualization: Li J, Wang Z, Wang J

Methodology, validation, formal analysis and data curation: Jiang M, Hu X, Xing C, Yang Z, Chen Y

Investigation, write-review and editing and visualization: Wang J, Jiang M, Hu X

Availability of data and materials

The data supporting the findings of this study are included in the Supplementary Materials.

Financial support and sponsorship

This work was supported by the National Key Research and Development Program (Grants: 2023YFC3010902). The authors also thank the High-Performance Computing Center of Northwestern Polytechnical University, China, for its computer time and facilities.

Conflicts of interest

Hu X is affiliated with Xi’an Rare Metal Materials Institute Co. Ltd., while the other authors have declared that they have no conflicts of interest.

Ethical approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Copyright

© The Author(s) 2024.

Supplementary Materials

REFERENCES

1. Zhang L, Li J, Wang Y, et al. Understanding the formation of equiaxed α during dynamic precipitation in titanium alloys by elastoplastic phase field simulation. J Mater Res Technol 2023;27:8181-96.

2. Liu G, Xie D, Wang S, Misra A, Wang J. Mesoscale crystal plasticity modeling of nanoscale Al–Al2Cu eutectic alloy. Int J Plasticity 2019;121:134-52.

3. Zhang J, Li H, Sun X, Zhan M. A multi-scale MCCPFEM framework: modeling of thermal interface grooving and deformation anisotropy of titanium alloy with lamellar colony. Int J Plasticity 2020;135:102804.

4. Zhang H, Liu J, Sui D, Cui Z, Fu M. Study of microstructural grain and geometric size effects on plastic heterogeneities at grain-level by using crystal plasticity modeling with high-fidelity representative microstructures. Int J Plasticity 2018;100:69-89.

5. Liu P, Wang Z, Xiao Y, et al. Integration of phase-field model and crystal plasticity for the prediction of process-structure-property relation of additively manufactured metallic materials. Int J Plasticity 2020;128:102670.

6. Kotha S, Ozturk D, Ghosh S. Parametrically homogenized constitutive models (PHCMs) from micromechanical crystal plasticity FE simulations: Part II: thermo-elasto-plastic model with experimental validation for titanium alloys. Int J Plasticity 2019;120:320-39.

7. Nguyen C, Cailletaud G, Barbe F, Marini B, Nguyen D, Phan H. Identification of crystal plasticity parameters for a non-irradiated and irradiated A508 bainite steel. Metall Res Technol 2021;118:204.

8. Li J, Romero I, Segurado J. Development of a thermo-mechanically coupled crystal plasticity modeling framework: application to polycrystalline homogenization. Int J Plasticity 2019;119:313-30.

9. Bandyopadhyay R, Prithivirajan V, Sangid MD. Uncertainty quantification in the mechanical response of crystal plasticity simulations. JOM 2019;71:2612-24.

10. Sedighiani K, Diehl M, Traka K, Roters F, Sietsma J, Raabe D. An efficient and robust approach to determine material parameters of crystal plasticity constitutive laws from macro-scale stress–strain curves. Int J Plasticity 2020;134:102779.

11. Zhang K, Holmedal B, Hopperstad O, et al. Multi-level modelling of mechanical anisotropy of commercial pure aluminium plate: crystal plasticity models, advanced yield functions and parameter identification. Int J Plasticity 2015;66:3-30.

12. Tu X, Shahba A, Shen J, Ghosh S. Microstructure and property based statistically equivalent RVEs for polycrystalline-polyphase aluminum alloys. Int J Plasticity 2019;115:268-92.

13. Azhari F, Wallbrink C, Sterjovski Z, et al. Predicting the complete tensile properties of additively manufactured Ti-6Al-4V by integrating three-dimensional microstructure statistics with a crystal plasticity model. Int J Plasticity 2022;148:103127.

14. Guery A, Hild F, Latourte F, Roux S. Identification of crystal plasticity parameters using DIC measurements and weighted FEMU. Mech Mater 2016;100:55-71.

15. Chakraborty A, Eisenlohr P. Evaluation of an inverse methodology for estimating constitutive parameters in face-centered cubic materials from single crystal indentations. Eur J Mech A Solid 2017;66:114-24.

16. Cauvin L, Raghavan B, Bouvier S, Wang X, Meraghni F. Multi-scale investigation of highly anisotropic zinc alloys using crystal plasticity and inverse analysis. Mater Sci Eng A 2018;729:106-18.

17. Cao B, Yang S, Sun A, Dong Z, Zhang T. Domain knowledge-guided interpretive machine learning: formula discovery for the oxidation behavior of ferritic-martensitic steels in supercritical water. J Mater Inf 2022;2:4.

18. Zhang T. New tool in the box. J Mater Inf 2021;1:1.

19. Lu T, Li M, Lu W, Zhang T. Recent progress in the data-driven discovery of novel photovoltaic materials. J Mater Inf 2022;2:7.

20. Guo C, Hu X, Han X, et al. Laser precise synthesis of oxidation-free high-entropy alloy nanoparticle libraries. J Am Chem Soc 2024;146:18407-17.

21. Debnath A, Krajewski AM, Sun H, et al. Generative deep learning as a tool for inverse design of high entropy refractory alloys. J Mater Inf 2021;1:3.

22. Liu P, Huang H, Antonov S, et al. Machine learning assisted design of γ’-strengthened Co-base superalloys with multi-performance optimization. npj Comput Mater 2020;6:334.

23. Hu X, Zhao J, Li J, Wang Z, Chen Y, Wang J. Global-oriented strategy for searching ultrastrength martensitic stainless steels. Adv Theor Simul 2022;5:2100411.

24. Hu X, Chen Y, Lu J, et al. Three-step learning strategy for designing 15Cr ferritic steels with enhanced strength and plasticity at elevated temperature. J Mater Sci Technol 2023;164:79-94.

25. Xue D, Balachandran PV, Hogden J, Theiler J, Xue D, Lookman T. Accelerated search for materials with targeted properties by adaptive design. Nat Commun 2016;7:11241.

26. Yuan R, Liu Z, Balachandran PV, et al. Accelerated discovery of large electrostrains in BaTiO3-based piezoelectrics using active learning. Adv Mater 2018;30:1702884.

27. Hu X, Zhao J, Chen Y, Li J, Wang Z, Wang J. Continually reactivating iterative-projection method for instantiating microstructure from two-point statistics. Acta Mater 2022;238:118230.

28. Jiang M, Hu X, Li J, Wang Z, Wang J. An interface-oriented data-driven scheme applying into eutectic patterns evolution. Mater Design 2022;223:111222.

29. Hu X, Zhao J, Chen Y, et al. Structure-property modeling scheme based on optimized microstructural information by two-point statistics and principal component analysis. J Mater Inf 2022;2:5.

30. Xue D, Xue D, Yuan R, et al. An informatics approach to transformation temperatures of NiTi-based shape memory alloys. Acta Mater 2017;125:532-41.

31. Zhang H, Fu H, Zhu S, Yong W, Xie J. Machine learning assisted composition effective design for precipitation strengthened copper alloys. Acta Mater 2021;215:117118.

32. Goswami S, Anitescu C, Chakraborty S, Rabczuk T. Transfer learning enhanced physics informed neural network for phase-field modeling of fracture. Theor Appl Fract Mec 2020;106:102447.

33. Samaniego E, Anitescu C, Goswami S, et al. An energy approach to the solution of partial differential equations in computational mechanics via machine learning: concepts, implementation and applications. Comput Method Appl M 2020;362:112790.

34. Wei S, Kim J, Tasan CC. Boundary micro-cracking in metastable Fe45Mn35Co10Cr10 high-entropy alloys. Acta Mater 2019;168:76-86.

35. Ebrahimi M, Liu G, Li C, et al. Characteristic investigation of trilayered Cu/Al8011/Al1060 composite: interface morphology, microstructure, and in-situ tensile deformation. Prog Nat Sci Mater Int 2021;31:679-87.

36. Yan D, Tasan CC, Raabe D. High resolution in situ mapping of microstrain and microstructure evolution reveals damage resistance criteria in dual phase steels. Acta Mater 2015;96:399-409.

37. Lebensohn RA, Kanjarla AK, Eisenlohr P. An elasto-viscoplastic formulation based on fast Fourier transforms for the prediction of micromechanical fields in polycrystalline materials. Int J Plasticity 2012;32-3:59-69.

38. Lebensohn RA, Tomé CN, Castañeda PP. Self-consistent modelling of the mechanical behaviour of viscoplastic polycrystals incorporating intragranular field fluctuations. Philos Mag 2007;87:4287-322.

39. Lebensohn RA, Brenner R, Castelnau O, Rollett AD. Orientation image-based micromechanical modelling of subgrain texture evolution in polycrystalline copper. Acta Mater 2008;56:3914-26.

40. Eisenlohr P, Diehl M, Lebensohn R, Roters F. A spectral method solution to crystal elasto-viscoplasticity at finite strains. Int J Plasticity 2013;46:37-53.

41. Michel JC, Moulinec H, Suquet P. A computational method based on augmented lagrangians and fast Fourier transforms for composites with high contrast. Comput Model Eng Sci 2000;2:79-88.

42. James G, Witten D, Hastie T, Tibshirani R, Taylor J. Linear regression. An introduction to statistical learning. Cham: Springer International Publishing; 2023. pp. 69-134.

43. Cule E, De Iorio M. Ridge regression in prediction problems: automatic choice of the ridge parameter. Genet Epidemiol 2013;37:704-14.

44. McDonald GC. Ridge regression. WIREs Comp Stats 2009;1:93-100.

45. Moguerza JM, Muñoz A. Support vector machines with applications. Statist Sci 2006;21:322-36.

46. Pisner DA, Schnyer DM. Chapter 6 - Support vector machine. In: Machine learning. Elsevier; 2020. pp. 101-21.

47. Auret L, Aldrich C. Interpretation of nonlinear relationships between process variables by use of random forests. Miner Eng 2012;35:27-42.

48. Rigatti SJ. Random forest. J Insur Med 2017;47:31-9.

49. Deringer VL, Bartók AP, Bernstein N, Wilkins DM, Ceriotti M, Csányi G. Gaussian process regression for materials and molecules. Chem Rev 2021;121:10073-141.

50. Schulz E, Speekenbrink M, Krause A. A tutorial on Gaussian process regression: modelling, exploring, and exploiting functions. J Math Psychol 2018;85:1-16.

51. Wang J. An intuitive tutorial to Gaussian process regression. Comput Sci Eng 2023;25:4-11.

52. Jiang M, Chen Y, Yang Z, et al. Crystal plasticity modeling of deformation behavior of Al–Al2Cu eutectics based on high-fidelity representative microstructures. J Mater Res Technol 2024;29:5259-70.

53. Reed R. Aluminium 2. A review of deformation properties of high purity aluminium and dilute aluminium alloys. Cryogenics 1972;12:259-91.

54. Pham HH, Williams ME, Mahaffey P, Radovic M, Arroyave R, Cagin T. Finite-temperature elasticity of fcc Al: atomistic simulations and ultrasonic measurements. Phys Rev B 2011;84:064101.

55. Eshelman FR, Smith JF. Single-crystal elastic constants of Al2Cu. J Appl Phys 1978;49:3284-8.

56. Roters F, Eisenlohr P, Hantcherli L, Tjahjanto D, Bieler T, Raabe D. Overview of constitutive laws, kinematics, homogenization and multiscale methods in crystal plasticity finite-element modeling: theory, experiments, applications. Acta Mater 2010;58:1152-211.

57. Segurado J, Lebensohn RA, Llorca J, Tomé CN. Multiscale modeling of plasticity based on embedding the viscoplastic self-consistent formulation in implicit finite elements. Int J Plasticity 2012;28:124-40.

58. Balachandran PV, Xue D, Theiler J, Hogden J, Lookman T. Adaptive strategies for materials design using uncertainties. Sci Rep 2016;6:19660.

59. Simmons G. Single crystal elastic constants and calculated aggregate properties. Southern methodist university press; 1965.

60. Mahata A, Asle Zaeem M. Effects of solidification defects on nanoscale mechanical properties of rapid directionally solidified Al-Cu Alloy: a large scale molecular dynamics study. J Cryst Growth 2019;527:125255.

61. Davidson CJ, Smith IO, Chadwick GA. Effect of heat treatment and interlamellar spacing on the tensile deformation of the aligned Al CuAl2 eutectic. Acta Metall 1980;28:61-73.

62. Dey BN, Tyson WR. Plastic deformation of CuAl2. Phys Stat Sol 1972;9:215-21.

63. Cheng J, Yun Y, Rui J. Enhanced tensile plasticity in ultrafine lamellar eutectic Al-CuBased composites with α-Al dendrites prepared by progressive solidification. Appl Sci 2019;9:3922.

64. Pattnaik A, Lawley A. Deformation and fracture in AI-CuAl2 eutectic composites. Metall Trans 1971;2:1529-36.

65. Zhang X, Li S, Guo X, Wang H, Yu Q, Wu P. Effects of texture and twinning on the torsional behavior of magnesium alloy solid rod: a crystal plasticity approach in comparison with uniaxial tension/compression. Int J Mech Sci 2021;191:106062.

66. Wang H, Zhang X, Wu W, et al. On the torsional and coupled torsion-tension/compression behavior of magnesium alloy solid rod: a crystal plasticity evaluation. Int J Plasticity 2022;151:103213.

67. Xu Y. A non-local methodology for geometrically necessary dislocations and application to crack tips. Int J Plasticity 2021;140:102970.

68. Yu X, Xu Y, Morales-espejel G, Dunne F, Dini D. On the importance of crystal plasticity finite element discretisation for the identification of crack initiation in RCF using energy-based criteria. Comput Mater Sci 2024;232:112651.

69. Garatti S, Bittanti S. A new paradigm for parameter estimation in system modeling. Adapt Control Signal 2013;27:667-87.

70. Sin G, De Pauw DJ, Weijers S, Vanrolleghem PA. An efficient approach to automate the manual trial and error calibration of activated sludge models. Biotechnol Bioeng 2008;100:516-28.

71. Shahmardani M, Vajragupta N, Hartmaier A. Robust optimization scheme for inverse method for crystal plasticity model parametrization. Materials 2020;13:735.

72. Depriester D, Goulmy J, Barrallier L. Crystal plasticity simulations of in situ tensile tests: a two-step inverse method for identification of CP parameters, and assessment of CPFEM capabilities. Int J Plasticity 2023;168:103695.

73. Tarantola A. Inverse problem theory and methods for model parameter estimation. 2005.

Cite This Article

Research Article
Open Access
Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning
Meirong Jiang, ... Jincheng Wang

How to Cite

Jiang, M.; Hu, X.; Xing, C.; Yang, Z.; Chen, Y.; Li, J.; Wang, Z.; Wang, J. Fast parameter optimization for high-fidelity crystal plasticity simulation using active learning. J. Mater. Inf. 2024, 4, 22. http://dx.doi.org/10.20517/jmi.2024.31

Download Citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click on download.

Export Citation File:

Type of Import

Tips on Downloading Citation

This feature enables you to download the bibliographic information (also called citation data, header data, or metadata) for the articles on our site.

Citation Manager File Format

Use the radio buttons to choose how to format the bibliographic data you're harvesting. Several citation manager formats are available, including EndNote and BibTex.

Type of Import

If you have citation management software installed on your computer your Web browser should be able to import metadata directly into your reference database.

Direct Import: When the Direct Import option is selected (the default state), a dialogue box will give you the option to Save or Open the downloaded citation data. Choosing Open will either launch your citation manager or give you a choice of applications with which to use the metadata. The Save option saves the file locally for later use.

Indirect Import: When the Indirect Import option is selected, the metadata is displayed and may be copied and pasted as needed.

About This Article

Special Issue

This article belongs to the Special Issue Application of Machine Learning to Crystal Structure Prediction
© The Author(s) 2024. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Data & Comments

Data

Views
244
Downloads
54
Citations
0
Comments
0
0

Comments

Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at support@oaepublish.com.

0
Download PDF
Share This Article
Scan the QR code for reading!
See Updates
Contents
Figures
Related
Journal of Materials Informatics
ISSN 2770-372X (Online)
Follow Us

Portico

All published articles are preserved here permanently:

https://www.portico.org/publishers/oae/

Portico

All published articles are preserved here permanently:

https://www.portico.org/publishers/oae/