Table 2.1 Thermal properties of some common materials.
Material | Thermal diffusivity α (m2 s−1) | Volume thermal capacity ρC s (Jm−3 K−1) | Thermal conductivity k (Jm−1s−1K−1) | Melting point (K) |
Aluminum | 8.5 × 10−5 | 2.7 × 106 | 229.0 | 933 |
Carbon steel | 9.1 × 10−6 | 4.5 × 106 | 41.0 | 1800 |
Austenitic stainless steel | 5.3 × 10−6 | 4.7 × 106 | 24.9 | 1773 |
Ti alloy | 9.0 × 10−6 | 3.0 × 106 | 27.0 | 1923 |
Copper | 9.6 × 10−5 | 4.0 × 106 | 384.0 | 1336 |
2.2.2.2 Rosenthal's Three‐Dimensional Equation
The analytical solution derived by Rosenthal for three‐dimensional heat flow in a semi‐infinite workpiece during welding, Figure 2.18, is as follows [25]:
(2.9)
where R is the radial distance from the origin, namely, (x 2 + y 2 + z 2 ) 1/ 2. Equation (2.9) can be used to calculate the steady‐state temperature T(x, y, z), with respect to the moving heat source, at any location in the workpiece (x, y, z), for instance, at x = 1 cm, y = 4 cm, and z = 0 cm, as shown in Figure 2.18. The temperatures at other locations along y = 4 cm can also be calculated, and the temperature distribution along y = 4 cm can thus be determined.
Figure 2.18 Three‐dimensional heat flow during welding of semi‐infinite workpiece.
2.2.2.3 Step‐by‐Step Application of Rosenthal's Equations
The procedure for calculating thermal cycles is illustrated as follows. Equation (2.9) for 3D heat flow is used, with V = 2.4 mm/s, Q = 3200 W and T o = 25 °C. Consider the temperature distribution along the straight line of y = 4 cm.
Step 1. Calculate T (1, 4). Let x = 1 cm, y = 4 cm, R = [(1)2 + 42 + 02]1/2.
Step 2. Calculate T from Eq. (2.9) using k and α of the workpiece material (e.g. steel).
Step 3. Repeat Steps 1 and 2 for x = −1, 0, 1, 2, 3, 4, 5 cm.
Step 4. Sketch the temperature distribution T (x, 4, 0).
Step 5. Convert temperature distribution T (x, 4, 0) into thermal cycle T (t) by dividing x by V.
Figure 2.19 explains how to convert the calculated temperature distribution in Figure 2.18 to a thermal cycle. It is assumed that the workpiece in Figure 2.18 is long enough such that, with respect to the moving coordinate system, the temperature field does not change. The T – x plot, i.e. the temperature distribution along the welding direction, can be converted into T – t plot, namely, the thermal cycle, by calculating time t using t = (x − 0)/V. For instance, assume the travel speed of the heat source is V = 4 mm/s. At the point x = 2 cm, y = 4 cm, and z = 0 cm, t = (20 mm − 0)/(4 mm/s) = 5 s.
Figure 2.19 Converting the calculated temperature distribution in Figure 2.18 to thermal cycle, assuming the temperature field is steady with respect to the moving coordinate system.
The shape of the weld pool can be calculated such as at the top surface of the weld pool, i.e. at z = 0 cm. Let T L be the liquidus temperature of the workpiece material, e.g. 1530 °C for steel. The liquidus temperature of an alloy is equivalent to the melting point of pure metal T m, above which the workpiece material is melted completely. The procedure for calculating pool shape is illustrated using Eq. (2.9) for 3D heat flow and V = 2.4 mm/s, Q = 3200 W and T o = 25 °C.
Step 1. Let T = T L in Eq. (2.9). Calculate R at x = 0 cm and z = 0 cm from Eq. (2.9) by substituting into the equation the values of T o , V, Q, and the physical properties of the workpiece materials (e.g. steel) k and α.
Step 2. Find the value of y from R = (02 + y 2 + 02)1/2.
Step 3. Repeat steps 1 and 2 for x = 1 cm, 2 cm, etc. Use the calculated y values to construct the pool shape at the workpiece surface.
Let T H be the temperature at which solid‐state phase transformation occurs (e.g. 780 °C). The shape of the HAZ can be calculated, such as at the top surface of the weld pool, i.e. at z = 0 cm, following the same procedure for calculating the pool shape except that