Abstract
Thermal accumulation is a characteristic feature of additive manufacturing (AM), which can lead to layer-dependent thermal histories, microstructural heterogeneity, and property variations along the build direction. To improve the consistency of melt pool depth in wire arc AM (WAAM), an artificial intelligence (AI)-assisted heat input design was proposed and evaluated through numerical simulation and experimental feasibility tests. A moving heat source model was first established to calculate the melt pool depth, and the simulated substrate penetration showed a mean absolute percentage error of 5.65% compared with the experimental measurements. The layer thickness and layer number in the numerical simulations were used as the inputs of the neural network, while the corresponding current was selected as the output. The correlation coefficient between the predicted and target current values was 0.99944. Compared with conventional WAAM under fixed heat input, where the melt pool depth varied from 2.72 to 6.67 mm, the proposed AI-assisted heat input design reduced the melt pool depth range to 1.92–2.72 mm under simulated WAAM conditions. In addition, WAAM deposition experiments using the AI-predicted layer-wise current schedule confirmed the practical implementability and macro-forming feasibility of the proposed strategy.
Keywords
Introduction
Additive manufacturing (AM) is quickly applied to different industries due to its high ability to be combined with computers.1–6 AM can be classified into different categories according to the uses of powder/wire and the applied heat source like laser, electron beam, or wire arc. Generally, laser-direct energy deposition, powder bed fusion (laser or electron beam), and wire arc AM (WAAM; electric arc or laser and electric arc) are used to fabricate metal components.7–9 Among all the AM technologies, WAAM is a type of AM using wire as feeding material and electric arc as heat source for the fabrication layer by layer.10–12 In comparison with the traditional manufacturing technologies, WAAM can reduce production time by 40–60%, with a rework time reduction of 15–20%, depending on the size and the complexity of the components. 13
In AM, thermal accumulation generally increases with build height and leads to layer-dependent thermal histories.14–16 As a result, the microstructure and related properties may vary along the build direction.17–19 A commonly observed tendency is that the lower region of the build, which is strongly affected by heat extraction through the substrate, exhibits a relatively different solidification condition from the middle region, where repeated thermal cycling and thermal accumulation become more pronounced. In the upper region, the final deposited layers experience a distinct thermal history because no subsequent layers are deposited above them. These differences in thermal conditions can result in variations in grain morphology, remelting behavior, and microstructural heterogeneity along the build direction.20–24 Such phenomena can be further understood from the perspectives of undercooling, nucleation, and thermal cycling.25–29 For clarity, the typical relationship between build height, thermal condition, and microstructural distribution is summarized in Table 1. These build-direction-dependent thermal and microstructural variations indicate that heat input regulation is essential for improving melt pool consistency and reducing microstructural heterogeneity during WAAM.
Typical Relationship Between Build-Direction Position, Thermal Condition, and Microstructural Characteristics in Additive Manufacturing
The evolution of the melting pool depth can be well visualized by numerical methods with different kinds of heat source models.30–34 The traditional heat source models like the Gaussian heat source model and double ellipsoid model used in welding simulations35–37 can also be directly applied to AM for the temperature visualizations.38–42 All temperature features, including maximum temperature rises, cooling rates, and the temperature gradients, can be well revealed by the numerical models.43–45 In particular, particle features46–48 and melt pool motions49–51 can be further incorporated into the heat source model to highlight key features of AM. The details of temperature variations and thermal accumulations are well shown by numerical models. However, further work is needed to determine how to control thermal accumulation and ensure melt pool consistency along the build height.
Recent developments in artificial intelligence (AI) provide possible solutions to the challenges in AM. Xie et al. 52 used a data-driven method to predict the melt pool geometries, in which simulation cases are used to form the database and Gaussian regression model is used to predict the keyhole sizes. Wang et al. 53 converted data like point clouds into control-field representations, and the map rules between control and property are set. Both the initial data and the mapping rules form the multi-information model to generate property data for industrial manufacturing. Mozaffar et al. 54 used neural networks for graph-based representation to obtain the thermal responses in AM, and the tests on directed energy deposition reveal the accuracy of the developed deep learning architecture on the prediction of thermal histories. Popova et al. 55 developed four-step data workflows to predict the process–structure relationships, in which data preprocessing, microstructure, dimensionality reduction, and extraction/validation are incorporated. The microstructures are predicted by the Potts Monte Carlo method, and the process parameters can be correlated with the final microstructures using the data-driven method. Herriott et al. 56 applied machine-learning and deep-learning methods to predict the microstructure-dependent mechanical properties of additively manufactured metals. Their study showed that AI-based models can effectively capture the relationship between microstructural features and mechanical properties, thereby demonstrating the potential of data-driven methods for process–microstructure–property analysis in AM. Xie et al. 57 used convolutional neural network integrating wavelet transforms to predict location-dependent mechanical properties based on the input temperature histories. The dominant mechanistic features like temperature ranges and thermal frequencies can be revealed by the developed data-driven computational framework. Nguyen et al. 58 presented a novel end-to-end AI-enhanced metal additive manufacturing lifecycle framework spanning the stages of design, build, post-processing, and end-of-life, which enables a more comprehensive evaluation of AI applications across the broader AM field. Representative applications of AI in AM are summarized in Table 2. Caggiano et al. 59 presented a process qualification-oriented, data-driven framework for WAAM, integrating qualification data, process monitoring, and feedback control. Xiong et al. 60 presented an adaptive control system for maintaining a constant distance between the nozzle and the top surface during GMAW-based layer AM. Love et al. 61 developed a vision-based monitoring system integrated with a reinforcement learning framework to enable intelligent in situ monitoring and control of WAAM. Xia et al. 62 proposed a deep-learning-based vision monitoring method for diagnosing different anomalies during the WAAM process. AI can be used to establish data-driven relationships between process parameters and process responses in AM.63–67 In this work, AI is employed to predict the current required under different layer-number and layer-thickness conditions for improving melt pool depth consistency in WAAM. Accordingly, the present study focuses on an AI-assisted offline heat input design for improving melt pool depth consistency in WAAM.
Representative Applications of Artificial Intelligence in Additive Manufacturing
AI, artificial intelligence; CNN, convolutional neural network; MAM, metal additive manufacturing.
Recent studies have demonstrated the strong potential of AI in AM for applications such as prediction, monitoring, and process–structure/property mapping. However, comparatively fewer studies have focused on the layer-by-layer design of heat input for mitigating thermal accumulation and improving melt pool depth consistency during WAAM deposition. This gap is important because thermal accumulation along the build direction leads to layer-dependent thermal histories, melt pool variations, and microstructural heterogeneity, which in turn affect the quality and stability of deposited components. To address this issue, the present study develops an AI-assisted heat input design using a backpropagation neural network for the offline layer-by-layer design of heat input in WAAM. The specific objective is to establish the relationship among layer number, layer thickness, and the required current to maintain a more consistent melt pool depth under thermal accumulation conditions. In this way, the present work provides a focused step toward AI-assisted heat input regulation in WAAM while also offering a basis for future experimentally validated control strategies.
Models and Methods
Model description
The finite element model is shown in Figure 1a. The substrate dimensions are 100 × 50 × 6 mm. The dimensions of the 20 deposited layers are 90 × 8 × 0.6 mm. The mesh is generated with linear hexahedral element DC3D8. The melt pool is shown in Figure 1b.

