Next Article in Journal
Removal of Methylene Blue from Aqueous Solution by Mixture of Reused Silica Gel Desiccant and Natural Sand or Eggshell Waste
Previous Article in Journal
Preparation of Anorthite/Mullite In Situ and Phase Transformation in Porcelain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Laser Welding Penetration Monitoring Based on Time-Frequency Characterization of Acoustic Emission and CNN-LSTM Hybrid Network

1
School of Materials Science and Engineering, Shanghai University of Engineering Science, Shanghai 201620, China
2
Shanghai Collaborative Innovation Center of Laser of Manufacturing Technology, Shanghai 201620, China
3
School of Materials Science and Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
4
State Key Laboratory of Advanced Welding and Joining, Harbin Institute of Technology, Harbin 150001, China
5
Department of Engineering, Lancaster University, Lancaster LA1 4YW, UK
*
Author to whom correspondence should be addressed.
Materials 2023, 16(4), 1614; https://doi.org/10.3390/ma16041614
Submission received: 7 January 2023 / Revised: 8 February 2023 / Accepted: 12 February 2023 / Published: 15 February 2023
(This article belongs to the Section Manufacturing Processes and Systems)

Abstract

:
In-process penetration monitoring of the pulsed laser welding process remains a great challenge for achieving uniform and reproducible products due to the highly complex nature of the keyhole dynamics within the intense laser-metal interactions. The main purpose of this study is to investigate the feasibility of acoustic emission (AE) measurement for penetration monitoring based on acoustic wave characteristics and deep learning. Firstly, a series of laser welding experiments on aluminum alloys were conducted using high-speed photography and AE techniques. This allowed us to in-situ visualize the complete keyhole dynamics and elucidate the generation mechanism of acoustic waves originating from pressure fluctuations at the keyhole wall. Then, an adaptive time-frequency technique namely VMD (Variational Mode Decomposition) was proposed to characterize the acoustic energy distribution among the nine subsignals with low-frequency and high-frequency components under different welding penetrations. Lastly, a novel hybrid model combing CNN (Convolutional Neural Network) and LSTM (Long Short Term Memory) was designed to deeply mine the spatial and temporal acoustic features from the extracted frequency components. Extensive experiments demonstrate that our proposed approach yields a remarkable classification performance with a test accuracy of 99.8% and a standard deviation of 0.21, which obtains a high recognition rate. This work is a new paradigm in the digitization and intelligence of the laser welding process and contributes to an alternative way of developing an efficient end-to-end penetration monitoring system.

1. Introduction

Aluminium alloys have the advantages of high specific strength, good corrosion resistance, and high thermal conductivity and are widely used in aerospace and automotive manufacturing [1]. 6XXX series aluminium alloys are widely used in the automotive manufacturing industry, usually in the fields of battery manufacturing and automotive chassis manufacturing by joining in a welding way.
However, due to the poor heat resistance and strong thermal conductivity of aluminium alloys, defects are very easy to form in welding [2]. Traditional arc welding for welding aluminium alloys has large welding deformation and low productivity, while electron beam welding of aluminium alloys requires a guaranteed vacuum environment, resulting in high costs and low efficiency. However, laser welding has a high energy density and high welding speed, which allow more precise control of the heat input. Therefore, laser welding in aluminium alloys has a high energy density, low heat input, low welding distortion, and a large depth to width ratio.
Therefore, laser welding processes with such advantages as high energy, high precision, and high efficiency have been widely applied in high-end industry fields such as aerospace, automotive, shipbuilding, etc. [1]. With its precise regulation of laser energy, the pulsed mode laser welding technique is particularly suitable for joining thin aluminium alloys with high reflectivity and high thermal conductivity [2]. In general, the combination of pulse frequency, pulse width, peak power, and other parameters can directly affect the welding process stability. Particularly, the intense interaction between the high-energy laser beam and aluminium alloy easily drills an extremely unstable vapor capillary (referred to as a “keyhole”) into the workpiece. As such, any tiny change of operation conditions can easily produce the instability phenomena of melt pool morphology and intense keyhole oscillation, which may easily cause the welding defects occurring in the welding process [3]. Therein, insufficient penetration weakens the mechanical properties of the welded parts, while excessive penetration can be characterized by sagging, holes, and bum-through.
To understand the keyhole dynamics and further improve the product quality, a series of in-situ monitoring approaches, including thermal [4,5], photoelectric [6], and optical sensors [7,8], have been proposed in industry and academia. The design idea of mainstream monitoring solutions lies in the utilization of the phenomenon of keyhole formation, which carries various types of valuable information about weld quality. It is worth noting that acoustic waves play a key role in driving pores near the keyhole tip far enough away from the large thermal gradient field around the keyhole during the laser melting process [9].
Therefore, a promising acoustic emission (AE) monitoring system has received much attention recently due to its cost-efficient, in-line, and non-invasive measurement [10]. It can capture the strong acoustic wave generated from the laser-material interactions, propagating inside the workpiece and atmosphere. Due to this nature, on-line AE measurement can detect the pressure fluctuation generated by the molten pool pulsation, plasma/metal vapor generation, thermal stress, and keyhole oscillation, which are closely related to the depth of weld penetration. However, the key challenge in the development of a reliable and cost-effective acoustic monitoring system is the intricate dependency between the process parameters, the acoustic signatures, and the weld quality. A detailed investigation of the generation mechanism of acoustic waves and its relation to welding quality has not yet been assessed conclusively.
(1)
Acoustic emission monitoring
To overcome this issue, Gu et al. [11,12] applied a frequency analysis and a statistical approach to deal with the airborne sound signals acquired in the CO2 laser welding process. The results suggest that the band of frequency components is well associated with the keyhole oscillation, and these frequencies are distinctive when the welds are fully penetrated. Huang et al. [13,14] used a spectral subtraction method to reduce the background noise in the acquired acoustic signals and then applied a power spectrum method (Welch-Bartlett) to analyze the frequency characteristics of the acoustic signals. The relationship between acoustic features and various penetration depths can also be established based on time- and frequency-domain processing. Recently, Yusof et al. [15,16] extracted the time-frequency characteristics of the collected AE signals based on synchrosqueezed wavelet analysis during pulsed laser welding experiments and statistically correlated them with two process regimes, including half-penetration and full-penetration status.
To sum up, the observation of distinct AE signals could provide a base for developing a diagnostic/analytical technique for laser welding processes, and it also underpins the existing high potential of AE monitoring in laser processes. Although previous studies have shown that acoustic signals can effectively monitor the laser welding of mild steel, stainless steel, and high-strength steels, there is limited literature on the study of acoustic signals during the pulsed laser welding of aluminum alloys. Moreover, pulse mode laser welding could possibly emit non-stationary period acoustic signals when the pulsed laser beam intermittently interacts with the workpiece over very short time intervals. How to accurately extract the useful information for reflecting the keyhole oscillation and welding stability independently remains a challenge and warrants further research. Consequently, a good prior knowledge of the laser welding process is required to explain the acoustic generation mechanism and acoustic emission trends.
(2)
Machine learning modeling
Based on the afore-cited research, it could be summarized that much of the research up to now has attempted to reveal the significant correlation between the weld quality and the acoustic signatures. Basically, an accurate and quantitative classification of the weld penetration is also imperative as it could directly give feedback to the control system in the automated welding process [17,18]. Wu et al. [19] pointed out that the key term “monitoring” is applied to indicate in-situ information collection, feature extraction, and process modeling. Multiple approaches have been proposed, which are mostly based on analytic/numerical modeling via AE measurements and process data. Due to the high complexity and dynamic nature of laser-metal interactions, a large number of factors affect the final result of repetitive pulsed laser irradiation, and a wide range of relevant physical processes complicate a detailed analysis. Moreover, establishing a complex theoretical or numerical model is rather time-consuming and not suitable for real-time monitoring and quality control.
Therefore, proper acoustic in-process monitoring associated with data-driven machine learning (ML) and feature extraction is needed to increase the productivity and repeatability of laser welding. To our knowledge, a common ML-based algorithm, namely artificial neural networks (ANNs), is capable of dealing with complex databases, and has been widely used to perform tasks such as nonlinear approximation [20], prediction, and classification [21] in a dynamic and complicated industrial process, especially in the field of laser welding. For instance, Huang et al. [14] proposed a neural network and multiple regression method to characterize the depth of weld penetration based on the acquired signatures. The results show that the acoustic signatures can predict the penetration depth well under different laser welding parameters. Lee et al. [22] presented an in-process monitoring scheme for pulsed Nd:YAG laser spot welding based on the AE technique. The amplitudes of three different frequency ranges (100–200, 200–300, and 300–500 kHz) are used as input vectors, and three different process conditions (unsuccessful, successful, and over-welding) are clearly distinguished by utilizing the ANN-based classification model. However, the shallow ANNs usually require feature extraction of acoustic signals in time or frequency domain, which could not completely reflect the dynamic changes of the melt pool and the keyhole under non-stationary and complex conditions of the AE events.
To overcome the above difficulty, a recently developed convolutional neural network (CNN) is emerging as a promising deep learning (DL) framework since it can automatically learn and characterize the non-linear AE signatures through the multi-layer stacks. For instance, Shevchik et al. [23] studied the feasibility of using acoustic emission (AE) and Spectral Convolution Neural Network (SCNN) for in-situ and real-time quality monitoring by using the energy in narrow-frequency acoustic features. Pandiyan et al. [24] extracted the AE spectrograms using wavelet transforms and applied two CNN architectures, including VGG-16 and ResNets-18, to predict the build quality, including balling, LoF pores, conduction mode, and keyhole pores, during the laser additive manufacturing.
However, the current CNN methods are limited in exploring the spatial features of the two-dimensional AE signal spectrogram. There are few studies investigating the temporal characteristics of the one-dimensional (1-D) AE signal during the dynamic laser welding process. The newly developed DL framework, called long-short term memory (LSTM), has several advantages in dealing with time information in time-series AE signals. Since the CNN can extract much greater high-level spatial features, recent research attempted to combine the CNN and LSTM models to in-depth mine the spatio-temporal characteristics of signal data, in order to enhance the predicting accuracy and reliability.
For example, Zhang et al. [25] proposed a dynamic model that combines CNN and LSTM networks to account for the detected AE signals, which can accurately monitor the surface quality during laser shock peening. Shi et al. [26] developed a novel deep learning algorithm (BiConvLSTM) that combines CNN and LSTM networks for planetary gearbox fault diagnosis under different operating conditions. The BiConvLSTM can automatically and simultaneously extract spatio-temporal features from the condition monitoring data collected by three accelerometers and one tachometer sensor. Song et al. [27] developed a CNN-LSTM bearing fault diagnosis model using the VMD (Variational Mode Decomposition) technique to reconstruct the bearing signals, which was then fed into a CNN-LSTM network to extract spatio-temporal features simultaneously and obtain a good diagnostic result. However, there are no investigations in the state of the art yet that correlate the acoustic characteristics to keyhole dynamics and weld penetration depth based on CNN-LSTM. Driven by the intrinsic relationship between acoustic characteristics and welding quality, our study develops a hybrid CNN-LSTM model for penetration prediction to provide a new approach for in-situ monitoring of the laser welding process. The main contributions to this work are as follows:
  • Through investigating the correlations between the transient nature of the subsurface keyhole and the acoustic signals, the AE technique could better characterize the keyhole dynamics and weld penetration. It will also provide important guidance for understanding the complicated interaction mechanisms between the pulsed laser radiation and aluminum alloy and can be extended to other laser-based manufacturing scenarios.
  • In contrast to traditional time- or frequency-domain processing for non-stationary AE signals, an adaptive time-frequency technique called VMD was proposed, which can accurately distinguish between low-frequency and high-frequency components for better describing the weld penetration.
  • A novel CNN-LSTM hybrid model was proposed to deeply mine the spatial and temporal acoustic features from the extracted frequency components. It can improve the penetration predicting performance of AE sensing and provide a potential and reliable monitoring technology for the dynamic laser welding process.
