Between 7000 and 12,000 spot welds hold a car body together depending on its size [
1] and 91.7 million cars were produced in the automotive industry worldwide [
2]. One major technological and economical advantage of resistance spot welding over other joining technologies is the ease of integration in automated production lines. Furthermore, it is a lucrative process in technological and organizational terms. The process setup is simple, and the process cycle is, in the order of milliseconds, very short. After stacking one metal sheet on the top of another, pneumatic driven cylinders move electrode caps to clamp the metal sheet stack. Subsequently, the metal sheets are pressed together and a load of up to several hundred MPa is applied (squeeze time). Then, the electrodes are connected to a voltage, and thus an electric current crosses the metal sheets, in which Joule heat initiates weld nugget growth and sheet metal fusion (welding time). At the end, the electrodes rest briefly upon the metal sheets (hold time) before they are moved to the next weld spot (off time) and the process cycle restarts. The most important quality criterion for spot welds is its tensile strength. It can be determined by destructive test methods such as chisel or tensile test. However, after the test the weld is destroyed, and further deployment is impossible. Alternatively, non-destructive test methods can evaluate the weld quality and warrant further use of the assembly after testing. A widespread non-destructive method is ultrasonic testing which aims at detecting the effective contact area size of the weld joint. If this area exceeds a minimal threshold, the joint weld is accepted and dismissed otherwise. The equipment is expensive, requires qualified staff, and ultrasonic testing results scatter broadly. It tends to underestimate the welding spot diameter by approximately up to 2.5 mm [
3]. Other non-destructive methods are numerical methods that simulate nugget growth and geometric parameters, such as nugget diameter and penetration depth, which—due to their correlation with weld strength—allow an indirect assessment of the weld joint quality. However, current RSW models are too slow for integration into real-time monitoring and control systems, as discussed in the state-of-the-art section. An online monitoring system for RSW bears the potential of adjusting welding process parameters from one manufactured spot weld to the next in the assembly line, and, as a consequence, reducing the number of NOK welds and save time, costs, and energy. In the early stages, refs [
4,
5] applied the finite-difference method (FDM) to electro-thermal models of RSW. With the assumption of a constant electric current, ref [
4] proposed a model to predict the temperature distribution as a function of time and space, allowing for variations in the mechanical properties of the sheet metal. Ref [
5] presents the temperature-dependent electrical potential distribution in the base metal and the interfaces for various electrical currents. In both studies, the nugget diameter and penetration depth derived from the computed temperature field showed good agreement with experimental data. Ref [
6] used a control volume formulation and central differences to model the dimensionless temperature field, the nugget growth for different welding currents, electrode tip shapes and thickness ratios of work pieces. The enthalpy-temperature relation was capitalized to account for the phase change. The simulation results in terms of weld nugget growth, nugget thickness and shape were consistent with experimental results. In [
7] the finite volume method (FVM) was adopted to build a complex RSW simulation model, which considered—among other aspects—the electric current density, the magnetic field intensity, the temperature, and the velocity field for work pieces with flat faced or truncated electrodes. The effects of the electrode face radius and cone angle on transport mechanism, for example, mass transfer, and various other non-linear phenomena were clearly demonstrated; simulation results agreed well with experimental data. In [
8], the mass, momentum, heat and species transport, as well as the magnetic field intensity, were discretized by a control-volume formulation to compute the dynamic electrical resistance during RSW. The simulation result suggest that the dynamic resistance of AISI 1008 steel can be divided into four distinct stages, in which the contact resistance and the bulk resistance contributions vary over time. Several years later, ref [
9] developed a control volume based finite difference code for the electrical and thermal field and combined it with a commercial code that provided the mechanical model. Based on this hybrid-approach, the computed nugget size deviated from the experimental data by merely three percent. Many simulations of RSW are based on the finite-element method (FEM), which derives model equations from integration over the finite-element domain. For example, the general-purpose simulation program ABAQUS© (Version 5.7) was used to conduct a parametric study on different electrode shapes, welding currents, and electrode forces for Al-alloys in [
1]. Ref [
10] analyzes the influence of electrode-water cooling on welding of aluminum alloys AA5182. It utilizes LS-DYNA© (R11.0) to build a thermal-electrical-mechanical model. The simulation results indicated that water cooling affects the temperature distribution in the sheets only slightly, and thus, does not influence nugget growth at all. However, it has a significant effect on the electrode cooling during hold-time. Another study on aluminum alloys for RSW processes was conducted in [
11], where a calibrated contact resistance model for AA5182 was presented. The underlying electric–thermal–mechanical FE model could reproduce weld nugget diameters deviating from real experiments by four percent. FEM-based, SORPAS© 3D is a special purpose simulation program with a multi-physics model for RSW [
12]. It was applied to investigate short-pulse welding on aluminum alloy AA6016-T4 to reduce the required energy for producing sound welds in thin sheets [
13]. The nugget formation was found to happen in two distinct phases: the nucleation, in which 60–80% of the final diameter evolves due to high contact resistance, and the growth stage, when further nugget growth is induced by heat conduction. Savings of approximately 50% regarding energy and time were achieved. SORPAS©.2D provides simulation results with high accuracy, but it requires approximately an hour to run an RSW-simulation with a resolution of 1000 finite elements on a conventional desktop computer [
14].
In view of the development labor in past decades, it becomes clear that in the framework of appropriately set model assumptions and on the fundament of suitable material and process data, the nugget growth of RSW can be simulated with sufficient accuracy. Hitherto, numerical analyses of RSW in academia and commercial special purpose programs paid attention to simulation accuracy rather than computational speed; none of the previously cited research papers indicated information on the simulation run-time. In terms of a model-based real-time monitoring and quality assessment system, two essential requirements can be formulated: it should be capable of differentiating between OK and NOK spot welds (sufficient model accuracy) and the time window for the quality assessment is to be shorter than the time between two consecutive spot welds (computation speed). In numerical simulation, these requirements contradict each other, that is, increased simulation accuracy comes along with increased computational time. In the sense of this optimization problem, it is reasonable to constrain the simulation model to physical phenomena, which are predominantly relevant to RSW, to reduce the computation cost. RSW is based on Joule heating, i.e., the transformation of electrical to thermal energy, to join metal stacks by fusion. Thus, the electric-thermal model is considered as the core of any multi-physics model for RSW. The finite difference method, for heat transfer in solids elaborately described in [
15], derives model equations from replacing partial derivative terms by finite differences. In this paper, this numerical method is applied to develop an electric-thermal model for RSW by means of the Peaceman–Rachford scheme. It leads to a set of linear equations with tridiagonal band matrices which are solved by the Thomas Algorithm rapidly. Based upon this solving algorithm, the electric-thermal model is verified and investigated on its suitability for real-time simulations.