Numerical simulation:
Both the substrate and the deposited materials are Ti-6Al-4V. The temperature-dependent material properties of Ti-6Al-4V are listed in Table 3. 68 The mechanical properties of the deposited material are assumed to be identical to the substrate. The current is 180 A, the voltage is 15 V, and the scanning speed is 300 mm/min.
Temperature-Dependent Material Properties of Ti-6Al-4V
Heat transfer
The transient temperature field T(x, y, z, t) can be determined by the solution to the 3D heat conduction equation along with the relevant boundary conditions. The energy balance equation is as follows:
A double-ellipsoid heat source model is selected to simulate the heat input in this WAAM model. The power density in the front ellipsoid is given by the following equation:
The power density in the rear ellipsoid is given by the following equation:

Double-ellipsoid heat source shape parameters.
Usually, the thermal input q0 is constant during the AM process. But the temperature can be increased with the increase in building height due to thermal accumulations. To solve this problem, a BP neural network method is used to determine the layer-wise current required for heat input adjustment. The heat input is adjusted by changing the current for different deposited layers. In this way, the heat input can be varied by adjusting the layer-wise current. The heat input is determined by voltage, current, and thermal input efficiency.
The temperature field was not directly measured in the experiment. Instead, it was modeled using a 3D transient heat conduction model with a Goldak double-ellipsoidal heat source and temperature-dependent material properties. Accordingly, the present study uses indirect thermal validation through comparison of fusion-zone geometry and substrate penetration, rather than direct in situ temperature-field measurement.
For the deposition process, the elements assigned to the deposited layers were initially inactive and were activated sequentially along the prescribed scanning path as the moving heat source advanced. At each time increment, the volumetric heat input was applied only to the currently active elements, thereby ensuring a consistent coupling between material addition, heat input, and the layer-by-layer deposition sequence.71,72
BP neural network
The BP network is a type of multilayer feedforward network with unidirectional propagation, with the backpropagation algorithm. The flowchart of the BP algorithm is shown in Figure 3a. It is composed of an input layer, a hidden layer, and an output layer, as illustrated in Figure 3b. The layer thickness (h) and the number of layers (n) are given as input parameters. The corresponding current (I) through each layer is selected to be the output parameter. The BP neural network was selected in this study as a practical nonlinear regression model. The mapping from layer number and layer thickness to the required current is highly nonlinear, particularly because of the sharp variation in the early layers and the gradual stabilization in the later layers. In such a case, polynomial fitting would require a predefined functional form, whereas the BP neural network can learn the nonlinear relationship directly from the dataset. Moreover, previous studies on AM have reported that neural-network-based models are effective for bead-geometry-related prediction and may outperform conventional regression-based models when the underlying process relationship is strongly nonlinear.73–76 Therefore, the BP neural network was used here to construct the layer-by-layer current prediction model.