The framework of this paper is organized as follows: Section 2 introduces the experimental setup and process analysis; Section 3 illustrates the AE frequency component extraction process using the VMD; Section 4 presents the spatio-temporal CNN-LSTM model for predicting the weld quality; Section 5 illustrates the performance results compared with other algorithms; the main conclusions and future study are stated in Section 6.

2. Methodical Approach

2.1. Experimental Setup and Data Acquisition

To reveal the relationship between AE signals and keyhole characteristics as well as the weld penetration during pulsed laser welding, a complete experimental platform was created, including: (i) a multi-information acquisition system, (ii) a pulsed mode laser welding system, and (iii) a traveling mechanism system, as illustrated in Figure 1. In our work, the laser welding system mainly contains a rectangular pulsed fiber laser device (TRUPULSE 556, λ = 1064 nm) with a 10 kW maximum output power. A laser head (Trumpf BEO D70, Ditzingen, Germany) with an inclination angle of 5° was applied to avoid the back reflection damaging the delivery optics, particularly for welding highly reflective aluminum alloys. Pure argon was used as the shielding gas, and a rectangular pulse type was selected with a beam diameter of 0.6 mm at the focal position. The traveling mechanism system was applied to drive the workpiece while the laser head remained fixed in the vertical position.
To gain a better understanding of the AE technique and to broaden its application for estimating different weld penetration, it is necessary to test its capability with a wide range of laser parameters. Thus, a series of different bead-on-plate experiments were conducted on aluminum alloy (6061) plates by modifying the laser power. Other parameters remained unchanged, as listed in Table 1. This procedure was repeated in order to create a data basis for the evaluation of the correlation between achieved penetration depth and AE with respect to the adjustable welding parameters. After laser welding, the weld’s quality was characterized postmortem by means of metallographic testing using an optical microscope setup.
The multi-information synchronous acquisition system consists of two subsystems: visual sensor and the ABAE sensor (BSWA TECH, Beijing, China, MA231 pPreamplifier). Therein, the visual sensing system contains a high-speed camera (PHANTOM VEO 710L) with a band-pass filter (810 nm) and an auxiliary illuminant (CAVILUX HF, 810 nm ± 10 nm), which could efficiently reduce the existing interferences of spatters and plasma plumes. The sampling frequency was set to 10,000 fps, and the exposure duration was 1 μs. The original lateral keyhole image was 512 × 384 pixels with a resolution of 0.022 mm/pixel. In addition, a butt-joint configuration consisting of an aluminum alloy (6061) with size of 100 mm (length) × 50 mm (width) × 2 mm (height), and transparent heat-resistant quartz glass with size of 100 mm (length) × 5 mm (width) × 2 mm (height), was constructed to directly observe the keyhole dynamics inside the melt pool through the high-speed imaging during the laser welding process. The ABAE signals (airborne sound pressure waves) were collected using a condenser microphone (BSWA TECH, MA231 Preamplifier) with a sampling frequency of 20 kHz. The microphone was oriented at an angle of 45° with the horizontal axis, and the distance between the laser-metal interaction and the microphone was 50 mm. The national standards for the chemical composition [28] and mechanical properties [29] of 6061 aluminium alloys are given in Table 2.

2.2. Detected AE Signals and Process Analysis

