This study investigates the innovative application of machine learning (ML) models to predict critical parameters—stress magnitude (SM), peak temperature (PT), and plastic strain (PS)—in ultrasonic welding of metallic multilayers. Extensive numerical simulations were employed to develop and evaluate three ML models: Radial Basis Function (RBF), Random Forest (RF), and Kernel Ridge Regression (KRR). According to the results, the KRR model demonstrated superior performance, achieving the lowest RMSE and highest R² values of 0.068 (R² = 0.941) for SM, 0.075 (R² = 0.929) for PT, and 0.071 (R² = 0.946) for PS, with fewer data samples required. KRR also exhibited low squared bias and variance values, ranging from 1×〖10〗^(-4)-3.2×〖10〗^(-4) for bias and 2.2×〖10〗^(-4)-3.6×〖10〗^(-4) for variance, indicating its precision in predicting the output targets. Moreover, the systematic categorization of input features, including material properties, geometrical factors, and welding parameters, highlighted their significant influence on predictive accuracy, particularly the crucial role of welding parameters at higher output values. Finally, a case study on ultrasonic welding of copper multilayers underscores the model's effectiveness in unraveling complex relationships, providing a robust tool for optimizing and advancing ultrasonic welding processes.