BP neural network:
The core training process of the BP neural network is divided into two stages: forward propagation and backward propagation. In forward propagation, the input data are processed layer by layer through the network’s weights, the bias calculations, and the activation functions to generate the necessary output data. In backward propagation, the error between the predicted value and target values is propagated backward from the output layer to the input layer. The weights and the bias are then updated using the gradient descent method to optimize the network parameters, reducing the error.
The neural network structure consists of an input layer, a hidden layer, and an output layer. In this model, h = (h1, h2, …, hn) and n = (n1, n2, …, nn) are considered to be the input data. I = (I1, I2, …, In) represents the data from the output layer. The parameter θ represents the weight vector connecting the input layer to the hidden layer. θij represents the weight connecting the i-th input node to the j-th hidden layer node. λ represents the weight vector connecting the hidden layer to the output layer, specifically corresponding to the weight coefficients from the hidden layer nodes to the output layer nodes. The number of neurons in the hidden layer is 3. For instance, λjk represents the weight connecting the j-th hidden layer node to the k-th output layer node.
The computed output of the output nodes is given by,
The function f in the above equation is the activation function being selected to be the form below:
The error E is generated and needs to be updated in the backward propagation:
Deposition process
Ti-6Al-4V is used as both the deposited material and the substrate material. The experiments were carried out using a common MIG welder. The wire diameter is 1.2 mm, with a wire feeding rate of 1 m/min. The scanning speed is set to 300 mm/min, and the welding current is maintained at 180 A. The welding current was directly prescribed in the WAAM system during deposition. In contrast, the arc voltage was not directly recorded in the experiment. Therefore, a constant nominal voltage was adopted as the model input for heat input calculation. The adopted value was selected to be consistent with comparable Ti-6Al-4V WAAM conditions reported in the literature, where average arc voltages of ∼15 V. 77 The efficiency η is assumed to be 0.6. The experiments were performed under continuous deposition conditions. The interval between adjacent layers was short, and no intentional interlayer cooling was applied during the build. The ambient temperature was set to 25°C, and the convection coefficient was taken as 20 W/(m2·K). The substrate dimensions are 200 × 100 × 6 mm. 1-layer, 3-layer, and 5-layer specimens are fabricated on the substrate, as illustrated in Figure 4a and b. The 1-layer specimen is measured approximately to be 80 mm in length, 8 mm in width, and 0.6 mm in height. The 3-layer specimen is 80 × 8 × 1.8 mm, while the 5-layer specimen is 90 × 8 × 3 mm. The maximum width and the total deposited height of the additively manufactured wall were measured from the deposited geometry. In the numerical model, the deposition domain was constructed according to the measured maximum width and total height. For simplification, a uniform layer width was assumed for all deposited layers, and the layer height for each activated layer was taken as the total deposited height divided by the number of deposited layers. The corresponding deposition elements were initially defined as inactive and were activated sequentially along the predefined scanning path during the transient analysis. The deposited layers are cut on the cross-section at a distance of 3 cm along the scanning direction to obtain the specimens, as shown in Figure 4c–e. The cross-sections of the specimens are subsequently grounded, polished, and etched.