It is demonstrated that the weld penetration is of importance to weld formation/joint quality, and the backside weld width is usually proposed to evaluate the penetration status quantitatively. According to the actual weld widths, the degree of weld penetration can be categorized into three types, including: (i) none-penetration (NP), (ii) partial-penetration (PP), and (iii) full-penetration (FP), as depicted in Table 1. Figure 2 displays the transverse cross sections (perpendicular to the laser path) of the fusion zone and the corresponding AE signals detected under different penetration statuses. We can find that the AE signals mainly change in the time-domain in terms of amplitude fluctuation. It also consists of a lot of impulses, referring to the pulsed AE events. The cycle of each pulsed AE event is 0.10 s because the laser pulse frequency applied in our experiment is 10 Hz. When the amplitude distribution of a pulsed AE event in the time-domain is uniform, the corresponding welding process is relatively stable, which creates a better appearance of cross-section morphology. When the laser heat input is always less than the energy threshold for forming a penetration hole, the intensity of the acoustic signal intensity remains small (less than 2 V) under the non-penetration status.
As the higher laser energy is gradually delivered to the lower surface of the workpiece, the full-penetration hole may not necessarily be formed in each pulse period, which refers to critical-penetration status. Some pulse periods of acoustic signal intensity increased to exceed 2 V, and other pulse periods of signal intensity were less than 2 V. As the laser heat input further increased to significantly exceed the energy threshold, the full-penetration status easily occurred and the acoustic signal intensity reached a maximum value (larger than 3 V) under a stronger laser-beam illumination.
In general, the pulsed AE events can effectively reflect the keyhole dynamics of the pulsed laser affecting the workpiece during the laser welding process. The process-based definitions of various penetration modes come directly from time-resolved imaging of keyhole morphology. Taking a full-penetration weld as an example, Figure 3 displays the dynamic evolution of the keyhole and a single pulse AE event extracted from the complete AE signal during pulsed laser mode. Under the high instantaneous energy of the laser beam in Phase 1, a concave quickly appeared within a very short time due to a strong recoil momentum from the vaporization of the liquid surface. The intensity of the AE signal also began to increase. With an increase in laser energy in Phase 2, the laser fluence is sufficient to vaporize the metal, generating a large recoil pressure that opens a deep and slender vapor depression (referring to a gas-vapor keyhole) based on multiple reflections, which substantially enhances the laser energy absorption by the surrounding melt pool. Meanwhile, the amplitude of the AE signal increased significantly to a maximum value, which indicates that the formation process of the keyhole releases intense AE energy. From Phase 3 to 4, the laser absorption and reflection occurred due to its interaction with the vaporized metal; only part of the laser energy was absorbed by the molten metal. Thus, the depth of the keyhole decayed with time, and the amplitude of the AE signal began to decrease accordingly. After the pulse laser termination, the keyhole disappeared quickly over a very short interval of time and reappeared until the next laser pulse started. Meanwhile, the intensity of the AE signal also degraded rapidly due to the absence of laser energy, as illustrated in Phase 5.
By discussing the correlations between the transient nature of the subsurface keyhole and the AE signal, it is indicated that the acoustic emission could provide a wealth of information about the keyhole dynamics under the complicated interaction between the repetitively pulsed laser radiation and the aluminum alloy workpiece. To fully implement the acoustic monitoring system and manufacture defect-free welded products, we need a more comprehensive understanding of the acoustic source and its internal relationship to the weld quality.

2.3. Generation Mechanism of Acoustic Wave

As depicted in Figure 4a, a high-energy Gaussian-profiled laser beam (I ≥ 1 MW/cm2) focusing on the surface of a workpiece creates a thin and narrow capillary called a keyhole in a moving material. According to the hydrodynamic theory [30], the dynamic pressure balance at the keyhole wall is a premise of the stable keyhole that exists in the deep-penetration laser welding process. Significantly, the slight variation of the keyhole shape and size (depth d and exit radius R) could easily cause the pressure to fluctuate at the keyhole wall.
To keep the keyhole open and stable, the various driving pressures at the keyhole wall should satisfy a balance equation of pressure Δ p = 0 , as follows [30]:
Δ p = p γ ( r ) + p a b l ( T ) + p l + p g
Surface   tension   pressure :   p γ = γ R
Ablation   pressure :   P a b l ( T s ) = I r L v ( κ R T s M )
Hydrodynamic   pressure :   p l = 1 2 ρ v 2
Hydrostatic   pressure :   p g = ρ g d
where Pγ denotes the surface tension pressure induced by the molten metal around the keyhole, γ the surface tension coefficient, Pabl the ablation pressure induced by the material evaporation, Lv the latent heat of vaporization, Ir the laser intensity for material vaporization, κ the thermal diffusivity, Ts the keyhole surface temperature, M the molar mass of aluminium alloy, ρ the fluid density and g the gravity acceleration, v the flow velocity of the liquid metal fluid, and d the depth of the keyhole.
It is demonstrated that the equilibrium values of the keyhole sizes depend on the Gaussian laser energy distribution on the keyhole wall and material properties. With the increase in temperature at the bottom of the vapor depression, the laser absorption also increases abruptly, disrupting the force balance at the local vapor-liquid interface. The keyhole then becomes unstable and grows rapidly in depth. Due to the instability of the keyhole itself, any small deviation of the keyhole depth and radius from their equilibrium values leads to a pressure fluctuation and keyhole oscillation. Significantly, the response of keyholes to pressure fluctuations has been predicted to involve the excitation of a number of normal modes at characteristic frequencies [31]. The eigenfrequency f0 of keyhole oscillation can be calculated by Equation (6):
f 0 2 = γ d 2 + 2 ρ g R d 3 3 γ R 2 ln 3 × δ ( 3 R 5 + d 2 R 2 )
It can be seen that the eigenfrequency f0 is closely related to the keyhole depth d and radius R under the dynamic balance of the keyhole. Meanwhile, the keyhole oscillations (e.g., radial, axial, and azimuthal) induced by the interaction between the different laser energy intensities and materials could generate shock acoustic waves in different ranges of amplitude/frequency, which could be reflected in the acoustic intensity and spectrum. In addition, the acoustic emission in the audible frequency range is also likely to be related to the vapor plume ejection from the keyhole outlet. The keyhole-induced vapor plume compresses the ambient atmosphere and generates a high air pressure. It will further produce an acoustic wave with a large overall amplitude, depending on the plume size, which will rely on the laser energy input. Particularly, the expanding plume dynamics are influenced by the keyhole fluctuations caused by distortion of the pressure balance. The vapor ejection and the plasma size are modulated by the keyhole oscillation under the dynamic pressure fluctuations. Thus, the oscillation of the keyhole itself is considered a potential point acoustic source, and the molten pool acts like a speaker diaphragm, generating acoustic waves in the surrounding atmosphere.
To verify the rationality of the acoustic source mechanism, the keyhole depth was extracted under different laser powers with a common imaging-processing algorithm according to our previous studies [32,33]. From a series of variation curves of keyhole depths in each laser pulse (see Figure 4b), we found that: (1) when the laser energy is lower (3 kW), the depth variation of keyhole increases from 0.7 mm to 1.5 mm. The average depth (see the solid red line) also changes gradually, which is close to 1.2 mm. Consequently, the acoustic signal intensity is lower under the none-penetration status, as depicted in Figure 2a. This is because the oscillation behavior of the keyhole is not intense under the relatively smooth changing pressure acting on the keyhole wall; (2) when the laser energy is further increased, the keyhole depth dynamically changes from 1.8 mm to 3.0 mm, and the acoustic signal intensity also increases significantly under the full-penetration status (see Figure 2c). It is indicated that the evolution process of a penetrated keyhole releases intense AE energy due to a more intense keyhole oscillation behavior. Based on the theoretical and experimental analysis, the AE characteristics are a strong function of various influencing factors, including the material properties, experimental conditions (laser power, pulse duration, and frequency), and acoustic waves induced by the keyhole oscillation.
For better correlating the acoustic wave characteristics with the various weld penetrations, revealing the amplitude and spectrum distribution of acoustic signals with time-frequency analysis is indeed crucial to enhance the potential adoption of the AE technique, which will be discussed in the next section.

2.4. Time-Frequency Analysis of AE Signals

Since the spectrum property of the AE signal is closely related to the keyhole oscillation and vapor plume expansion generated by the interaction between the pulsed laser beam and materials, the frequency-domain analysis is a more accurate and efficient analytical approach, providing additional interesting insights compared to the time-domain. Additionally, the frequency-domain features on the specific frequency band are not sensitive to the environment and the machine noises from other frequency bands. In our study, a Fast Fourier Transform (FFT) algorithm was applied to extract the frequency characteristics of the AE signal segments, as shown in Equation (7):
X ( k ) = n = 0 N 1 x ( n ) ( cos 2 π k n N j sin 2 π k n N ) ( k = 0 , 1 , 2 , , N 1 )
where X(k) is the calculated data after FFT processing, x(n) is the recorded acoustic signal, and N is the sample points (N = 20 k). To directly reflect the distribution of frequency-domain characteristics over time, the acoustic signals belonging to different penetration states were transformed into a 2-D spectrogram using the short-time Fourier transform (STFT). The corresponding “time-frequency-amplitude” information is investigated in combination with the keyhole behavior. By shifting the window function h(t) on the time axis, the acoustic signal was analyzed segment by segment to obtain a local spectrum of the acoustic signals. The STFT of the raw signal x(t) is defined as:
S T F T t , f = + x τ h τ t e j 2 π f τ d τ
Figure 5 depicts the formation and dynamic evolution of the keyhole and the corresponding depth variation, as well as the 2-D spectrograms under different weld penetrations. It is found that the time-frequency spectrograms transformed by the STFT algorithm can visually display an obvious pulse component describing the pulse laser action over time. Moreover, the spectrum intensity of a low-frequency band (less than 2 kHz) under various weld penetration is almost equal, while the spectrum intensity of a high-frequency band under PP and FP states are larger than in the NP state. It may be due to of the increased laser energy input; the keyhole (see Figure 5) and the metal vapor ejections are becoming more intense, exhibiting much stronger fluctuations (larger amplitude and frequency). While the NP state does not have a deep enough keyhole, the laser energy is relatively small, melting the metal more slowly, making the metal vapor and plasma ejection less intense and the keyhole depression slower, so the AE event occurs in the low-frequency range.
Therefore, the results of spectrum distribution further indicate that the keyhole and its surrounding liquid molten pool layer act as a frequency-selective amplifier for pressure fluctuations induced by changes in the interaction of pulsed laser radiation with the keyhole wall. The experimental data show a strong correlation between the acoustic emission spectrum at 1~10 kHz and keyhole behavior as well as weld penetration. These acoustic spectra undoubtedly reflect a variety of quasi-periodic phenomena that are characteristic of the laser-metal interaction during laser welding. It is worth noting that an individual frequency component at a given frequency has no clear physical significance, which is not sufficient to be used as an efficient acoustic feature.

3. VMD-Based Frequency Component Extraction

To investigate more about the feasibility of distinguishing between low-frequency and high-frequency components for better describing the weld penetration, an efficient frequency decomposing technique is crucial. In this research, we apply the VMD (Variational Mode Decomposition) method, which is an adaptive signal processing method that decomposes a none-stationary AE signal f into K modal component sub-signals uk with each component having a defined finite bandwidth and a central frequency ωk. The method assumes that any signal consists of a series of sub-signals with a specific center frequency and finite bandwidth, i.e., an IMF (intrinsic mode function) [34].
Based on the classical Wiener filter, the modal function is obtained by solving the variational problem to obtain the central frequencies and bandwidth limits and finding the effective components corresponding to each central frequency in the frequency domain. The variational problem for the constraint is:
min u k ω k k = 1 K | | t δ t + j π t u k t e j ω k t | | 2 2 s . t . k = 1 K u k = f
where: {uk} = {u1,…, uK}, {ωk} = {ω1,…, ωK} are all modes and their central frequencies, respectively; t is the partial derivative with respect to t, δ(t) is the Dirac distribution, and * denotes the convolution.
To solve Equation (9), a quadratic penalty factor α and a Lagrange multiplier operator λ can be introduced to convert Equation (9) into an unconstrained variational problem, the augmented Lagrangian function is expressed as in Equation (10):
L u k , ω k , λ = α k = 1 K | | t δ t + j π t u k t e i ω k t | | 2 2 + | | f t k = 1 K u k t | | 2 2 + λ t , f t k = 1 K u k t
By using the alternating directional multiplier (ADMM) iterative algorithm combined with the Parseval/Plancherel, Fourier isometric transform, the modal components, and central frequencies were all optimally obtained, and the saddle points of the extended Lagrange function were searched for. The expressions for uk, ωk, and λ after the alternating search iteration are as follows [34]:
u ^ k n + 1 ω = f ^ ω i k u ^ i ω + λ ^ ω 2 1 + 2 α ω ω k 2
The same method of updating the central frequency is understood to be:
ω k n + 1 = 0 ω u ^ k ω 2 d ω 0 u ^ k ω 2 d ω
where μ ^ k n + 1 ω   is equivalent to the Wiener filter of the current residual f ^ ω i k u ^ i ω ; ω k n + 1 is the center of gravity of the current eigenmode power spectrum.
When using VMD for signal processing, a K value that is too large or too small will result in spurious or missing components, so the parameter K needs to be optimized. Using the instantaneous frequency averaging method [35], the instantaneous frequency of the jth sample point is obtained as fij, assuming that the ith modal component of the original signal has M samples. The original signal was pre-decomposed using the VMD method, and the modal classification at K = 2~9 was Hilbert transformed to obtain the resolved signals. By applying the instantaneous frequency averaging method:
f i = 1 N j f i j
Where N is the number of instantaneous frequencies of a modal component, and the penalty factor α was set to 2000. The normalized instantaneous frequency averaging curves for VMD modal components 2 to 9 were analyzed, as depicted in Figure 5.
It can be seen that when the modal component increases to a certain value, the curve bends distinctly from the neighboring points, and the instantaneous frequency mean of the modal component at that point changes abruptly, indicating that the different sub-signals produce abrupt changes in frequency during the decomposition. If the curve does not bend significantly but forms a gentler horizontal or vertical line with the adjacent points, it indicates that the different sub-signals have generated frequency blending during the decomposition of the signal. According to the instantaneous frequency averaging transformation curves from K2 to K9 in Figure 6, the result of the VMD decomposition layer at K9 is a relatively smooth folded line. There is no sudden change in the frequency components of each mode, which means that the VMD method at K9 can better decompose the raw AE signals into different IMF components.
Taking a full-penetration AE signal as an example, the signal was decomposed into 9-layer IMF components using the VMD method, and the penalty factor α was set to 2000 based on the discrimination method of center frequency observation. Figure 7a illustrates the various reconstitution sub-signals of each layer of the IMF. Then the FFT-based spectral analysis was performed for each IMF layer to visualize the nine IMF components of the raw AE signal, as shown in Figure 7b.
It can be seen that the frequency distribution of the IMF sub-signal derived from the raw AE signal after 9-layer decomposition is more reasonable, which is consistent with the interpretation of the instantaneous frequency averaging method. To calculate whether the VMD method separates signals with high energy in the low-frequency band, the time domain energy calculation method is introduced:
E n x = 0 N | x t | 2 d t
N is the maximum number of time series; N is set to 800,000. The energy percentage of each sub-signal is calculated as shown in Figure 8. It indicates that the 9th layer IMF component is in the lowest frequency band with the highest energy, and the remaining 8 IMF components are in higher frequency bands with lower energy. The comparison result demonstrates the feasibility and veracity of the VMD technique in reflecting the acoustic energy distribution among the low-frequency and high-frequency components of AE signals.