Wire arc additive manufacturing:
Results and Discussions
Experimental verification
The numerical and experimental results were correlated through the comparison of cross-sectional fusion-zone geometry and substrate penetration. Specifically, the experimentally measured and simulated substrate penetration values were compared for the 1-layer, 3-layer, and 5-layer specimens at the selected cross-section. The resulting mean absolute percentage error (MAPE) was used to quantify the agreement between the numerical model and the experimental observations.
The cross-sections of the deposited specimens are shown in Figure 5. The overall fusion depths for the 1-layer, 3-layer, and 5-layer specimens are 2.26, 3.66, and 4.94 mm, respectively. The overall fusion depth increases with the number of deposited layers. The substrate penetrations are 1.71, 2.11, and 2.12 mm, respectively. The substrate penetration remains essentially unchanged between the 3-layer and 5-layer specimens.

Cross-sections of the experimental specimens:
The numerical results on the cross-sections at different scan lengths (Ls) are shown in Figure 6. The melt pool in the simulation is determined according to the melting point of the deposited material, that is, 1660°C. At a scan length of 1 cm, the melt pool depths for the five layers are 2.34, 2.83, 3.17, 3.50, and 4.07 mm, respectively. When the scan length is increased to 5 cm, these values increase to 2.45, 3.25, 3.80, 4.17, and 4.47 mm. For substrate penetration, the values are 1.74, 1.63, 1.37, 1.10, and 1.07 mm at 1 cm, increasing to 1.85, 2.05, 2.00, 1.77, and 1.47 mm at 5 cm. Both the melt pool depth and the substrate penetration increase with increasing scan length.

Numerical simulation of cross-sections at different scan lengths:
Figure 7 shows a comparison of substrate penetration between experimental and simulated results at a scan length of 3 cm. For layers 1, 3, and 5, the experimental values are 1.71, 2.11, and 2.12 mm, respectively, while the simulated values are 1.83, 2.01, and 2.01 mm. The MAPE of 5.65% indicates good agreement between the simulation and experimental results.

Comparison of substrate penetration between experimental and simulated results:
To verify the credibility of the numerical model, experimental measurements were compared with simulation results, as shown in Figure 8. In the experiments, the widths of layers 1, 3, and 5 range from 8.21 to 8.52 mm, while the simulated width is consistently 8 mm. The measured heights of layers 1, 3, and 5 are 0.65, 1.98, and 3.18 mm, respectively, whereas the corresponding simulated values are 0.6, 1.8, and 3.0 mm. The comparison shows that the MAPE for width and height is 3.94% and 7.48%, respectively.

Comparison between experiment and simulation:
Thermal accumulation under fixed heat input
Building upon the previous experimental cases, the numerical analysis was extended to a 20-layer simulation case to investigate the cumulative effect of thermal accumulation at a larger build height. Similar approaches have been used in previous WAAM simulation studies, where models validated on simpler experimental cases were further applied to analyze multilayer builds beyond the experimentally tested layer numbers.78,79
Figure 9 shows the melt pool depth evolution along the scanning path within a layer. When the scan length increases from 1 to 2 cm, the melt pool depth of layer 2 increases from 2.83 to 3.11 mm, representing a 9.89% increase. The melt pool depth of layer 4 increases from 3.50 to 3.96 mm, with a 13.14% increase. Layer 20 exhibits an increase from 4.91 to 5.53 mm, corresponding to a 12.63% increase. When the scan length increases from 2 to 8 cm, the melt pool depth of layer 2 changes by 0.17 mm (5.48), that of layer 4 changes by 0.27 mm (6.90%), and that of layer 20 changes by 0.40 mm (7.42%). From 1 to 2 cm along the scanning path, the melt pool depth increases noticeably. This indicates that the local thermal state is still evolving in the early stage of deposition, because the material at these positions has experienced different elapsed deposition times and different degrees of prior thermal accumulation. In contrast, from 2 to 8 cm, the variation in melt pool depth becomes much smaller, suggesting that the local thermal state along the scanning path has approached a relatively stable condition. Therefore, the melt pool depth at the midpoint of the scan path (4.5 cm) is selected as the representative value for each layer.

Melt pool depth evolution along the scanning path within a layer.
Figure 10 shows the melt pool depth evolution with increasing layer number. When the layer number increases from 1 to 2, the melt pool depth at a scan length of 1 cm increases from 2.34 to 2.83 mm, representing a 20.94% increase. At a scan length of 2 cm, the melt pool depth increases from 2.42 to 3.11 mm, with a 28.51% increase. At a scan length of 8 cm, the melt pool depth increases from 2.43 to 3.10 mm, corresponding to a 27.57% increase. When the layer number increases from 19 to 20, the melt pool depth at a scan length of 1 cm increases from 4.88 to 4.91 mm, with an increase of 0.61%. At a scan length of 2 cm, the melt pool depth increases from 5.51 to 5.53 mm, representing a 0.36% increase. At a scan length of 8 cm, the melt pool depth increases from 5.38 to 5.39 mm, with an increase of 0.19%. Thermal accumulation is relatively weak at a scan length of 1 cm. With the increase in the layer number, the melt pool depth still increases, but the increase rate slows down, indicating progressively weaker thermal accumulation effects.

Melt pool depth evolution with increasing layer number.
Figure 11 shows the variation in the maximum temperatures of the substrate and the deposited layers with increasing layer number. For layer 1, the local peak temperature of the substrate-side elements adjacent to the fusion boundary reaches 2244°C. As the layer number increases, the local peak value decreases, reaching 1306°C at layer 20. It should be noted that these values represent model-predicted transient local peaks near the fusion region, rather than the average temperature of the substrate or the interpass temperature. The simulation was carried out under a continuous-deposition condition without active interlayer cooling, which promotes stronger heat accumulation and repeated reheating during the build.80–83 In contrast, for layer 1, the maximum temperature of the deposited layer is 2321°C. For layer 2, the maximum temperature of the deposited layer rises to 2488°C, representing a 7.20% increase. For layer 19, the maximum temperature of the deposited layer reaches 2820°C. For layer 20, it further rises to 2822°C, with a 0.07% increase. Heat accumulates in the deposited layer,84–86 causing the maximum temperature to increase with the layer number, but the rate of increase slows down.

Variation in the maximum temperatures of the substrate and deposited layers with increasing layer number.
The temperature cycles extracted at a fixed point located at the midpoint of the scanning path on the top surface of layers 1, 5, 10, and 15 are shown in Figure 12. For layer 1, the temperature reaches 2301°C during its deposition. Owing to subsequent thermal accumulation, the temperature further increases to 2313°C during the deposition of layer 2 and then decreases. Layer 1 undergoes eight remelting cycles. For layer 5, the temperature reaches the maximum of 2663°C during its deposition and subsequently decreases. Layer 5 undergoes nine remelting cycles. Layers 10 and 15 reach their peak temperatures during their respective depositions, at 2760°C and 2790°C. Layer 10 undergoes nine remelting cycles. Since a total of 20 layers are deposited, the number of remelting cycles decreases for higher layers. Layer 15 undergoes five remelting cycles. Variations in the thermal cycles across different layers result in different melt pool depths. These variations lead to differences in grain morphology from layer to layer, which ultimately influences the mechanical properties. Therefore, controlling the melt pool depth of each layer is essential.

Temperature cycles of different layers.
The melt pool depth for different deposited layer thicknesses is shown in Figure 13. For layer thicknesses of 0.5, 0.6, 0.7, 0.8, and 0.9 mm, the melt pool depth increases from layer 1 to layer 2 by 0.82, 0.83, 0.86, 0.89, and 0.93 mm, respectively. In contrast, from layer 19 to layer 20, the corresponding increases are 0.03, 0.03, 0.03, 0.01, and 0.01 mm, respectively. The melt pool depth increases, but the rate of increase slows down. From layer 1 to layer 20, the melt pool depth increases by 3.11, 3.33, 3.53, 3.70, and 3.89 mm, respectively. With increasing layer thickness, the melt pool depth increases more significantly, suggesting enhanced thermal accumulation.