4. The Proposed Spatio-Temporal CNN-LSTM Model

During the pulsed laser welding process, affected by the high transient characteristics of the laser beam, the AE signal presents a highly dynamic, complex, and changeable characteristic. In particular, the energy distribution of different components in AE signals over time would vary under different welding conditions based on VMD pre-processing. Therefore, a deep learning algorithm based on the fusion of spatio-temporal information was proposed to make full use of the time-dependent relationship and spatial characteristics of the time-frequency information in a time-varying AE signal. In this section, a 1-D CNN model was firstly applied to deeply explore the spatial feature differences between the low-frequency high-energy components and the high-frequency low-energy components. After extracting the spatial features of the AE signal, an LSTM model was connected to mine the temporal features of this time series. Then, a hybrid CNN-LSTM model was constructed to predict the laser welding penetration status.

4.1. CNN Framework

In this study, the 1-D CNN architecture was used, and the decomposed AE sub-signals using VMD were used as the model input. Due to the advantage that 1-D CNN only convolves the filter width without the filter height [36], it can effectively convolve the 9-layer sub-signals at each time-series sample to deeply explore the energy distribution feature between low-frequency and high-frequency components. As shown in Figure 9, the filter width and height were set to 9 and 1, and the width and height of each sequence vector were 9 and 480,000. The filter sliding step was 1 from the start position along the height direction sliding to the final position, for a total of 479,999 sliding times. The second sequence vector was obtained after one filter convolution operation with the width and height changing, and was calculated as follows:
H 2 = H 1 F H 1 S 1 + 1
W 2 = W 1 F W 1 S 1 + 1
Taking the height calculation as an example, H1 and H2 are the heights of the first and second sequence vectors, respectively; FH1 is the height of the first convolutional layer; and S1 is the sliding length. The width was also calculated in the same way. Two layers of 1-D CNN were used to extract the sequence vector features in turn. Nine filters were set in the first layer to extract the width features of the sequence vectors, resulting in sequence vectors of 9 × 480,000. The second layer was set with six filters to extract the width features of the sequence vectors, resulting in sequence vectors of 6 × 480,000. Each 1-D CNN layer was then connected with one RELU layer as an excitation function, and a dropout layer was used at the end of the two convolution layers to prevent over-fitting.

4.2. LSTM Framework

LSTM (Long Short Term Memory) has a more accurate information transfer mechanism and can automatically store and remove temporal state information compared to RNN (Recurrent Neural Network). In addition, it can effectively solve the gradient disappearance problem and dependency problems caused by long input sequences [37]. As shown in Figure 10, the basic structure of the LSTM model uses input gates, output gates, and forgetting gates to protect and control the state of the storage unit.
Firstly, the input Xt is used in conjunction with the previous output ht−1, in order to obtain the output ft of the forgetting gate; the output range is used to filter whether the output of the previous unit is remembered, or not. The input gate is then responsible for selectively remembering the current input information and outputting it, and together they update the state information in the long-term memory chain Ct. Finally, ht−1 and Xt are passed through the output gate and multiplied by the activation function to obtain the output ht for that cell. The flow of information in the LSTM can be described as the following equation:
f t = σ W f h t 1 + X t + b t
i t = σ ( W i h t 1 + X t + b i )
C t = t a n h W c h t 1 + X t + b C
C t = f t C t 1 + i t C t
o t = σ ( W o h t 1 + X t + b o )
h t = o t t a n h C t
where Wf, Wi, Wo, Wc are weighting coefficients and bt, bi, bc are bias coefficients.
After extracting the deep spatial features of the temporal vectors using two 1-D CNN layers, the time-series size of the temporal vector became 6 × 480,000. Then, a LSTM layer was connected to explore the temporal features of this time series. As shown in Figure 11, one pulse period of the pulsed laser is 10 ms, and the sampling frequency was set to 20 kHz, i.e., one pulse period covers 2000 sampling points. With one pulse cycle as one frame, there are 240 frames within 480,000 sequence points. Therefore, the number of hidden units in the LSTM layer was fixed to 240, and each hidden unit processed one pulse cycle, which greatly speeds up the efficiency of the LSTM layer.

4.3. Establishment of the CNN-LSTM Hybrid Model

After adaptively decomposing the raw AE signal into 9-layer sub-signals with the VMD technique, a two-layer 1-D CNN model fully extracts the spatial features characterizing the time-frequency energy distribution of each sub-signal, while the LSTM model extracts the temporal state information reflecting the keyhole dynamics. Thus, the combination of CNN-LSTM hybrid models can describe the complex variation trend of acoustic waves under different welding conditions and obtain more accurate penetration predicting results. After combining with CNN and LSTM layers, a dropout layer was connected to prevent over-fitting with a dropout value of 0.5, followed by a FC (fully connected) layer, a Softmax activation function layer, and a classification layer including three types of penetration states (NP, PP, and FP). Figure 12 depicts the basic architecture of the CNN-LSTM hybrid model, which will be trained in order to verify the model’s performance in the next section.

5. Results and Discussion

5.1. Performance Evaluation of the Constructed Model

In this subsection, the predicting accuracy of different penetration states was evaluated using the constructed deep learning model. Due to the large number of parameters involved in CNN-LSTM, many hyperparameters (e.g., number of CNN or LSTM layers, batch size, and CNN filter number) directly affect the model recognition accuracy. Based on a great deal of experimentation, the optimized hyperparameters are provided in Figure 13; the rest of the parameters are defaults. A total of 800,000 samples of data were available for each of the three penetration states. 60% of the sample data from each state was used as the training set, 20% as the validation set, and the remaining 20% as the testing set. All experiments were deployed based on the Matlab 2022a environment and validated, using a computer equipped with an Intel(R) core(TM) i7-10750H processor and a NVIDIA RTX2060 GPU. The performances of implemented methods were assessed and compared within terms of recognition, accuracy, and corresponding loss values.
According to Figure 13, although there is a decline in accuracy after 40 iterations, the proposed model has a higher overall recognition accuracy of 99.8% after 150 training iterations. Meanwhile, the training and validation of model loss converged to zero after completing iterating. In addition, the normalized confusion matrix was calculated to clearly display the classification and misclassification results of each penetration status under different feature representations. As depicted in Figure 14, only one sample was misclassified as an FP state in an NP state classification, and the recognition accuracy reached about 100%. In the PP state classification, 183 samples were misclassified as NP, and the recognition accuracy reached 99.9%. In the FP state classification, 539 samples were misclassified as PP state, and the recognition accuracy reached 99.7%. The CNN-LSTM hybrid model can obtain a very high classification accuracy of 99.8%.

5.2. Performance Comparison between Different Models

To further demonstrate the superiority of the CNN-LSTM model, the performance comparison with other models was verified by examining the test accuracy. The effects of different network structures on the model’s performance were also tested. To guarantee the reliability of the proposed approach, 5-fold cross-validation was applied, and the average of those five times was calculated as the final accuracy of each model. The standard deviation (STD) was also calculated for each model’s accuracy to judge the model’s stability. The performance statistics of the replicate experiments for each model were depicted and compared in Figure 15, the parameters of the different networks are detailed in Table 3, and the following conclusions can be drawn:
  • No. 1 and No. 3 network models only compared the effect of different hidden layer nodes in LSTM on the model performance when other parameters remained unchanged. As the number of hidden nodes increases from 120 to 240, the average recognition accuracy also increases from 96.43% to 98.50%, which indicates that a lager node number in the LSTM hidden layer has a higher classification performance;
  • Since the 1st layer of CNN was applied to extract the spatial features of 9-layer sub-signals through VMD, the filter number of the 1st layer was fixed at 9, thus No. 2–No. 5 network models compared the effect of the filter number of the 2nd CNN layer in turn. It is worth noting that the No. 4 model has a best classification performance of 99.85% when the filter number of 2nd layer is set to 6. Meanwhile, it has a minimum accuracy fluctuation when STD = 0.21, confirming that the optimal model is stable. Although other models (No. 2 and No. 5) can also reach a larger recognition accuracy exceeding 98%, there exists a strong accuracy fluctuation at 50 or 100 iterations, indicating that these models are not stable;
  • As shown in Figure 15b, the No. 4 and No. 6 compared the effect of the VMD-based frequency components on the model’s performance. It can be found that the input data of the No. 6 network was not decomposed using the VMD method, which led to a significant accuracy (67.83%) decrease compared to the No. 4 network with the VMD method;
  • Lastly, the No. 4 and No. 7 compared the effect of the CNN-based spatial feature extraction on the model performance after the VMD processing. The final recognition classification of No. 4 (99.85%) is higher than No. 7 (95.04%), and the standard deviation of No. 4 (0.21) is also lower compared to No. 7 (0.28). Meanwhile, the No. 4 structure achieves the global optimum in fewer iterations, indicating that the CNN layer contributes greatly to the penetration and prediction of results.
Table 3. Detailed hyperparameters and performance statistics for different networks.
Table 3. Detailed hyperparameters and performance statistics for different networks.
No.1234567
VMD methodYesYesYesYesYesNoYes
Filter number
of 1st CNN layer
99999--
Filter number
of 2nd CNN layer
31369--
LSTM hidden nodes120240240240240240240
Max Acc(%)97.4598.5098.8299.9999.9970.3296.25
Min Acc(%)96.3497.7898.1099.3197.3165.3394.16
Avg Acc(%)96.4398.1698.5099.8598.3767.8395.04
Standard deviation0.282.330.250.212.320.320.28
Figure 15. (a) Comparison results between No. 1~No. 5 network models. For instance, the network structure of 9-3-120 means it has 9 filters in the 1st layer, 3 filters in the 2nd layer of the CNN model, 120 hidden nodes in the LSTM layer, and so on. (b) Comparison results between No. 4, 6, and 7 network models. No. 6 and No. 7 networks reflect the model only having a LSTM layer without a CNN layer, namely Origin-LSTM. No. 6 network deals with the raw AE signals without VMD, while the No. 7 network deals with the decomposed AE signals with VMD.
Figure 15. (a) Comparison results between No. 1~No. 5 network models. For instance, the network structure of 9-3-120 means it has 9 filters in the 1st layer, 3 filters in the 2nd layer of the CNN model, 120 hidden nodes in the LSTM layer, and so on. (b) Comparison results between No. 4, 6, and 7 network models. No. 6 and No. 7 networks reflect the model only having a LSTM layer without a CNN layer, namely Origin-LSTM. No. 6 network deals with the raw AE signals without VMD, while the No. 7 network deals with the decomposed AE signals with VMD.
Materials 16 01614 g015
In conclusion, the combination of the detected AE signals, VMD processing, and CNN-LSTM approach (No. 4 network) can obviously contribute to achieving excellent performance for predicting the penetration status during pulsed laser welding. The reason can be explained: the VMD technique can efficiently decompose the non-stationary AE signal into low-frequency and high-frequency components. After that, the proposed approach integrates the excellent spatio-temporal feature extraction capabilities of CNN and LSTM for dealing with the extracted frequency components. Compared to traditional feature extraction and machine learning methods used in other studies, such as in [15,16], the present method does not require the reduction of noise collected by an acoustic emission sensor and exhibits better adaptive feature extraction properties as well as higher prediction accuracy and a lower standard deviation, which is undoubtedly a great improvement in the model. Thus, the hybrid model provides a more accurate dynamic modeling approach for monitoring welding penetration based on AE sensing while being able to adequately characterize the keyhole dynamics during pulsed laser welding.

6. Conclusions and Future Works

This work developed a novel acoustic in-process penetration monitoring method based on the acoustic wave characteristics of a CNN-LSTM hybrid model for aluminum alloys in pulsed laser welding. The following conclusions can be drawn:
  • Combining with the high-speed photography and AE measurements, the characteristics of the AE signal are closely related to the keyhole oscillation and various weld penetrations. Based on the proposed mechanism of an acoustic source, the keyhole oscillation under the dynamic pressure fluctuation is considered a potential point acoustic source, and the surrounding molten pool acts like a speaker diaphragm, generating the spherical acoustic waves propagating in the workpiece;
  • According to the STFT time-frequency analysis, the acoustic spectra undoubtedly reflect a variety of quasi-periodic phenomena that are characteristic of the laser-metal interaction during laser welding; then the proposed VMD technique adaptively decomposed the raw AE signal into nine distinct frequency components, which can precisely characterize the acoustic energy distribution among the low-frequency and high-frequency components, under different welding penetrations, and improved frequency domain identifiability;
  • Finally, a novel hybrid model combing CNN and LSTM was designed to deeply mine the spatial and temporal acoustic features from the extracted frequency components. Extensive experiments demonstrate that our proposed approach yields a remarkable classification performance with a test accuracy of 99.8% and a standard deviation of 0.21, which obtains the best recognition performance compared with other deep learning methods.
The current research focuses mainly on the penetration status classification for the fully penetrated butt-joint welds during pulsed laser welding. Aiming at non-penetrated welds (lap or fillet joints), our future study will apply the acoustic emission and CNN-LSTM techniques for penetration depth regression forecasting, which can further verify its reliability and broaden its industrial application domain. In addition, future works will also focus on the combination of acoustic emission with thermal and photoelectric sensors and give full play to the potential advantages among monitoring techniques. It will provide the possibility to detect internal macroscopic seam defects (cracks and pores) and deliver better performance prediction during highly dynamic laser welding.

Author Contributions

Z.L.: investigating, writing-original draft. D.W.: conceptualization, methodology, reviewing, and funding acquisition. P.Z.: data curation and modeling. X.Y.: software, validation. H.S.: experimental testing. X.C.: supervision and editing. Y.T.: writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Class III Peak Discipline of Shanghai Materials Science and Engineering (High-Energy Beam Intelligent Processing and Green Manufacturing) and in part by the State Key Laboratory of Advanced Welding and Joining under Grant AWJ-23M05.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw/processed data and modeling codes required to reproduce these findings cannot be shared at this time as the data also forms part of an ongoing study.

Conflicts of Interest

The authors declare that they have no known competing financial interest, or personal relationship, that could influence the work reported in this paper.