Melt pool depth for different layer thicknesses.
AI-assisted heat input design
To maintain the melt pool depth across different layers within a relatively stable range, an AI-assisted heat input design is proposed. The AI-assisted heat input design is implemented through a layer-by-layer decreasing current strategy, in which the current is assigned individually to each layer to keep the melt pool depth at approximately twice the deposited layer thickness. For layer thicknesses ranging from 0.5 to 0.9 mm, the designed current values are summarized in Table 4 and used as the training set for the neural network. The neural network is used to establish the relationship among layer number, layer thickness, and the required current. The proposed method is an offline layer-by-layer heat input design based on the current adjustment. Given the layer number and layer thickness, the AI-assisted heat input design predicts the corresponding current values for the deposited layers.
Current (A) per Layer of the Artificial Intelligence-Assisted Heat Input at Different Layer Thicknesses
Each combination of layer thickness and layer number corresponds to one current value, yielding a total of 100 samples. Among them, 70% is used as the training set, 15% as the validation set, and 15% as the test set. The layer thickness and layer number are input variables, while the current value is taken as the output. To improve the training performance, both input and output vectors are normalized prior to training. After normalization, the original data are transformed into dimensionless values ranging from −1 to 1, which can be expressed as:
To investigate whether the training set size is sufficient, the relative error (e) and correlation coefficient (ρ) between the true and identified current values for 10 different sizes (10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, and 100% of the total size) are analyzed. Training performance improves as ρ approaches 1 and e approaches 0.
Let It and Ii represent the true and identified current values in the test set, respectively; the correlation coefficient ρ(It, Ii) and relative error e(It, Ii) can be described as:
The correlation coefficients and relative errors for different training set sizes are shown in Figure 14. As the size of the training set increases, ρ gradually increases, while e correspondingly decreases. When the sample size exceeds 60, ρ stabilizes at approximately 0.997, approaching 1, and e stabilizes at about 0.026, approaching 0. Therefore, a dataset size of 100 is sufficient.

Correlation coefficient and relative error under different training set sizes.
The error performance curve during training is shown in Figure 15a. The training process runs for a total of 162 epochs. As the number of epochs increases, the training, validation, and test errors all decrease. The best validation performance is achieved at the 162nd epoch, with a mean squared error of 4.06 × 10−4.

BP neural network:
The training regression results are shown in Figure 15b. The apparent two-region distribution in Figure 15 originates from the layer-dependent current values in the dataset. Specifically, under all five thickness conditions, the first layer requires a substantially higher current than the subsequent layers. Consequently, the first-layer samples are grouped in the high-current region, while the majority of the later-layer samples fall into the low-current region, which makes the intermediate region less populated. Most data points are close to the target line. The correlation coefficient between the predicted and actual values of the dataset is 0.99944, which is close to 1 and indicates a strong linear relationship. This demonstrates that the prediction model has good data fitting capability for the training samples and achieves satisfactory training performance.
Results under AI-assisted heat input
The temperature cycles of different layers under the AI-assisted heat input are shown in Figure 16. Compared with fixed heat input, the AI-assisted heat input reduces the peak temperatures for all layers. For layer 5, the maximum temperature is 1772°C, which is 891°C lower than that under fixed heat input, and the number of remelting cycles decreases from nine to one. For layer 10, the maximum temperature is 1751°C, representing a reduction of 1009°C, and the remelting cycles decrease from nine to one. For layer 15, the maximum temperature is 1748°C, which is 995°C lower, and the remelting cycles decrease from five to one. The use of the AI-assisted heat input effectively controls heat input and maintains consistency of melt pool depth across layers. The remelting behavior across layers becomes more uniform, reducing differences in thermal history. This homogenized thermal cycling strongly influences grain growth behavior and microstructural evolution, thereby improving microstructural uniformity and the stability of mechanical properties.

Temperature cycles of different layers using the artificial intelligence-assisted heat input at a layer thickness of 0.6 mm.
To further verify the feasibility of the prediction model, the AI-assisted heat input design was used to calculate the current values for a layer thickness of 1.0 mm. Layer-wise current values are listed in Table 5. The melt pool morphologies obtained using fixed heat input and AI-assisted heat input at a layer thickness of 1.0 mm are shown in Figure 17. The melt pool depth varies from 2.72 to 6.67 mm in traditional WAAM and can be adjusted to 1.92 to 2.72 mm using the proposed AI-assisted heat input design. Melt pool depth consistency is improved using the AI-assisted heat input design. The results show that the MAPE, mean absolute error, root mean square error, and coefficient of variation of the root mean square error are 2.8158%, 0.0563%, 0.0600%, and 3.0022%, respectively. These metrics describe the prediction error from the perspectives of relative deviation, average absolute deviation, overall error level, and normalized error dispersion, respectively, indicating the predictive accuracy of the proposed method.
Layer-Wise Current Values Predicted by the Artificial Intelligence-Assisted Heat Input Design for the 1.0 mm Layer-Thickness Case