References

  1. Svenungsson, J.; Choquet, I.; Kaplan, A. Laser Welding Process—A Review of Keyhole Welding Modelling. Phys. Procedia 2015, 78, 182–191. [Google Scholar] [CrossRef] [Green Version]
  2. Jia, Z.; Zhang, P.; Yu, Z.; Shi, H.; Liu, H.; Wu, D.; Ye, X.; Wang, F.; Tian, Y. Effect of pulse shaping on solidification process and crack in 5083 aluminum alloy by pulsed laser welding. Opt. Laser Technol. 2021, 134, 106608. [Google Scholar] [CrossRef]
  3. Gao, X.; Chen, Y.; You, D.; Xiao, Z.; Chen, X. Detection of micro gap weld joint by using magneto-optical imaging and Kalman filtering compensated with RBF neural network. Mech. Syst. Signal Process. 2017, 84, 570–583. [Google Scholar] [CrossRef]
  4. Xiao, X.; Liu, X.; Cheng, M.; Song, L. Towards monitoring laser welding process via a coaxial pyrometer. J. Mater. Process. Technol. 2020, 277, 116409. [Google Scholar] [CrossRef]
  5. Zalakain-Azpiroz, A.; Rodríguez, N.; de la Yedra, A.G.; Piccini, J.; Angulo-Vinuesa, X. A calibration tool for weld penetration depth estimation based on dimensional and thermal sensor fusion. Int. J. Adv. Manuf. Technol. 2021, 119, 2145–2158. [Google Scholar] [CrossRef]
  6. Eriksson, I.; Powell, J.; Kaplan, A. Signal overlap in the monitoring of laser welding. Meas. Sci. Technol. 2010, 21, 105705. [Google Scholar] [CrossRef]
  7. You, D.; Gao, X.; Katayama, S. Monitoring of high-power laser welding using high-speed photographing and image processing. Mech. Syst. Signal Process. 2014, 49, 39–52. [Google Scholar] [CrossRef]
  8. Tang, Z.-j.; Liu, W.-w.; Zhang, N.; Wang, Y.-w.; Zhang, H.-c. Real–time prediction of penetration depths of laser surface melting based on coaxial visual monitoring. Opt. Lasers Eng. 2020, 128, 106034. [Google Scholar] [CrossRef]
  9. Zhao, C.; Parab, N.; Li, X.; Fezzaa, K.; Tan, W.; Rollett, A.; Sun, T. Critical instability at moving keyhole tip generates porosity in laser melting. Science 2020, 370, 1080–1086. [Google Scholar] [CrossRef]
  10. Ju, S.; Li, D.; Jia, J. Machine-learning-based methods for crack classification using acoustic emission technique. Mech. Syst. Signal Process. 2022, 178, 109253. [Google Scholar] [CrossRef]
  11. Gu, H.; Duley, W. Resonant acoustic emission during laser welding of metals. J. Phys. D-Appl. Phys. 1996, 29, 550–555. [Google Scholar] [CrossRef]
  12. Gu, H.; Duley, W. A statistical approach to acoustic monitoring of laser welding. J. Phys. D-Appl. Phys. 1996, 29, 555–560. [Google Scholar] [CrossRef]
  13. Huang, W.; Kovacevic, R. Feasibility study of using acoustic signals for online monitoring of the depth of weld in the laser welding of high-strength steels. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2009, 223, 343–361. [Google Scholar] [CrossRef]
  14. Huang, W.; Kovacevic, R. A neural network and multiple regression method for the characterization of the depth of weld penetration in laser welding based on acoustic signatures. J. Intell. Manuf. 2009, 22, 131–143. [Google Scholar] [CrossRef]
  15. Yusof, M.; Ishak, M.; Ghazali, M. Classification of weld penetration condition through synchrosqueezed-wavelet analysis of sound signal acquired from pulse mode laser welding process. J. Mater. Process. Technol. 2020, 279, 116559. [Google Scholar] [CrossRef]
  16. Yusof, M.; Ishak, M.; Ghazali, M. Weld depth estimation during pulse mode laser welding process by the analysis of the acquired sound using feature extraction analysis and artificial neural network. J. Manuf. Process. 2021, 63, 163–178. [Google Scholar] [CrossRef]
  17. Cai, Y.; Luo, Y.; Wang, X.; Yang, S.; Zhang, F.; Tang, F.; Peng, Y. Physical mechanism of laser-excited acoustic wave and its application in recognition of incomplete-penetration welding defect. Int. J. Adv. Manuf. Technol. 2022, 120, 6091–6105. [Google Scholar] [CrossRef]
  18. Wu, D.; Huang, Y.; Zhang, P.; Yu, Z.; Chen, H.; Chen, S. Visual-Acoustic Penetration Recognition in Variable Polarity Plasma Arc Welding Process Using Hybrid Deep Learning Approach. IEEE Access 2020, 8, 120417–120428. [Google Scholar] [CrossRef]
  19. Wu, D.; Zhang, P.; Yu, Z.; Gao, Y.; Zhang, H.; Chen, H.; Chen, S.; Tian, Y. Progress and perspectives of in-situ optical monitoring in laser beam welding: Sensing, characterization and modeling. J. Manuf. Process. 2022, 75, 767–791. [Google Scholar] [CrossRef]
  20. Xie, S.; Zhang, Y.; Chen, C.; Zhang, X. Identification of nonlinear hysteretic systems by artificial neural network. Mech. Syst. Signal Process. 2013, 34, 76–87. [Google Scholar] [CrossRef]
  21. Movsessian, A.; Cava, D.G.; Tcherniak, D. An artificial neural network methodology for damage detection: Demonstration on an operating wind turbine blade. Mech. Syst. Signal Process. 2021, 159, 107766. [Google Scholar] [CrossRef]
  22. Lee, S.; Ahn, S.; Park, C. Analysis of Acoustic Emission Signals During Laser Spot Welding of SS304 Stainless Steel. J. Mater. Eng. Perform. 2013, 23, 700–707. [Google Scholar] [CrossRef]
  23. Shevchik, S.; Le-Quang, T.; Meylan, B.; Farahani, F.; Olbinado, M.; Rack, A.; Masinelli, G.; Leinenbach, C.; Wasmer, K. Supervised deep learning for real-time quality monitoring of laser welding with X-ray radiographic guidance. Sci. Rep. 2020, 10, 3389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Pandiyan, V.; Drissi-Daoudi, R.; Shevchik, S.; Masinelli, G.; Le-Quang, T.; Logé, R.; Wasmer, K. Semi-supervised Monitoring of Laser powder bed fusion process based on acoustic emissions. Virtual Phys. Prototyp. 2021, 16, 481–497. [Google Scholar] [CrossRef]
  25. Zhang, Z.; Qin, R.; Li, G.; Du, Z.; Li, Z.; Lin, Y.; He, W. Deep learning-based monitoring of surface residual stress and efficient sensing of AE for laser shock peening. J. Mater. Process. Technol. 2022, 303, 117515. [Google Scholar] [CrossRef]
  26. Shi, J.; Peng, D.; Peng, Z.; Zhang, Z.; Goebel, K.; Wu, D. Planetary gearbox fault diagnosis using bidirectional-convolutional LSTM networks. Mech. Syst. Signal Process. 2022, 162, 107996. [Google Scholar] [CrossRef]
  27. Song, R.; Jiang, Q. Application of VMD Combined with CNN and LSTM in Motor Bearing Fault. In Proceedings of the 2021 IEEE 16th Conference on Industrial Electronics and Applications (ICIEA), Chengdu, China, 1–4 August 2021; pp. 1661–1666. [Google Scholar]
  28. GB/T 3190-2020; State Administration of Market Regulation and the Standardization Administration. China Non-Ferrous Metals Industry Association: Beijing, China, 2020.
  29. GB/T 3880.2-2012; General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China. Standardization Administration of the People’s Republic of China: Beijing, China, 2012.
  30. Klein, T.; Vicanek, M.; Kroos, J.; Decker, I.; Simon, G. Oscillations of the keyhole in penetration laser beam welding. J. Phys. D-Appl. Phys. 1994, 27, 2023. [Google Scholar] [CrossRef]
  31. Klein, T.; Vicanek, M.; Simon, G. Forced oscillations of the keyhole in penetration laser beam welding. J. Phys. D-Appl. Phys. 1996, 29, 322. [Google Scholar] [CrossRef]
  32. Wen, X.; Wu, D.; Zhang, P.; Liu, S.; Luo, Z.; Jia, Z.; Ye, X.; Shi, H.; Yu, Z. Influence mechanism of the keyhole behavior on penetration depth by in-situ monitoring in pulsed laser welding of aluminum alloy. Optik 2021, 246, 167812. [Google Scholar] [CrossRef]
  33. Liu, S.; Wu, D.; Luo, Z.; Zhang, P.; Ye, X.; Yu, Z. Measurement of pulsed laser welding penetration based on keyhole dynamics and deep learning approach. Measurement 2022, 199, 111579. [Google Scholar] [CrossRef]
  34. Dragomiretskiy, K.; Zosso, D. Variational Mode Decomposition. IEEE Trans. Signal Process. 2014, 62, 531–544. [Google Scholar] [CrossRef]
  35. Zan, T.; Pang, Z.; Wang, M. Early Fault Diagnosis Method of Rollong Bearing Based on VMD. J. Beijing Univ. Technol. 2019, 45, 103–110. [Google Scholar]
  36. Kiranyaz, S.; Ince, T.; Gabbouj, M. Real-Time Patient-Specific ECG Classification by 1-D Convolutional Neural Networks. IEEE Trans. Biomed. Eng. 2016, 63, 664–675. [Google Scholar] [CrossRef]
  37. Hochreiter, S.; Schmidhuber, J. Long short-term memory. Neural Comput. 1997, 9, 1735–1780. [Google Scholar] [CrossRef]