Melt pool morphologies using fixed heat input and artificial intelligence-assisted heat input.
Additional experimental deposition was carried out to further examine the practical implementability of the layer-wise current schedule listed in Table 5, as shown in Figure 18. Figure 18a shows the 20-layer WAAM wall deposited using the AI-predicted layer-wise current schedule. The wall could be continuously formed without severe macroscopic collapse or discontinuity, indicating that the predicted current range was practically implementable under the present deposition conditions. Figure 18b shows the cross-section of the deposited wall, and Figure 18c presents the stitched metallographic cross-section. The cross-section shows a continuous deposited structure and metallurgical bonding between the deposited wall and the substrate. The enlarged regions in Figures 18d–f further show the local microstructural features at different build-height positions. It should be noted that the individual melt pool boundaries of each deposited layer cannot be clearly distinguished in the final cross-section, mainly because repeated remelting and thermal cycling during multilayer WAAM obscure the original layer-wise fusion boundaries. Therefore, this experiment is not used as a direct quantitative validation of transient layer-wise melt pool depth. Instead, it provides macro-forming and process-feasibility evidence for the AI-assisted heat input design.

Additional experimental deposition using the layer-wise current schedule listed in Table 5:
Limitations and future research
Although the proposed AI-assisted heat input design improves the simulated consistency of melt pool depth, some limitations should be noted. First, the present study still mainly relies on numerical simulation for evaluating the transient, layer-wise melt pool depth. The experimental results validate the thermal model indirectly through substrate penetration and deposited geometry, and the additional 20-layer deposition experiment in Figure 18 demonstrates the macro-forming feasibility of the AI-predicted layer-wise current schedule. However, the transient, layer-wise melt pool depth was not directly measured. Second, the thermal model includes several simplifications, such as constant nominal arc voltage, assumed heat input efficiency, simplified boundary conditions, and nominal layer geometry. These assumptions may affect the quantitative accuracy of the simulation under real WAAM conditions. Third, the proposed method is an offline heat input design method based on layer-wise current adjustment, not a real-time closed-loop control strategy. The predicted current schedule does not yet include feedback from melt pool geometry, arc signals, bead height, or interlayer temperature.
Future work will focus on direct experimental validation using in situ melt pool monitoring, thermal imaging, high-speed vision, or pyrometer-based temperature measurement. In addition, current adjustment should be combined with wire feed rate, travel speed, interlayer cooling, and bead-height monitoring to move toward closed-loop WAAM process control.
Conclusions
This study developed an AI-assisted heat input design using a BP neural network to improve the layer-to-layer consistency of melt pool depth in WAAM. The main findings are summarized as follows:
Thermal accumulation becomes more pronounced with increasing build height in WAAM, resulting in layer-dependent thermal histories and an increase in melt pool depth under fixed heat input conditions. A relationship among layer number, layer thickness, and the required current was established using a BP neural network, thereby enabling offline layer-by-layer heat input design for improving melt pool depth consistency. For the additional validation case, the predicted melt pool depth showed an error of 2.82% relative to the target value, supporting the effectiveness of the proposed AI-assisted heat input design under simulated WAAM conditions. The proposed AI-assisted heat input design improves melt pool depth consistency, narrowing the depth variation from 2.72–6.67 mm in traditional WAAM to 1.92–2.72 mm.
Authors’ Contributions
H.Z.: Writing—original draft, data curation, and formal analysis. Z.Z.: Writing—reviewing and editing, supervision, and project administration.
Data Availability Statement
The raw/processed data required to reproduce these findings cannot be shared at this time because of technical or time limitations but are available from the authors upon reasonable request.
Footnotes
Acknowledgments
The authors would like to acknowledge the help from Prof. Hongyang Wang at Dalian University of Technology with the experimental validation.
Author Disclosure Statement
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this article.
Funding Information
This work was supported by the