Figure 1. View of the pulsed laser welding experimental platform and multi-information acquisition system: (a) schematic diagram; (b) experimental configuration.
Figure 1. View of the pulsed laser welding experimental platform and multi-information acquisition system: (a) schematic diagram; (b) experimental configuration.
Materials 16 01614 g001
Figure 2. Weld cross-section and top morphology as well as AE signals under different weld penetrations.
Figure 2. Weld cross-section and top morphology as well as AE signals under different weld penetrations.
Materials 16 01614 g002
Figure 3. Relationship between the keyhole behavior and its AE signal under full-penetration state during the repetitively pulsed laser welding: (a) formation and collapse dynamics of a gas-vapor keyhole; (b) waveform of a single pulsed AE event; (c) waveform of a raw AE signal; (d) time-dependent distribution of pulsed laser power.
Figure 3. Relationship between the keyhole behavior and its AE signal under full-penetration state during the repetitively pulsed laser welding: (a) formation and collapse dynamics of a gas-vapor keyhole; (b) waveform of a single pulsed AE event; (c) waveform of a raw AE signal; (d) time-dependent distribution of pulsed laser power.
Materials 16 01614 g003
Figure 4. Generation mechanism of acoustic waves during the pulsed laser welding process: (a) the pressure balance inside the keyhole under the laser beam radiation; (b) dynamic variation of keyhole depth under different laser powers.
Figure 4. Generation mechanism of acoustic waves during the pulsed laser welding process: (a) the pressure balance inside the keyhole under the laser beam radiation; (b) dynamic variation of keyhole depth under different laser powers.
Materials 16 01614 g004
Figure 5. Comparison of time-frequency characteristics with a 2-D spectrogram and keyhole depth variation under different weld penetrations.
Figure 5. Comparison of time-frequency characteristics with a 2-D spectrogram and keyhole depth variation under different weld penetrations.
Materials 16 01614 g005
Figure 6. Standardized instantaneous frequency averaging curves of each modal from decomposition layer K2 to K9.
Figure 6. Standardized instantaneous frequency averaging curves of each modal from decomposition layer K2 to K9.
Materials 16 01614 g006
Figure 7. 9-layer IMF component decomposition using VMD and FFT-based spectral visualization. (a) shows the variational mode decomposition of the original signal; (b) shows the frequency distribution of each components according to the original signal.
Figure 7. 9-layer IMF component decomposition using VMD and FFT-based spectral visualization. (a) shows the variational mode decomposition of the original signal; (b) shows the frequency distribution of each components according to the original signal.
Materials 16 01614 g007
Figure 8. Energy percentage of each IMF sub-signal using the VMD method.
Figure 8. Energy percentage of each IMF sub-signal using the VMD method.
Materials 16 01614 g008
Figure 9. 1-D Convolutional Neural Network structure for deeply exploring the energy distribution feature between low-frequency and high-frequency components (Structure includes 1-D convolution kernel, i.e., filter, and sequence vectors. Width × Height of the filter and sequence vectors are 9 × 1 and 9 × 480,000 respectively).
Figure 9. 1-D Convolutional Neural Network structure for deeply exploring the energy distribution feature between low-frequency and high-frequency components (Structure includes 1-D convolution kernel, i.e., filter, and sequence vectors. Width × Height of the filter and sequence vectors are 9 × 1 and 9 × 480,000 respectively).
Materials 16 01614 g009
Figure 10. The detailed architecture of the LSTM model.
Figure 10. The detailed architecture of the LSTM model.
Materials 16 01614 g010
Figure 11. The time-series sample determination for the LSTM model.((A) shows the orininal signal of the ABAE sensor; (B) shows the amplified part of the red box in (A); (C) shows the power of the pulsed laser machine mode.).
Figure 11. The time-series sample determination for the LSTM model.((A) shows the orininal signal of the ABAE sensor; (B) shows the amplified part of the red box in (A); (C) shows the power of the pulsed laser machine mode.).
Materials 16 01614 g011
Figure 12. The basic architecture of the CNN-LSTM hybrid model.
Figure 12. The basic architecture of the CNN-LSTM hybrid model.
Materials 16 01614 g012
Figure 13. Optimized hyperparameters of the proposed CNN-LSTM model and the final training/validation recognition accuracy after 150 iterations.
Figure 13. Optimized hyperparameters of the proposed CNN-LSTM model and the final training/validation recognition accuracy after 150 iterations.
Materials 16 01614 g013
Figure 14. Confusion matrix of the VMD-CNN-LSTM model.
Figure 14. Confusion matrix of the VMD-CNN-LSTM model.
Materials 16 01614 g014
Table 1. Detailed welding parameters used in experiments and corresponding welding quality.
Table 1. Detailed welding parameters used in experiments and corresponding welding quality.
CaseLaser
Power
Pulse
Width
Pulse FrequencyWelding SpeedDefocusing
Distance
Penetration
Status
15 kW10 ms10 Hz20 mm/s0 mmFP
24 kW10 ms10 Hz20 mm/s0 mmPP
33 kW10 ms10 Hz20 mm/s0 mmNP
(Full-penetration: FP, Partial-penetration: PP; None-penetration: NP).
Table 2. The national standards for the chemical composition and the mechanical properties of 6061 alloys.
Table 2. The national standards for the chemical composition and the mechanical properties of 6061 alloys.
AlloySi (%)Fe (%)Cu (%)Mn (%)Mg (%)Cr (%)Zn (%)Ti (%)Tensile Strength
(MPa)
Yield Strength
(MPa)
Elongation after Fracture (%)
60610.4–0.80.70.15–0.40.150.8–1.200.04–0.350.250.15≥290≥2407
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luo, Z.; Wu, D.; Zhang, P.; Ye, X.; Shi, H.; Cai, X.; Tian, Y. Laser Welding Penetration Monitoring Based on Time-Frequency Characterization of Acoustic Emission and CNN-LSTM Hybrid Network. Materials 2023, 16, 1614. https://doi.org/10.3390/ma16041614

AMA Style

Luo Z, Wu D, Zhang P, Ye X, Shi H, Cai X, Tian Y. Laser Welding Penetration Monitoring Based on Time-Frequency Characterization of Acoustic Emission and CNN-LSTM Hybrid Network. Materials. 2023; 16(4):1614. https://doi.org/10.3390/ma16041614

Chicago/Turabian Style

Luo, Zhongyi, Di Wu, Peilei Zhang, Xin Ye, Haichuan Shi, Xiaoyu Cai, and Yingtao Tian. 2023. "Laser Welding Penetration Monitoring Based on Time-Frequency Characterization of Acoustic Emission and CNN-LSTM Hybrid Network" Materials 16, no. 4: 1614. https://doi.org/10.3390/ma16041614

APA Style

Luo, Z., Wu, D., Zhang, P., Ye, X., Shi, H., Cai, X., & Tian, Y. (2023). Laser Welding Penetration Monitoring Based on Time-Frequency Characterization of Acoustic Emission and CNN-LSTM Hybrid Network. Materials, 16(4), 1614. https://doi.org/10.3390/ma16041614

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop